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Abstract 



After a review of the problems associated with the conventional treatment of parti- 
cle oscillations, an oscillation formula is derived within the framework of quantum 
field theory. The oscillating particle is represented by its propagator and the ini- 
tial and final states by wave packets. It is obviously relativistic from the start and 
moreover applies both to stable (neutrinos) and unstable particles (K and B mesons, 
unstable neutrinos). CPLEAR and DAFNE experiments are studied as examples, 
with special attention directed to CP violation. The problems resulting from equal 
energies/momentum/velocities prescriptions are analyzed and solved. Oscillations of 
associated particles are found to be nonexistent. The relativistic generalization of 
the Wigner-Weisskopf equation is also derived. 
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Introduction 



Bien que la theorie des champs serve principalement au calcul de processus 
physiques microscopiques, elle permet aussi l'etude de phenomenes macroscopiques. 
Par macroscopique j'entends une distance ou un temps tels que leur inverse est 
beaucoup plus grand en unites naturelles (h = 1, c = 1) que l'energie typique 
impliquee dans le processus. Le phenomene le plus fondamental de ce type est la 
propagation d'une particule entre deux points de l'espace-temps, separes par une 
distance et/ou un temps macroscopiques. 

En physique classique, la probabilite de detecter une particule stable est cons- 
tante au cours du temps, par definition meme de la stabilite, tandis que la pro- 
babilite de detecter une particule instable decroit au cours du temps comme une 
exponentielle : V(t) — e~*/ r , ou r est le temps de vie moyen. Ce comportement 
decouvert par Rutherford est derive en physique classique bien que la disintegration 
soit un phenomene quantique. En physique quantique, cette image fonctionnant 
merveilleusement bien du point de vue experimental est modifiee sous deux aspects. 

En premier lieu, la brique de base de la physique quantique n'est plus la pro- 
babilite de detection mais l'amplitude de probabilite, dont la norme au carre est 
en correspondance avec la probabilite classique de detection. Cette amplitude suit 
une evolution deterministe. La mecanique quantique nous dit que l'amplitude de 
detection d'une particule stable de masse m oscille dans le repere au repos de la 
particule comme une exponentielle : A(t) ~ exp(— imt). La probabilite de detection 
d'une particule stable libre reste done constante au cours du temps. La realite de 
ce concept quantique se manifeste le plus clairement a travers 1' interference obser- 
vable entre plusieurs amplitudes, dont des exemples verifies experiment alement sont 
1' experience a deux trous pour des electrons, l'effet Aharonov-Bohm et les oscilla- 
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tions dans le systeme K°K°. 

En deuxieme lieu, les principes fondamentaux de la theorie quantique interdisent 
une forme exactement exponentielle pour la loi de disintegration. La raison en est 
que les produits de disintegration peuvent se recombiner pour reformer l'etat initial 
HI . On peut le voir en divisant l'intervalle de temps £ en deux sous-intervalles t\ et 
£2. L'amplitude de probabilite doit satisfaire a l'equation 

A(h + £ 2 ) = A(h) A(t 2 )+ < % nst \e~ iHt2 1 Vdesih) > , 

ou ^mst symbolise l'etat instable et ^des{ti) represente l'etat des produits de disinte- 
gration au temps t\. En ne gardant que le premier terme du membre de droite 
de cette equation, on obtient la formule classique de disintegration exponentielle. 
Le second terme, purement quantique, modifie ce comportement : les produits de 
disintigration peuvent retourner a l'itat initialQ. Bien que l'existence thiorique de 
ces corrections a la loi de Rutherford soit connue depuis longtemps, elle n'a pas pu 
etre testie expirimentalement jusqu'a prisent 0. L'ordre de grandeur de ces correc- 
tions peut etre estimi en micanique quantique, mais il vaut mieux recourir tout de 
suite a la thiorie des champs, puisqu'elle constitue le seul cadre thiorique cohirent 
pour la description des particules instables. Dans ce formalisme, les contributions 
exponentielles sont en correspondance avec les poles complexes du propagateur de la 
particule |6|, §], et l'amplitude de propagation macroscopique oscille dans le repere 
au repos de la particule comme une exponentielle : A(t) ~ exp(— imt — Tt/2). La 
partie imaginaire de l'exposant s'interprete comme le produit de la masse de la par- 
ticule et du temps tandis que la partie rielle de l'exposant s'interprete comme le 
produit de l'inverse du temps de vie moyen par le temps. Les contributions non 
exponentielles proviennent de l'existence de seuils de production d'itats a plusieurs 
particules ainsi que par des facteurs de seuil figurant dans les amplitudes de pro- 
duction et de ditection de la particule se propageant f|, [5]. On peut virifier dans ce 
cadre que ces corrections sont nigligeables pour toute propagation macroscopique. 
Nous y reviendrons. 

1 D'autres demonstrations de l'existence de corrections non exponentielles existent, tant en 
mecanique quantique ||, ||, qu'en theorie des champs [EL |l ||. 
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Dans ce travail, nous allons nous interesser principalement a V oscillation de par- 
ticules. Si les particules en question sont quasiment degenerees en masse et melangees 
par une interaction, des effets d'interference apparaitront sous la forme d'oscillations 
spatiales de la probabilite de detection des particules. Ces effets dependent directe- 
ment de la difference de masse entre les etats se propageant. Cette difference doit 
etre tres petite pour qu'une oscillation spatiale macroscopique soit observable. II 
faut done verifier que ce qui a pu etre neglige a l'echelle de la masse dans l'analyse 
de la propagation d'une seule particule peut encore etre neglige a l'echelle de la 
difference de masse entre des etats quasiment degeneres. Dans le cas du systeme 
K°K°, la difference de masse causant l'oscillation est du meme ordre que les cor- 
rections habituellement negligees en T/m : Am/m ~ T s /m ~ I0~ 14 ! D'autres 
problemes se posent a propos de l'interface classique/quantique. En effet, bien que 
les oscillations dans la probabilite de detection soient de nature purement quantique, 
leur observation n'est possible que si les particules sont bien localisees dans l'espace, 
comme des particules classiques. Nous developperons ces sujets ulterieurement. 

Pourquoi les oscillations de particules suscitent-elles tant d'attention? Deux 
motivations principales sous-tendent cet interet : d'une part, l'etude de la violation 
CP dans les systemes oscillants K°K° et B°B°, et d'autre part la mesure des masses 
des neutrinos au moyen de leurs oscillations. 



Bien que la prediction || et l'observation |TU] des oscillations des kaons datent 
de plus de quarante ans, le phenomene voit son interet sans cesse renouvele. La 
premiere raison est que le systeme des kaons neutres fut jusqu'a l'annee passee le 



seul endroit ou une violation CP se manifestait. L'observation |TT] de cette vio- 
lation dans les oscillations des kaons, dite violation CP indirecte, est maintenant 
doublee par l'observation [O] d'une violation CP dans leurs disintegrations, dite 
violation CP directe. Si la violation indirecte est vue comme une perturbation d'ordre 
e d'un systeme K°K° respectant CP, les experiences montrent que e ~ I0~ 3 et 
que la violation CP directe est de l'ordre de e 2 . II est done important de disposer 
d'une description coherente du systeme a l'ordre 0(e 2 ). Ce n'est pourtant pas le 
cas du formalisme couramment utilise dans la litterature, qui est 1' approximation 
de Wigner-Weisskopf [f[3| . Dans cette approche, revolution temporelle des etats 
instables obeit a une equation effective de type Schrodinger dont 1'hamiltonien non 
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hermitien permet les disintegrations. Par consequent, les matrices de diagonalisation 
de l'hamiltonien sont en general non unitaires, les etats correspondants ne sont pas 
orthogonaux et leur normalisation ne peut se faire sans ambiguite. Ce probleme est 
souvent neglige, car, dans la limite ou la symetrie CP est respect ee, la matrice de 
diagonalisation est unitaireQ. II en result e que la non-orthogonalite des etats propres 
de masse est de l'ordre de e et que l'ambiguite de normalisation est encore plus petite, 
de l'ordre de e 2 . C'est pourquoi elle est generalement negligee pour la description de 
la violation CP indirecte. On a vu ci-dessus qu'il ne peut plus en etre question dans 
l'etude de la violation CP directe qui est du meme ordre de grandeur. 

Une seconde raison pour s'interesser aux oscillations de kaons est la transposi- 
tion immediate du formalisme approprie aux kaons a d'autres systemes meson/anti- 
meson. Le plus prometteur est le B°B°, dont les oscillations ne sont observees que 
depuis 1987 [tL5) tandis qu'une premiere mesure de la violation CP dans ce systeme 



a ete effectuee recemment [16|. Bien que le formalisme de Wigner-Weisskopf n'ait 
pas ete jusqu'a present en contradiction avec les donnees experimentales, il n'est pas 
certain que cette description possede le meme degre de validite pour modeliser la 
violation CP du systeme B°B°, ou les corrections relativistes sont beaucoup plus 
importantes. 

Une derniere raison d'etudier ces systemes de meson/antimeson est l'etude 
experimentale de correlations quantiques macroscopiques (effet EPR) de deux kaons 



ou mesons B oscillant simultanement [17 



L'autre motivation principale de s'interesser aux oscillations spatio-temporelles 
de la probabilite de detection concerne les neutrinos. Historiquement, les premieres 
donnees experimentales suggerant des neutrinos massifs proviennent de l'anomalie 
du flux des neutrinos solaires, qui est bien inferieur aux predictions des modeles so- 
laires les plus sophistiques |18, [19|, [2(|. En 1968, Bruno Pontecorvo attribue le deficit 



2 Les symetries CPT et CP imposent toutes les deux l'egalite des elements diagonaux mais 
la symetrie CP impose en plus l'egalite des elements non diagonaux de l'hamiltonien |C3 ; ces 
contraintes permettent une diagonalisation d'une matrice non hermitienne par une matrice unitaire. 
En toute generalite, la condition necessaire et suffisante pour diagonaliser une matrice par une 
transformation unitaire est que la matrice soit normale, c'est-a-dire qu'elle commute avec son 
hermitienne conjuguee. 
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en neutrinos electroniques a des oscillations entre neutrinos de differents nombres 
leptoniques plfl . Plus tard, des solutions recourant a des neutrinos instables sont 



aussi proposees p2| . Plus recemment, on assiste a un retour en force du concept 
de neutrinos massifs en raison de l'observation d'une anomalie angulaire dans le 
flux des neutrinos atmospheriques P3fl . Les deux modeles les plus simples resolvent 
ce probleme en attribuant une masse aux neutrinos. Le premier modele explique 



l'anomalie par une oscillation de neutrinos [24 tandis que le deuxieme recourt a un 



neutrino instable [[25j . Une derniere raison pour soupgonner une masse aux neutrinos 
provient de l'experience LSND | 2"6|| . Quelle que soit l'explication de ces differentes 
anomalies, il parait difficile d'eviter le recours a des neutrinos massifs, ce qui n'est de 
toute fagon interdit par aucun principe fondamental. Les anomalies cities impliquent 
des masses pour les neutrinos extremement faibles par rapport aux energies typiques 
des processus concernes et done des oscillations macroscopiques sur de tres grandes 
distances. Des lors, il s'avere crucial de disposer d'une description coherente de la 
propagation de neutrinos massifs, puisque les effets d'oscillations resultent de la 
compensation presque totale d' energies- impulsions des differents etats oscillants, en 
raison de la petitesse des differences de masses par rapport a l'energie de la parti- 
cule. Notons enfin que la masse des neutrinos rend possible une violation CP dans 
le secteur leptonique. Etant donne que l'explication de l'anomalie des neutrinos at- 
mospheriques par une disintegration reste toujours possible, il me parait interessant 
d'etablir un formalisme oil les deux phenomenes d'oscillation et de disintegration 
coexistent. Cet objectif nous ramene au systeme K°K°, ou la disintegration et 
l'oscillation vont de pair. Dans ce travail, les analyses des propagations des kaons et 
des neutrinos seront regroupees dans un seul formalisme derive a partir des principes 
de la theorie des champs. 

Dans le premier chapitre, le concept d'oscillation est defini en mecanique quan- 
tique et le traitement traditionnel de ces oscillations est presente. Les difhcultes 
soulevees par cette approche sont passees en revue, avec pour conclusion que seul 
un calcul en theorie des champs peut clarifier ces problemes. Dans le deuxieme 
chapitre, un modele simplifie de propagation d'une particule est etabli. II nous con- 
duit a l'etude de la transformee de Fourier du propagateur, a l'obtention de la tradi- 
tionnelle loi d'evolution en exponentielle ainsi que des corrections non exponentielles 
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negligeables pour des propagations macroscopiques. Le troisieme chapitre est con- 
sacre a l'application de ce formalisme a la propagation de particules en melange : 
apres diagonalisation de leur propagateur, la probabilite de propagation et d'oscil- 
lation de ces particules est calculee dans le modele simplifie du deuxieme chapitre. 
On constate alors que ce modele ne suffit pas a repondre a toutes les interrogations 
initiates, bien qu'il eclaircisse deja une serie de points delicats. Dans le chapitre sui- 
vant, un modele de propagation plus realiste est developpe, toujours dans le cadre 
de la theorie des champs. Les etats entrants et sortants y sont modelises par des 
paquets d'ondes. Le calcul de la probabilite de detection d'une particule instable 
fournit la reponse attendue. Les corrections non exponentielles sont aussi recalculees 
et restent negligeables. II ne reste plus qu'a appliquer ce modele aux oscillations de 
particules, ce qui est l'objet du cinquieme chapitre. Avant de deriver la probabilite 
d'oscillation, on verifie que les corrections non exponentielles propres aux melanges 
de particules sont negligeables. Une probabilite d'oscillation dependant uniquement 
de la distance est ensuite obtenue et analysee terme par terme. Les differents fac- 
teurs y figurant representent la decroissance exponentielle due a l'instabilite de la 
particule, l'oscillation due au melange, la disparition de l'oscillation a grande dis- 
tance due a la decoherence ainsi que les conditions d'observabilite des oscillations 
dues aux conditions de production et de detection. Le chapitre se termine par une 
comparaison de nos resultats avec les differentes approches existantes. Au sixieme 
chapitre, on etablit des prescriptions de calcul simplifiant l'application de nos for- 
mules a des cas concrets. Ensuite, certaines quantites observables des experiences 
CPLEAR et DA<3>NE sont calculees dans notre modele. Enfin, le travail se termine 
par une analyse phenomenologique rigoureuse de la violation CP dans le systeme 
meson/antimeson. 

La nouveaute de ce travaille reside en l'application de la theorie des champs a 
la description du melange de particules instables, dans le but de traiter de facon 
unifiee les cas des mesons K et B et des neutrinos, avec l'espoir d'en tirer des 
legons interessantes pour chaque cas particulier. Comme nouveau resultat, nous 
obtenons une formule relativiste et directement applicable aux experiences, puisque 
ne dependant que de la distance de propagation et de l'observation des etats initiaux 
et finaux. Le recours a des concepts classiques exterieurs au formalisme est evite. Le 
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traitement de la particule oscillante comme etat intermediate permet de montrer 
l'inanite des questions suivantes, infiniment ressassees dans de nombreux articles : 
quelle est l'energie-impulsion de la particule oscillante?, comment definir des etats 
physiques non orthogonaux?, comment definir un etat propre de saveur? (notre 
reponse a cette derniere question n'est pas originale). Plusieurs controverses concer- 
nant la longueur d'oscillation et l'oscillation eventuelle des particules associees a la 
production de la particule oscillante trouvent enfin une solution claire dans notre 
formalisme, qui fournit sans ambiguite la generalisation relativiste de la formule de 
Wigner-Weisskopf pour un melange de particules. 

Les conditions d'observabilite des oscillations (necessite d'une incertitude sur 
l'impulsion de la particule oscillante, existence d'une longueur de coherence) sont 
etudiees pour la premiere fois dans le cas des particules instables. Comme sous- 
produit de notre analyse, nous avons obtenu une estimation, nouvelle en theorie des 
champs, des corrections non exponentielles a la propagation d'un melange de parti- 
cules. Comme nouvelles applications, nous avons traite rigoureusement la phenome- 
nologie de la violation CP a partir de nos formules. Le cas d'une oscillation double 
(ffi-^KK) est aussi etudie pour la premiere fois en theorie des champs. 
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Chapitre 1 

Oscillations en mecanique 
quantique 



1.1 Oscillations temporelles en mecanique 
quantique 

L 'etude de l'oscillation de la probability de detection pour un melange de par- 
ticules se fait habituellement dans le cadre de la mecanique quantique. Nous allons 
travailler avec des particules scalaires instables. Le spin des neutrinos est neglige 



dans ce traitement. On verra a la section [5J| que le spin peut etre neglige si les 
masses des particules melangees sont petites par rapport a leurs impulsions. Avant 
de definir la notion d'oscillation, il faut discuter du melange de particules, qui est 
une des conditions sine qua non d'apparition des oscillations. 



1.1.1 Le melange de particules 

Le concept de melange d'etats releve des principes fondamentaux de la theorie 
quantique. Le principe de la mesure implique que la mesure d'une observable pro- 
jette l'etat quantique initial sur un vecteur propre de l'operateur associe a cet ob- 
servable. II est done necessaire de choisir une base d'etats physiques constitute de 
vecteurs propres de l'operateur en question. Par consequent, le choix de mesurer 
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une caracteristique particuliere d'un etat quantique determine automatiquement le 
choix de la base des etats physiques. Les differentes bases associees aux differents 
operateurs sont reliees par des transformations lineaires ; un etat developpe dans 
une base donnee peut etre vu comme un melange des etats propres de la base en 
question. Un exemple bien connu de melange est donne par l'experience de Stern- 
Gerlach : si l'etat initial a son spin up selon l'axe x, la mesure de son spin selon 
l'axe y nous oblige a considerer cet etat initial dans la base de l'operateur mesurant 
le spin selon l'axe y, c'est-a-dire que l'etat initial est decompose en un melange d'un 
etat ayant son spin up selon y et d'un etat ayant son spin down selon y. Un autre 
exemple est celui du systeme K°K° ou l'on neglige la violation CP : un K° initial 
se desintegrant en etats propres sous CP (deux ou trois pions), l'identification des 
produits de disintegration equivaut a une mesure de la CP-parite. Le K° initial 
doit done etre decompose en un melange d'etats de CP-parite positive et d'etats de 
CP-parite negative. 

Le concept d' oscillation presuppose, d'une part, un melange de particules de 
masses differentes : il doit etre en principe possible de mesurer les caracteristiques 
des particules dans deux bases differentes, dont une est constitute d'etats propres 
de masse. D'autre part, la presence d'une oscillation dans la probabilite de detection 
interdit qu'une mesure soit effectuee dans la base d'etats propres de masse. Pourquoi 
cette base intervient-elle ? La raison en est qu'elle se distingue de toutes les autres 
bases par le fait qu'elle est la base d'etats propres de l'operateur devolution dans 
le temps H. Elle intervient done comme base intermediaire dans tout processus 
ou une evolution temporelle est presente, sans qu'il soit necessaire d'effectuer une 
mesure de la masse. Les differents etats propres de masse evoluent differemment, 
mais tant qu'une mesure ne les distingue pas, le principe quantique de superposition 
lineaire dicte que l'amplitude totale du processus est la somme des amplitudes des 
processus indistinguables. Dans ces circonstances, une oscillation dans la probabilite 
de detection de la particule apparaitra en raison de l'interference entre les amplitudes 
partielles correspondant aux differents etats propres de masse. En conclusion, la 
probabilite de detection d'une particule se propageant presentera des oscillations 
spatio-temporelles si l'etat initial de la particule et son etat final ne sont pas des 
etats propres de masse et s'il est impossible de determiner quel etat propre de masse 
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a contribue comme etat intermediate (ce qui implique que les differences de masse 
doivent etre petites). 

Precisons les deux bases pertinentes au phenomene d'oscillation. L'hamiltonien 
doit pouvoir etre subdivise en un hamiltonien de propagation H propag decrivant la 
propagation libre des particules et un hamiltonien d'interaction Hi nt decrivant les 
interactions produisant ces particules et eventuellement leur detection. Ces deux 
parties de l'hamiltonien pourront etre distinguees s'il existe un nombre quantique, 
appele saveur, preserve par H int mais viole par H propag . La base appropriee a la 
production de particules est constitute de vecteurs propres simultanes de H int et de 
l'operateur associe au nombre quantique de saveur. Elle est appelee base de saveur. 
La base appropriee a la description de la propagation des particules est celle dans 
laquelle H propag est diagonal. Elle est appelee base de propagation ou base de masse. 
Elle ne coincide pas avec la base de saveur puisque l'operateur de saveur ne commute 

paS aveC H propag . 

Le cas le plus simple est celui des neutrinos pris comme particules stables, ou 
l'hamiltonien H propag contient la matrice de masse tandis que l'hamiltonien Hi nt 
contient les interactions faibles des neutrinos avec les bosons W et Z. La saveur 
est dans ce cas le nombre leptonique qui est electronique, muonique ou tauique et 
eventuellement sterile. II est conserve par l'interaction faible, mais est viole par la 
matrice de masse qui melange les neutrinos de saveur definie. 

Un cas plus complique est celui des particules instables, ou une partie des in- 
teractions est incorporee dans H propag pour decrire la disintegration de la particule. 
Considerons par exemple les kaons neutres, qui peuvent etre produits par les inter- 
actions fortes mais ne peuvent se desintegrer que par interaction faible. Dans une 
premiere etape, l'hamiltonien totalQ contenant toutes les interactions des kaons est 
exprime comme la somme d'un hamiltonien Hp contenant les interactions fortes 
et d'un hamiltonien Hf contenant les interactions faibles. L'hamiltonien Hf peut 
etre vu comme une perturbation de l'hamiltonien Hp. La base de saveur est definie 
comme l'ensemble des vecteurs propres de Hp. La saveur est ici le nombre quan- 



1 II s'agit en fait d'un hamiltonien effectif qui peut etre construit dans le cadre de la theorie 
chirale perturbative. On definira plus precisement ce que Ton entend par interactions fortes et 
faibles dans ce contexte lorsque Ton abordera les melanges en theorie des champs au chapitre [| 



12 



CHAPITRE 1 . OSCILLATIONS EN MECANIQ UE Q UANTIQ UE 



tique d'etrangete, preserve par les interactions fortes mais viole par les interactions 
faibles. II vaut S = +1 pour le K° et S = — 1 pour le K°. A l'ordre zero de la 
perturbation de Hp par Hf, les kaons neutres ont la meme masse, sont stables et 
n'interagissent pas entre eux. Dans une seconde etape, avec le but de decrire la prop- 



agation et la disintegration de ces particules, la methode de Wigner-Weisskopf 
est utilisee pour construire un hamiltonien effectif non hermitien H propag restreint 
aux interactions entre les kaons neutres. L'hamiltonien H propag contient d'une part 
la matrice de masse de ces particules dans la limite ou Hf est neglige par rapport 
a Hp et d'autre part les interactions faibles de Hf qui permettent les transitions 
entre kaons neutres dans la base de saveur. La non hermiticite de H propag exprime le 
fait que le systeme des kaons n'evolue pas en isolation en raison des disintegrations 
possibles des kaons vers les autres etats propres de Hp. L'hamiltonien H int contient 
les interactions fortes des kaons avec les autres etats propres de Hp. La construction 
de Hp r0 p ag et de H int peut se decrire schematiquement comme 

Htotal = Hp + Hf 

Hint ~\~ H masse -\- Hf > Ht t a i^ effectif Hj, n t -\- H propa g . 

La relation entre la base de propagation et la base de saveur peut s'ecrire 

3 

ou V est la matrice diagonalisant H propag si cet hamiltonien a ete ecrit au depart 
dans la base de saveur. Les |^(0) > sont les etats propres de propagation au temps 
t = 0, de masse et largeur defmies, appartenant a la base de propagation tandis que 
les |z/ Q (0) > sont les etats propres appartenant a la base de saveur. 

La forme de la matrice V depend du type d'hamiltonien H propag decrivant la 
propagation de la particule. Dans le premier cas, ou les particules oscillantes sont 
stables et leur lagrangien est connu, la matrice V provient de la diagonalisation de 
la matrice de masse du lagrangien et est unitaire (une discussion plus elaboree de 
cette diagonalisation pour les fermions figure a la section |5.6|) . Dans le deuxieme 
cas, ou les particules oscillantes sont instables, nous avons vu ci-dessus que l'hamil- 
tonien H propag est necessairement non hermitien pour engendrer une evolution non 
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unitaire, c'est-a-dire pour permettre la disintegration des particules. La matrice V 
diagonalisant cet hamiltonien n'est en general pas unitaire, bien qu'elle puisse l'etre 
si certaines symetries imposent sufhsamment de contraintes sur 1'hamiltonien, de 
sorte que la matrice correspondante soit normale (c'est-a-dire qu'elle commute avec 
son hermit ienne conjuguee). 

II existe encore une autre origine possible de la matrice V, qui consiste a tenir 
compte dans l'liamiltonien des interactions de la particule avec le milieu dans lequel 
elle se propage (effet MSW P7fl). On obtient dans ce cas un hamiltonien effectif 
satisfaisant a une equation de type Schrodinger. Nous ne discuterons pas de ce 
mecanisme qui est tout a fait different, puisque la particule oscillante ne se propage 
pas librement. 

Pour calculer une amplitude, nous devons pouvoir calculer le produit scalaire de 
deux vecteurs de base. Qu'en est-il des relations d'orthogonalite dans ces deux bases ? 
Dans la base de saveur, les relations d'orthogonalite suivantes sont satisfaites : 

<^(0)|^(0)>=5^, (1.2) 

puisque H int est hermitien et qu'il preserve la saveur. Dans la base de propagation, 
H P ropag n'est hermitien que si les particules sont stables. Le produit scalaire entre les 
elements de cette base doit done etre defini en general a partir de la base de saveur : 

<u k (0)\u J (0)>=J2Vr 1 V; k 1 ^. (1.3) 

7 

Ces produits scalaires deviennent orthogonaux si la matrice V est unitaire, ce qui 
est le cas pour des particules stables. Comme on l'a dit plus haut, la matrice V peut 
aussi etre unitaire pour des particules instables si une symetrie impose suffisamment 
de contraintes sur H propag . Dans le systeme K°K°, la symetrie CP sufht pour que 
V soit unitaire. Comme la violation de cette symetrie est de l'ordre du millieme, le 
membre de droite de l'equation ( |1.3|) est aussi de l'ordre du millieme pour les kaons 
neutres. 

Notons cependant que les relations d'orthogonalite dans la base de saveur ne 
sont valables que si tous les etats propres de masse sont energetiquement permis 



28 1 . Par exemple, supposons qu'il existe quatre neutrinos stables tels que rrii = 
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pour i = 1, 2, 3 et m 4 > 1 GeV. Pour des neutrinos d'energie inferieure a 1 GeV, on 
aura 

Ce probleme est une premiere indication des dimcultes liees a la definition de la base 
de saveur. Par exemple, est-il possible d'associer une particule observable a chaque 
etat de saveur ? Remettons a plus tard le cas des particules instables : elles ne corres- 
pondent en toute rigueur a aucun etat de l'espace de Hilbert des etats physiques 
mais apparaissent en theorie des champs comme des poles de la matrice S. Par 
contre, a chaque particule stable est associee en mecanique quantique relativiste une 
representation du groupe de Poincare caracterisee par une masse et un spin definis. 
Dans ce cadre, les etats de saveur ne correspondent pas a des particules, a moins 
d'etre degeneres. Si Ton examine le probleme en theorie des champs, on remarque 
que les operateurs de creation et d'annihilation d'etats propres de masse sont bien 
definis. Par contre, le changement de base ( |1.1| ) applique a ces operateurs ne fournit 
pas des operateurs de creation et d'annihilation d'etats de saveur satisfaisant aux 
relations de commutation canoniques [j30| . On pourrait se demander quelle realite 
attribuer aux etats de saveur apparaissant dans l'equation ( |1 . 1|) . II se fait que dans 
un processus ou des oscillations sont presentes dans la propagation, les particules 
oscillantes sont uniquement identifiees a la production et a la detection dans une 
superposition d'etats propres de masse formellement equivalente k un etat de saveur. 
Dans ce sens-la, les etats de saveur sont observes dans ce type d'experience tandis 
que les etats propres de masse ne le sont pas. 



1.1.2 Oscillations de la probability de detection 

Partant d'un etat propre d'interaction de saveur a, l'amplitude de detecter au 
temps t un etat de saveur (3 est donnee selon les regies de la mecanique quantique 
par 

A(a^P,t) = <vp(0)\exp(-iH propag t)\v a (0)> 

= < M Q ) I E Vaj exp (-i tyt - Tjt/2) | ^(0) > 
j 
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= < ^(0) I E V*d ex P H E jt ~ T J t / 2 ) Vr 1 1 z/ 7 (0) > . 

Les symboles Ej, rrij et designent les energie, masse et largeur de l'etat propre 
de masse | Vj(t)>. L'amplitude peut etre reecrite comme 

A{a -/3,*) = E V aj exp {-i Ejt - r,-t/2) V~ p l . (1.4) 
j 

On a utilise la relation d'orthogonalite ( |1.2|) dans la base de saveur, puisque la base de 
propagation ([L3|) n'est en general pas orthogonale. Cette non orthogonalite implique 
que les etats 1 1^(0) > peuvent etre normalises de differentes fagons. Cela n'a pas d'im- 
portance tant qu'ils n'apparaissent pas dans le resultat final du calcul de l'amplitude, 
comme on le voit clairement ci-dessus ou la formule d'oscillation est calculee sans 
meme poser la question de cette normalisation. Cependant, cette question reapparait 
lorsque Ton tente d'ecrire des amplitudes impliquant un etat propre de masse comme 
etat initial ou final. Par exemple, l'amplitude M.(Kl-^-kti) =< nn \H tota i\ Kl > 
depend de la normalisation choisie pour l'etat | Kl >■ On verra au chapitre |] la 
maniere correcte de definir de telles amplitudes en theorie des champs. 

La formule ( |1.4j ) de l'amplitude d'oscillation est souvent ecrite |TL| en termes 
des <u k (0) | z^(0)> : 

A(a /3, t) = E V aj exp {-i Ejt - Tjt/2) V k \ < u k {0) | ^(0) > , (1.5) 

mais cette expression est equivalente a la formule ( |1.4|) . En effet, le remplacement 
dans l'equation ci-dessus des <z/ fc (0) | z/,(0)> par leurs expressions ( |1.3| ) donne 

A(a - ft f ) = E Vaj exp (-i Ejt - Tjt/2) V h \ Vff V~? f 

j,k,l 

= E V *i ex P H Ejt - Tjt/2) V- p l , 
j 

et l'on retrouve la formule ( |1.4| ). Ce n'est pas surprenant, puisqu'a travers l'equation 
(|1.3|), ce sont les relations ( |1.2| ) d'orthogonalite dans la base de saveur qui sont 
utilisees. 

Le resultat est egalement identique a celui obtenu par l'utilisation de deux 
bases de saveur |z/ Q ,(0)>j n et ou t<v a (0) |, reliees par l'operateur de renversement 
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dans le temps T |32|, [33], [J3Tj. Les ket appartiennent a la base in tandis que les 
bra appartiennent a la base out. On ne se sert pas des vecteurs hermitiens con- 
jugues. De nouvelles bases \vj(0) ^in 6t out 

<z/j(0) | d'etats propres de masse sont 
definies telles qu'elles diagonalisent rhamiltonien H propag , c'est-a-dire que les nom- 
bres ou t <Uj(0) \H propag \ Uj(0) >i n forment une matrice diagonale : 



| ^a(O) >in = V <*3 I ^i( ) >in , 

j 

out<Va(0)\ = out < Vj (0) | 



;i.e) 



mais les bases d'etats propres de masse ne sont pas reliees par l'operateur T. Les 
etats propres de masse in et out sont orthogonaux : 



out 



<fj(0) I Vi(0)> in = 5, 



13 ■ 



La formule d'oscillation obtenue par cette methode est identique a l'equation ( |1.4j ). 
De nouveau, ce n'est pas etonnant puisque les nouvelles bases in et out sont definies 
a partir de la base de saveur. Ce procede peut sembler artificiel : quelle est la 
signification de cette nouvelle base ? Peut-on prendre le produit scalaire entre ces 
deux bases qui ne sont meme pas duales sous l'operateur T ? Pourquoi introduire une 
nouvelle base alors que les etats propres de propagation ne sont jamais directement 
observes ? 

Notons comme cas particulier que si H pTopa g est hermitien, alors V est unitaire 
et Fj = 0. L'amplitude s'ecrit dans ce cas comme 

A(a -+P,t)=Y, V *i ex P H Ejt) . 



Cette expression est couramment utilisee pour les melanges de neutrinos stables. 

Dans le cas general d'un hamiltonien H propag non hermitien, la probabilite de 
detection est donnee par la norme au carre de l'amplitude d'oscillation (|1.4j ) et s'ecrit 

2 



3 

j,k 



;i-7) 
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oil Tj k = (Tj + T k )/2 et AEj k = Ej — E k . On peut aussi la reecrire comme 
V(a^P,t)=Y,e- r i t \V aj \ 2 \Vri\ 2 



+ 2 £ e- T ^Ve (v a3 V^v: k (V- l )ip) cos(AE jk t) 

j<k 

+ 2 £ e~ r ^lm (V aj V^V: k {V- x T kp ) sin(AE jk t) , 



Cette probabilite de detection est une somme de termes qui decroissent tous comme 
des exponentielles, soit en e~ Tjt , soit en e~ Tjkt . Certains de ces termes oscillent en 
cos(A£ j k t) et d'autres en sm(AEj k t). C'est a ces comportements oscillatoires que 
Ton pense lorsque Ton parle d'oscillation de la probabilite de detection. 

Quant a l'oscillation des antiparticules, le theoreme CPT applique a l'amplitude 
(|1.4j) donne la relation 

A(a -> P,t) = A{(3 -+ct,t), 

avec la consequence que A(a = A(a — > a,t). La violation CP dans les os- 

cillations ne peut done etre mesuree que s'il y a changement de saveur (a ^ (3), 
puisqu'elle se manifeste par une difference entre \ A(at — > (3, t)\ 2 et \A(a — > j3, t)\ 2 . Si 
la matrice de diagonalisation V est unitaire, seul le terme Tm(V Q j V* k (V" -1 )^) 



(invariant de Jarlskog P5| ) viole la symetrie CP dans l'expression de la probabilite 
(|1.8|) . En effet, le passage de V{a — > (3, t) a V(a P,t), qui revient a echanger a et 
P, est equivalent a remplacer V par V* si V est unitaire . Cette substitution change 
le signe de 2m{V a j V* k (V^ 1 )^) mais laisse les autres termes de la probabilite 
invariants. 

Les deux cas dans lesquels la formule ( |1.7|) est utilisee, - mesons K et B, neutri- 
nos -, necessitent un traitement relativiste, que n'offre pas cette formule, puisque le 
facteur AEj k t n'est pas un invariant sous les transformations de Lorentz. Cette non 
invariance resulte du fait que l'on a neglige la dependance spatiale des etats dans 
la derivation de la probabilite de detection. Pourtant, la difference entre le lieu de 
production et de detection des particules introduit une phase supplementaire dont 
il faut en principe tenir compte lorsque Ton etudie des phenomenes d'interference. 
Void les solutions couramment apportees a ce probleme. Le calcul de la probabilite 
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( |1.7j ) est en tout cas valable si la particule oscillante est considered dans son repere 
au repos. On remplace alors la phase — i Ejt par la phase correspondante dans le 
repere au repos, c'est-a-dire —i rrijT, oil r est le temps propre de la particule. Re- 
marquons que si la particule a une tri-impulsion p^- bien definie, cette phase peut 
aussi s'ecrire 

— i rrijT = —i Ejt + i pj • x . (1.9) 

Le membre de droite peut aussi etre calcule en resolvant l'equation de Schrodinger 
pour une onde plane, c'est-a-dire en tenant compte de la dependance spatiale des le 
depart dans le calcul de l'amplitude ( |1.4| ). Traditionnellement, on utilise le membre 
de gauche de l'equation ( |1.9|) pour les kaons et le membre de droite pour les neutrinos. 

Dans le cas des mesons K et B, on passe dans le repere du laboratoire par la 
relation t = / -/t = Er/m, oil m est choisi arbitrairement comme m = (rrij + m,fc)/2 
et E est l'energie de la particule. Si Ton note la phase du terme d'interference par 
$, la probability de detection ( |L7p contient des termes oscillant en 



exp (— «$) = exp (—i(irij — m^) r 

? 

E 



m 



exp ( —i(mj — rrik) — t 



( m 2 - — ml N 

exp — i — t 

y \ 2E 



et ces termes decroissent en 



exp (-T jk r) = exp y-T jk ^ t 

Comme les experiences ne mesurent que la dependance spatiale de la probability 
de detection, il faut passer a une expression ne dependant que de la distance. On 
substitue done dans ces expressions la formule de la trajectoire classique, t = EL/ P, 
oil P est l'impulsion de la particule. Definissant la longueur d 'oscillation pour le 
melange jk par 

47T P 

I™ EE -= . , (1.10) 

3h m] - ml 

la probability de detection ( |1.7| ) contient des termes oscillant en 

exp (-z$) = exp (^-2iri-^J , 
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et ces termes decroissent en 

m 



exp ( -T jk — L 



Les oscillations sont observables si L/l° s k c « 0(1) et si 0(l°f. c ) < P/(mTj k ), ou en- 
core si 0(rrij — rrik) > 0(Tjf.)- Cette derniere condition signifie banalement qu'on 
ne verra pas d'oscillations si la particule se desintegre avant d'osciller notab lenient. 

On peut deja se rendre compte des problemes lies a ce type d'extension rela- 
tiviste. Si Ton avait defini un temps propre pour chaque etat propre de masse, la 
phase $ serait devenue 

$ = m . T . _ mkU = ( m | _ m J) t/E j 

ou Ton a pris la meme energie pour les etats j et k, faute de mieux. La longueur 
d'oscillation correspondante est alors deux fois plus courte |57| ! Le probleme reside 



dans la question de l'existence d'un repere au repos d'une particule qui n'a pas 
de masse bien definie. II existe d'autres possibilites de se tromper en prenant des 
energies differentes, ou des temps differents, ou des vitesses differentes, ou des im- 
pulsions differentes pour les etats j et k. Les arguments utilises pour justifier le 
bon choix EM HD] ne sont pas vraiment convaincants. Pour des processus du 



type 0(1020) -^K°K°, ou deux particules oscillent simultanement, les ambigui'tes se 
multiplient ! 

La situation ne s'ameliore pas si Ton travaille avec le membre de droite de 
l'equation ( |1.9|) , ce qui est la procedure couramment suivie dans le cas des neutri- 
nos. En effet, il faut se donner des prescriptions pour manipuler l'energie-impulsion 
(Ej,pj), le temps t et la distance x. La solution couramment adoptee fHf , tout a 



fait ad hoc, est de supposer l'egalite des tri-impulsions, pj = p k = p, de sorte que 



2 2 

m A — m 



E 3 = ^p* + m] = E + ^ E 



ou m = (rrij + m&)/2, E = ^p 2 + m' 2 et Ton a suppose que \rrij — m k \ <^ E. La 
probabilite ( |1.7| ) oscille alors en 



exp (-iAEj k t) = exp f -% 3 — h t 
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Le raisonnement se poursuit alors comme pour la particule instable, avec la meme 
definition de la longueur d'oscillation. line deuxieme ecole f|2|, |43], £|4| utilise des 
energies et tri-impulsions differentes pour les differents etats propres de masse, et 
les calcule par conservation de l'energie-impulsion a la production. Le calcul de la 
phase $ du terme d'interference donne la meme longueur d'oscillation que l'equation 
fll.lUp si le temps t est remplace par t = EL/ P quel que soit l'etat propre de masse. 
Les ambiguites lors du passage d'une dependance temporelle a une dependance 
uniquement spatiale persistent dans ces deux traitements. On pourrait par exemple 
argumenter que les etats propres de masse ont des vitesses differentes et arrivent done 
au detecteur en des temps tj et tk legerement differents. Si Ton pose l'egalite des tri- 
impulsions, la longueur d'oscillation est alors deux fois plus courte que l'expression 
(|f .1 Op . Enfin, une troisieme possibility |4l| consiste a egaler les energies et a calculer 
les tri-impulsions par la formule p| = E 2 — mh. Dans ce cas, on obtient directement 
une dependance spatiale dans le terme d'oscillation, mais la prescription d'egalite 
des energies parait aussi ad hoc que la prescription d'egalite des tri-impulsions. 

II est pratique d'ecrire la longueur d'oscillation pour les neutrinos sous la forme 

!™~2 1?. m P(MeV) 

qui permet d'evaluer facilement la difference de masse necessaire pour une longueur 
d'oscillation donnee. 



1.2 Critique de la derivation traditionnelle de la 
formule d'oscillation 

La derivation expliquee a la section precedente de la probability de detection 
d'une particule en melange est incorrecte pour toute une serie de raisons. On peut 
emettre les objections suivantes : 

f. L 'interference entre etats de masses differentes est interdite en mecanique 
quantique non relativiste. En effet, l'invariance de l'equation de Schrodinger 
sous le groupe de Galilee fixe la loi de transformation d'un etat : il est multiplie 
par un facteur de phase qui depend de la masse de l'etat, de la position et du 
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temps. Des etats de masses differentes se transforment done differemment et 
la phase relative d'une superposition d'etats n'est pas conservee sous le groupe 
de Galilee. Une superposition coherente de ces etats est done exclue (regie de 
superposition de Bargmann fl6|). Pour la meme raison, les particules instables 
ne peuvent etre decrites en mecanique quantique non relativiste car il n'existe 
pas de transition entre etats de masses differentes. 

2. II faut de toute fagon tenir compte des corrections relativistes pour decrire 
les experiences. Habituellement, on part de la formule non relativiste dans 
le repere au repos de la particule et on effectue un boost pour passer dans 



le repere du laboratoire comme explique ci-dessus |3q| . Cette procedure est 
ambigue lorsque Ton examine l'interference entre etats de masses et/ou de 
largeurs differentes : existe-t-il un repere au repos de l'etat oscillant, ou en 
faut-il deux, un pour chaque etat propre de masse ? Cette ambiguite a suscite 



des controverses sur la question des oscillations [[37], |39[. II vaut mieux utiliser 
un formalisme relativiste des le depart, e'est-a-dire recourir a la theorie des 
champs, puis effectuer des approximations non relativistes si necessaire. 

3. Les etats propres d'interaction (ou de saveur) sont mal dermis. II n'existe pas 
d'espace de Fock pour des etats qui n'ont pas de masse bien definie, e'est-a-dire 
qu'il n'existe pas d'operateurs de creation et d' annihilation d'etats de saveur 
satisfaisant aux relations de commutation canoniques |3t| . Une tentative a ete 
faite en theorie des champs non perturbative, a l'aide de transformations de 
Bogoliubov mais elle mene a de serieux problemes d'interpretation et brise la 



covariance sous les transformations de Lorentz [f47f| . 

Pour les particules instables, un probleme supplementaire apparait puisque 
l'hamiltonien effectif n'est pas hermitien. Comme on l'a vu ci-dessus, les etats 
propres de masse ne sont plus orthogonaux avec leurs hermitiens conjugues et 
des problemes de normalisation se posent. Une solution possible (voir equation 
(|l.(j|)) consiste a utiliser deux bases in et out a la place d'une seule base (et son 



hermitienne conjuguee) |32], [33], |34[. Neanmoins, cette methode peut paraitre 
artificielle. La question de la normalisation a de l'importance car dans le 
systeme K°K°, l'incertitude sur la normalisation des etats est du meme ordre 
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que la violation CP directe H8I . 

5. Les particules instables n'apparaissent pas comme etats asymptotiques : ceux- 
ci ne sont done pas contenus dans l'espace de Hilbert des etats physiques. 
En fait, elles ne peuvent etre traitees de fagon coherente que comme etats 
intermediaires dans le cadre de la theorie quantique des champs fl49"fl . 

6. Les etats oscillants etudies jusqu'a present (mesons K et B, neutrinos) ne 
sont pas observes directement car ils sont electriquement neutres. Seuls les 
particules avec lesquelles ils interagissent (ou dans lesquelles ils se desintegrent) 
par interaction faible sont observes. Dans une approche realiste, on ne devrait 
calculer que des amplitudes de transition entre etats observables. 

7. Les experiences etudient la propagation spatiale, et non temporelle. L'utilisa- 
tion de la formule classique L = vt dans un calcul eminemment quantique pour 
remplacer le temps par la distance laisse songeur. D'ailleurs, ne faudrait-il pas 
utiliser des vitesses v differentes pour les differents etats propres de masse ? Ce 
probleme est lie a la question mentionnee ci-dessus du passage a un formalisme 
relativiste. Cette ambiguite mene a des differences du simple au double dans 



le calcul de la frequence d'oscillation |37|, 

8. La derivation traditionnelle implique une energie-impulsion parfaitement de- 
finie, e'est-a-dire que les etats sont des ondes planes. Par le principe d'incer- 
titude, leur localisation spatio-temporelle est alors inconnue. Comment parler 
dans ce cas de propagation sur une distance definie ? 

9. II n'y a pas de raison de penser que les energies ou les impulsions des etats 
oscillants soient egales. Cette question est d'ailleurs ambigue : peut-on parler 
de l'energie ou de l'impulsion d'un etat propre de masse lorsqu'il est superpose 
lineairement a un autre etat propre de masse ? L'hypothese qui consiste a egaler 
les tri-impulsions des etats propres de masse des neutrinos est done injustifiee, 
comme d'ailleurs aussi l'hypothese opposee qui consiste a egaler leurs energies. 
D'ailleurs, la conservation de l'energie- impulsion a la production n'impose-t- 
elle pas l'egalite de l'energie-impulsion des differents etats propres de masse ? 
Ou alors, faut-il considerer que l'energie-impulsion des autres particules emises 
a la production de la particule oscillante varie selon l'etat propre de masse 
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emis ? Cette question n'est pas purement academique. Tout d'abord elle remet 
en question la formule d'oscillation. Ensuite, l'utilisation d'une conservation 
stricte de l'energie-impulsion a la production mene a considerer les particules 
associees a la particule oscillante (par exemple le muon dans Tc—^fiu) comme 
une superposition d'etats d'energies-impulsions differentes, avec le resultat 
qu'ils oscillent eux aussi [50[ |5T||. Nous verrons au chapitre |6| que ce n'est 
pas le cas. 

10. Les criteres d'observabilite des oscillations (a part L « Z°| c ) n'apparaissent pas 
explicitement dans la formule d'oscillation ( |1 . 7| ) . Par exemple, on s'attend a 
ce que les oscillations soient inobservables 

- si la longueur d'oscillation est de l'ordre de grandeur ou plus petite que l'in- 
certitude sur la position de la source et/ou du detecteur. 

- si l'energie et l'impulsion des etats oscillants peuvent etre mesurees avec une 
precision suffisante pour determiner quel etat propre de masse se propage. 
Pour que la superposition quantique soit possible, il faut done simultanement 
une localisation sumsamment precise de la source et du detecteur et une incer- 
titude suffisante sur l'energie et/ou l'impulsion a remission et a la detection 

521. 



Toutes ces objections menent a abandonner un traitement direct des etats oscil- 
lants. Le processus d'oscillation va etre considere dans son ensemble : les etats oscil- 
lants deviennent des etats intermediaries, non observes directement, se propageant 
entre une source et un detecteur. lis sont represented en theorie des champs par 
leur propagateur. La forme de celui-ci, consistante avec la relativite restreinte et la 
causalite, determine revolution spatio-temporelle ainsi que la probabilite de disin- 
tegration si la particule est instable. Etant donne que la covariance sous les trans- 
formations de Lorentz est implicite dans cette approche, le recours a des reperes 
particuliers et des boosts est rendu inutile. Cette methode a ete proposee pour la 
premiere fois par Jacob et Sachs || pour les particules scalaires instables et Giunti, 



Kim et Lee |53[ pour les neutrinos. 

Dans un premier temps, nous allons presenter un modele simplifie suffisant a 
decrire la propagation d'une seule particule stable. Dans cet exemple, l'etude de la 
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dependance spatiale de l'amplitude globale de propagation se reduit a l'analyse de 
la transformed de Fourier du propagateur par rapport a sa tri-impulsion. L'esprit de 
cette methode est proche de Schwinger [Q qui etudie la dependance temporelle de 
l'amplitude. On appliquera aussi ce modele a la description de particules instables, 
bien qu'en toute rigueur il ne convienne pas a cette situation. On estimera les con- 
tributions non exponentielles a l'amplitude de propagation, et l'on calculera aussi 
une probability d'oscillation dependant uniquement de la distance. Ses insuffisances 
vis-a-vis des particules instables mises a part, le modele simplifie ne permet pas non 
plus de repondre a toutes les questions posees ci-dessus a propos de la derivation tra- 
ditionnelle des oscillations en mecanique quantique. Je pense d'une part au probleme 
de l'egalite ou non des energies-impulsions des etats oscillants et d'autre part aux 
conditions d'observabilite de ce phenomene. 

Dans un second temps, nous proposerons un modele plus sophistique. Pour 
tenir compte des restrictions experiment ales sur les oscillations, les etats entrants et 
sortants a la source et au detecteur seront representes par des paquets d'ondes. On 
pourra ainsi deriver les conditions d'observabilite des oscillations. La propagation 
de particules instables pourra aussi etre decrite de fagon coherente. Enfin, toutes 
les questions soulevees ci-dessus concernant les oscillations trouveront une reponse. 
Les formules theoriques seront appliquees dans une serie de cas concrets concernant 
des experiences actuelles. Une parametrisation coherente de la violation CP sera 
aussi etablie. Tout au long de ce travail, la symetrie CPT sera supposee exacte. 
Parametriser une violation eventuelle de cette symetrie dans un modele comme le 
notre base sur la theorie des champs est incoherent puisque personne n'a reussi 
jusqu'a present a construire une theorie des champs relativiste violant CPT. Autant 
parametriser cette violation dans le cadre de la mecanique quantique, comme cela 
a deja ete fait de nombreuses fois (toutefois, la premiere parametrisation, a ma 
connaissance, de la violation de CPT dans le systeme K°K° a ete effectuee en 
theorie des champs ! ) . 



Chapitre 2 

Le propagateur complet 



II apparait dans ce chapitre que le propagateur complet constitue l'element cle 
de toute description d'un phenomene de propagation macroscopique en theorie des 
champs. Remarquons tout d'abord que l'amplitude de propagation d'une particule 
depend en general des processus de production et de detection. Sous certaines con- 
ditions analysees dans le chapitre [|, il est possible de factoriser l'amplitude globale 
en trois parties : production, propagation et detection. La propagation d'une par- 
ticule peut alors etre etudiee a partir du seul propagateur, sans prendre en compte 
les etats entrants et sortants. Dans ce cas, l'utilisation de la theorie quantique des 
champs n'est pas plus compliquee que celle de la mecanique quantique et fournit 
directement une formule d'oscillation spatiale, ce qui est deja un plus par rapport a 
la formule ( |1.7| ) derivee en mecanique quantique. Notons cependant que nos calculs 
ne s'appliqueront pas a la propagation de bosons de jauge charges : le propagateur 
n'etant pas invariant de jauge, il ne peut etre considere independamment. Un tel 
processus ne peut etre factorise qu'en sous-processus invariants de jauge |5|J . 

Apres le passage en revue des hypotheses du modele simplifie, nous etudierons 
la representation spectrale du propagateur complet et les consequences de celle-ci 
sur l'amplitude de propagation spatiale, tant pour une particule stable que pour 
une particule instable. En particulier, les corrections non exponentielles dues aux 
seuils de production de plusieurs particules reelles seront examinees. Le chapitre se 
termine par un exemple de calcul de propagateur complet. 



25 



26 



CHAPITRE 2. LE PROPAGATEUR COMPLET 



2.1 Propagation : le modele simplifie 

Dans cette section, nous presentons un modele simplifie tel que l'amplitude de 
propagation se factorise et l'etude de la propagation se resume a celle du propaga- 
teur. Le processus de propagation d'une particule scalaire est symbolise par 

Pj{q) ^ P F (k) + i/(p) 

\ 

u(p) + D I (q') {tD -^ ] Dp(k') 

Pi represente l'ensemble des particules arrivant dans la region de production 
centree autour du point (t P , x P ), d'energie-impulsion totale q tandis que Pp represente 
l'ensemble des particules issues de la region de production et d'energie-impulsion 
totale k, a l'exception de la particule intermediaire v dont on etudie la propaga- 
tionQ. Dj, Dp et (t D ,x D ) ont des definitions similaires mais concernent le point de 
detection. Le point d'interaction a la production est note par x et le point d'inter- 
action a la detection est note x'. On pose que l'etat intermediaire se propageant de 
x en x' a les nombres quantiques d'une particule, pas d'une antiparticule. Toutes 
les particules du processus sont supposees scalaires. Les conditions experiment ales 
(energies-impulsions de Pi,f, D if ) sont choisies telles que le processus ci-dessus, 
ou une particule quasi-reelle v se propage macroscopiquement, est selectionne parmi 
tous les processus possibles. En d'autres termes, les energies-impulsions des etats ini- 
tiaux et finaux sont telles que l'element de matrice S correspondant a ce processus 
est evalue pres de la singular ite du propagateur de la particule intermediaire. 

Ce processus pourrait aussi etre symbolise par le diagramme suivant, sur lequel 
les regions de production et de detection sont indiquees par des cercles en pointille 
et les fleches indiquent la direction des impulsions : 



lr rout au long de ce travail, les vecteurs en 3 dimensions seront notes en gras tandis que les 
vecteurs en 4 dimensions seront notes normalement. Les produits scalaires s'ecrivent respectivement 
p • q et p • q. La metrique utilisee est g^ v = diag{\, —1, —1, —1). 
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27 



/ \ 




Les hypotheses simplificatrices de ce modele sont les suivantes : 

- tous les etats entrants et sortants sont stationnaires, c'est-a-dire que leurs 
energies ont des valeurs determinees. 

- un etat au point de production et un etat au point de detection sont localises 
spatialement a une precision infinie. 

- les autres etats ont des tri-impulsions determinees a une precision infinie, c'est- 
a-dire qu'ils sont des ondes planes dans l'espace de configuration. 

Pour fixer les idees, supposons que Pi et Di soient chacun constitues d'un seul 
etat stationnaire designe respect ivement dans l'espace de configuration par 

^ (x-x P ) e'^ et fo, (x'-x D ) e~ iE ^ . 

On ne pourra parler de propagation macroscopique que si les interactions a la pro- 
duction et a la detection sont suffisamment localisees par rapport a la distance de 
propagation. C'est pourquoi on a pose l'hypothese extreme que les fonctions d'onde 
des etats Pi et Di sont des distributions delta de sorte que les interactions de pro- 
duction et detection soient localisees spatialement a une precision infinie : 

0^ (x-xp) ~ 5 (3) (x-x P ) et fo, (x'-x D ) ~ 5 (3) (x'-x D ) . 

Bien sur, le prix a payer est une incertitude complete sur la tri-impulsion de la 
particule intermediaire. Quant aux etats Pp et D F , on suppose qu'ils ne contiennent 
chacun qu'une particule, qui est representee par une onde plane : 

<Mx,t) ~ e -^* +ik - x et <h F (x,t') ~ e-^f-Mk'-x' _ 



28 



CHAPITRE 2. LE PROPAGATEUR COMPLET 



L'amplitude du processus s'ecrit dans l'espace de configuration comme 
A ~ M P M D J d*x J d 4 x'^ (x-Xp) e* Ep i t fo I (x'-x D ) e"^*' 



j4 

x gl^t-ik-Xgl^t'-lK'-X' I _?LP_ e -ip-(x'-x) Q( p 2^ (2.1) 



(27T) 

ou Mp, Md sont les couplages, supposes constants et k, k' sont les impulsions totales 
des etats Pp et D F . Le propagateur est symbolise par G{p 2 ) ou par G(p°,p) selon 
les besoins. Cette formule est justifiee rigoureusement au chapitre f|, dans le cadre 
du modele sophistique. L'integration sur t et i! donne 

A ~ MpMp) J d 3 x J d 3 x f J -^5(E Pl -Ep F -p ) 5(e Di -E Df + p ° 



(2ny 

x e- ik - x - ik ' x '0 P/ (x-xp) D/ (x'-x D ) e ip - (x '- x) G(p 2 ) 

d 3 p 

xe- !k ' x - ik '' x 'S(x-x P ) &T,(x'-x D ) e^^'^G^p), 



M p M d 2tc5{E Di - E Df + Epj- E Pp ) J d 3 x J d 3 x' J 



ou E = E Pl — E Pp >0. L'integration sur x et x' donne, apres insertion des fonctions 
d'onde delta 



A~M p M d 2k8{E Di -E Df +E Pi -E Pf ) J ^e* pL G(£,p 



(2.2) 



ou L = x D -x P . 

On voit que dans ce modele simplifie, l'amplitude de propagation est directe- 
ment proportionnelle a la transformee de Fourier du propagateur sur sa tri-impulsion 
et aura done une dependance spatiale en L. Notons que Ton etudie habituellement 
la transformee de Fourier du propagateur sur l'energie |4|, || ce qui donne une am- 
plitude dependant du temps. Outre le fait qu'aucun modele n'est donne pour jus- 
tifier ce choix, il est plus logique d'etudier la dependance spatiale du propagateur 
puisqu'aucune experience n'observe directement la dependance temporelle. 

Ce modele-ci presente cependant plusieurs defauts : 

- il ne s'applique pas a la disintegration des particules instables, pour lesquelles 
il ne peut etre question d'etat stationnaire lors de leur disintegration. 



2.2. DEUX FORMES DE REPRESENTATIONS SPECTRALES 



29 



l'approximation, soit d'etats stationnaires, soit d'ondes planes, est trop forte 
car elle ne permet pas d'etudier l'inffuence des conditions de production et de 



detection. De plus, comme on le verra a la section £374], elle est trop imprecise 
pour que Ton puisse repondre par la suite a la question d'egalite ou non des 
energies-impulsions des etats oscillants. 

- l'amplitude est independante de la direction de L. 

Ces defauts ne deviennent preoccupants que pour l'etude de phenomenes tels 
que les oscillations de particules, qui resultent de la compensation quasi totale de 
differentes phases. Nous y reviendrons dans le prochain chapitre et developperons 
un modele plus complet par la suite. Une autre imperfection de ce modele releve de 
la question de la possibility d'une propagation libre. Dans une chambre a bulles, la 
trajectoire d'une particule chargee est signalee par des bulles resultant des interac- 
tions electromagnetiques de la particule avec le milieu de propagation |]]. II faudrait 
en principe tenir compte de ces mesures successives dans le formalisme quantique. 
Ce probleme ne sera pas examine ici puisque l'on vise a decrire la propagation de 
particules neutres (kaons, neutrinos). La question de principe subsiste cependant, 
car ces particules interagissent encore faiblement avec le milieu de propagation. Le 
modele MSW [E7|| est un exemple ou ces interactions jouent un role determinant. 

Pour revenir au modele simplifie, il convient maintenant d'etudier la forme 
generale du propagateur et sa transformed de Fourier. 



2.2 Deux formes de representations spect rales 

Un champ relativiste decrit une excitation locale qui produit un spectre d'ener- 
gies ou de masses possedant une limite inferieure caracteristique du type de champ. 
Par exemple, la limite theorique pour n'importe quel boson ou lepton electriquement 
charge est la masse de l'electron, tandis que pour les champs neutres de ces deux 
types, le spectre commence en principe en zero. Pour les baryons, cette limite est 
la masse du proton. Ces limites sont fixees par les masses des particules stables 
restantes apres que toutes les disintegrations possibles ont eu lieu. 

La propagation de ces excitations est decrite par la fonction de Green a deux 
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points ordonnee dans le temps, qui pour un champ bosonique scalaire est definie par 

G{x-y)=<U\T{<P{x) <j>*(y))\0> . 

Elle est aussi appelee propagateur complet. En vertu des principes fondamentaux 
de la theorie relativiste des champs, cette fonction prend la forme suivante, appelee 
representation spectrale ou de Kallen-Lehmann p5[ : 



ou 



/•CO 

G(x — y) = I ds p(s) Ap(x — y, s) 
Jm 2 

d 4 p 



A F (x-y,s 



(27T) 



-ip(x-y) 



P 



s + it 



est le propagateur de Feynman du champ scalaire libre. La fonction spectrale p(s) 
est definie par 



p(s = p 2 ) = (2vr) 3 2 ^ (4) (P " Pn) I < 10(0) \n>\ 2 , 



(2.3) 



ou la sommation sur les etats \n > a n particules comprend aussi une integration 
sur leurs tri-impulsions. Cette fonction peut etre interpretee comme etant la densite 
d'etats de masse s. M est la masse invariante de l'etat a plusieurs particules le plus 
leger en interaction avec le champ. On a 

p(s) = pour s < M 2 . 

La representation spectrale du propagateur complet dans l'espace des impulsions 
s'obtient par transformee de Fourier et s'ecrit 



(2.4) 



ou l'on reconnait une sommation ponderee selon la masse sur le propagateur de la 
particule libre. 

Le propagateur G(p 2 ) peut etre vu comme la valeur limite de la fonction G(z) 
definie dans le plan complexe : 

r°° i 
G(p 2 )= lim G(z)= lim / ds p(s) . 

z^+ie z^pl+ie JM 2 Z — S 
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La fonction G(z) est analytique partout sauf sur une coupure le long de l'axe reel 
positif a partir de z = M 2 . 

Pour que les integrates defmies ci-dessus existent, il faut que la fonction p(s) 
decroisse lorsque s tend vers l'infini. Par consequent, le propagateur a un comporte- 
ment asymptotique en G(p 2 ) ~ ip~ 2 lorsque p 2 — > oo. Cette propriete caracte- 
ristique d'un champ local peut etre utilisee pour etablir une representation spec- 
trale equivalente du propagateur ||], qui est particulierement pratique pour etudier 
les particules instables et faire le lien avec l'approche perturbative des diagrammes 
de Feynman. Observons d'abord que G(z) n'a pas de zeros complexes puisque 

■ \ ■ f°° , ( x - s )p( s ) f°° , P( s ) 

G(x + iy)=i ds 4P~h +y ds 7 ^7 9 = 

Jm 2 [x — s) z + y z Jm 2 [x — s) z + y z 

implique y = 0. La fonction P(z) = iG~ 1 (z) — z est done reguliere dans le plan 
complexe, a l'exception de la coupure sur l'axe reel positif deja mentionnee. La 
fonction z~ l P{z) tend vers zero a l'infini dans le plan complexe et presente un pole 
a l'origine. Elle peut s'ecrire 

z Jm 2 z — s 

ou a(s) est une fonction dont on va determiner les proprietes. On peut encore ecrire 

iG-\z) = z -M 2 -z r ds^-. 

Jm 2 z - s 

En comparant avec la formule precedente pour G(z) (equation (|2.4|) ), on a 

iG(0) = M n 2 = I" ds ^ > . 

Jm 2 s 

De plus, on verifie que 

TieG{x + ie) = - Tie G(x - % e) = vr p(x) > , (2.5) 
TmG(x + ie) = ImG(x — ie), 



done 



G~ x (x + ie) -G _1 (x-ie) > pour x > M 2 , 
G~\x + ie) - G-\x -ie) = pour x < M 2 , 
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On verifie aussi que 



G (x + ie) — G (x — ie) = 2tt x a(x) , 

done a(x) > avec a(x) = pour x < M 2 . Ces deux conditions, Mq > et 
a(x) > , garantissent que G(z) n'a pas de pole pour z reel negatif. En conclusion, 
la deuxieme representation spectrale du propagateur s'ecrit 



(2.6) 




ou 



U(z) =z I"* da 

Jm 2 z - s 



Remarquons que 



(2.7) 



{2.8) 



XmH{p 2 ) = lim Xmli{z) = — np 2 u(p 2 ) < . 

Cette forme du propagateur se retrouve aussi a partir des diagrammes de Feyn- 
man. Soit la fonction —iU(p 2 ), appelee energie propre, designant la somme des 
diagrammes irreductibles a une particule (1PI), 



-iU(p 2 ) = 




+ 




+ 




+ 



-1PI - 



Le propagateur complet est la somme de toutes les insertions possibles de l'energie 
propre dans le propagateur : 

G(p 2 )= + — 




1PI 



1PI 



et peut s'exprimer comme une serie geometrique en — % H(p 2 ) : 
G{p2) = W^M 2 + p^M 2 H 11 ^) W^Ml 



+ 



p 2 - M 2 - U(p 2 ) ' 



+ 



(2.9) 



2.3. PARTICULE STABLE : REPRESENTATION SPECTRALE 
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On retrouve la deuxieme forme de la representation spectrale. La fonction H(p 2 ) 
etant en general divergente, cette serie ne converge pas et il faut renormaliser 
la masse Mq par une constante infinie qui absorbe la divergence de H(p 2 ). Cette 
procedure faite, la serie converge pour un couplage faible et loin des poles eventuels 
de G(p 2 ). Soit M la masse renormalisee et f(p 2 ) la partie restante de n(p 2 ) apres 
renormalisation de la masse (son expression exacte depend du point de renormali- 
sation). Le propagateur complet devient 



La fonction f(p 2 ) acquiert une partie imaginaire pour p 2 > M 2 euil , qui est l'energie- 
impulsion minimale pour creer des particules intermediaries reelles (c'est-a-dire sur 
leur couche de masse). En prolongeant G(p 2 ) dans le plan complexe en G(z), on 
retrouve bien la coupure sur l'axe des reels, pour p 2 > M 2 euil , presente dans la 
representation spectrale. 

Remarquons enfin que le pole du propagateur complet en p 2 = m 2 apparait dans 
la deuxieme representation spectrale comme un zero du denominateur, c'est-a-dire 
comme solution de 



On va voir immediatement que ce pole peut etre interprets comme la masse physique 
de l'etat a une particule cree par le champ. Cette masse physique m n'est pas 
necessairement identique a la masse renormalisee M. 

2.3 Particule stable : representation spectrale 

Supposons que le champ cree un etat a une particule de masse m, ainsi que des 
etats a plusieurs particules. La fonction spectrale (|2.3j ) s'ecrit dans ce cas 





p 2_M 2 -f(p 2 )=0. 




+ J2W5^(p-Pn) | <0 |0(O)| Pn > 



n>2 



ou encore 




(2.11) 
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ou Z = | <0 |0(O)| q>\ 2 est une constante dependant uniquement de q 2 = m 2 par 
invariance de Lorentz. Elle peut s'interpreter comme la probability que le champ 
cree un etat a une particule a partir du vide. La fonction p mu iti(s) est une fonction 
continue et represente la contribution des etats a plusieurs particules. Ces etats ne 
contribuent que si l'energie disponible est superieure au seuil de production de ces 
etats sur leur couche de masse. Le seuil le plus bas est note M, 



2 

euil' 



Pmuiu(s) = pour s < M 2 euil . 

M seu u doit etre superieur a m sinon la particule de masse m se desintegrerait dans 
les particules associees a M seu u. Dans ce cas, cette particule n'apparaitrait pas dans 
les etats asymptotiques et la representation spectrale ne serait plus valide. On verra 
plus tard comment modifier cette representation pour y inclure la description des 
particules instables. 

En y inserant 1 'expression de la fonction spectrale ( |2.11| ), le propagateur complet 
(g2f ) devient 

G( P 2 ) = - iz 2 . + r d s^^- 

p z — m z + te Jm 2 ., v — s + te 

r 1 semi " 1 

= G lpart (p 2 ) + G multl (p 2 ). (2.12) 

Le pole du propagateur complet s'interprete done comme la masse physique de 
l'etat a une particule cree par le champ. La constante Z est le residu au pole du 
propagateur complet. 



2.4 Particule stable : propagation spatiale 



Rappelons que le but est d'etudier l'amplitude de propagation (|2.2|) , ce qui 
revient a analyser la dependance spatiale du propagateur en calculant sa transformee 
de Fourier par rapport a la tri-impulsion. Celle-ci sera notee G(p°, L) : 

^e^G(p 2 )= f- P 
(2tt) 3 [P ] J (2v 

A partir du propagateur (|2.12|) d'un champ creant une particule stable, on obtient 



G(p°,L) = J ^r 3 e^G( P 2 ) = J 7^e^ L (G lpart (p 2 ) + G multl {p 2 



r (r) o j, f d 3 p fp . L i_Z 

^lpart{P , Li) — / , . e -= 5 . 

J {Zn) 6 p l — m z + ie 



2.4. PARTICULE STABLE : PROPAGATION SPATIALE 
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Jo Jo Jo {ZTY)* p° — p 2 — m 2 + ie 



|p| d\p 



P^lPl i\p\ Lease < z 

(2vr) 3 
gi |p|i _ g— * IpI i 



47r 2 L 7o p° 2 — p 2 — m 2 + ie 

Z r°° e iqL 
qdq ■ 



47r 2 L j-oo p° 2 — q 2 — m 2 + ie 

~' e iPmL , (2.13) 



4ttL 



ou L = |L| et p m = yp° 2 — m 2 + ie. La derniere etape a ete calculee par integrate 
de contour en fermant le contour par le haut. 



Si p° 2 — m 2 < 0, p m = i ym 2 — p° 2 et l'amplitude decroit exponentiellement en 
fonction de la distance. 

Si p° 2 — m 2 > 0, la decroissance est en L~ l qui est le facteur geometrique 
attendu, de sorte que la probability de detection sur un element de sphere de rayon 
L decroisse en L~ 2 . 

De meme, 

G mulu {p°, L) = - — - / ds Pmulu (s) e^ L , (2.14) 
An L Jm 2 ., 



ou p s = \Jp° 2 — s + ie. 



En rassemblant les deux termes (|2.13j ) et ( [2.14|) , on obtient 
G(p°, L) = -p- (ze ip - L + r dsp multl (s) e^ L 

4:71 L V J M ., 



,s e u 1 1 



Par le lemme de Riemann, le deuxieme terme tend vers zero lorsque la distance L 
tend vers l'infini. En effet, on peut effectuer la decomposition suivante : 



q2 

ip 3 L 1_ [°° f^r-V^L 



Gwtt(p°,L) = — — - / ds p mu i ti (s) e WsL - —— / dsp mu i ti {s)e 
An L Jm 2 ., 47r L J P o 2 

= Gmulti,l(j> , L) + G mu iuji(p , L) . 

Notons que G mvl u,i = si p° 2 < M 2 euil puisque p mu iti{s) = pour s < M 2 euil . Sous 
le changement de variable s — > uj = p s , G mu iuj devient 



G mulu AP°, L) = - j^r f du f(u) e^ L , 
47T L Jo 
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oil f(u) = 2up mu i ti (p° - uj 2 ) et p seuU = yp° 2 - M 2 euil . Une integrale par parties 
donne 



KAttLJ XL/ Jo duj 
oil Ton a utilise le fait que f(p seu u) = 2p seui ip mu i ti (M 2 euil ) = 0. De meme, sous le 



changement de variable s — > a = y s — p° 2 , G mu iti,ii devient 

G mu iti,ii(p ,L) = — — - / dag(a)e~ aL , 

47T L Jo 

ou gr(a) = 2ap mu iu(a 2 + p° 2 ). Une integrale par parties donne 

1\ /•» dg{a)^ aL 



ou Ton a utilise le fait que g(0) = et que p mu iu(s) ne diverge pas a l'mfini. En 
repetant le processus J fois, on obtient 

Ce precede est repetable autant de fois que les fonctions f(u) et p(a) (c'est-a-dire 
la fonction p mu iti) sont derivables. La decroissance en L est done plus rapide que 
n'importe quelle puissance inverse de L si ces fonctions sont infiniment derivables. 
Mais p m uiti n'est pas infiniment derivable aux seuils s N , qui sont les points ou les 
etats a N particules commencent a contribuer. La valeur de J correspondant au 
seuil s N se calcule en observant que pres de ce seuil ||, 

Pmulti(s) ~ (y/s - Sjv) 

En developpant p s autour de ce seuil, on calcule que la contribution a G mu i ti (p°, L) 
due au seuil s N a un comportement asymptotique pour L grand en 

pan h1~ ,; t 

G Npartlc (p , L) ~ — — (jmL)s W-*> exp(ip SN L) , (2.15) 

oil 



Ps N = \IP° 2 - s N + ie. 

Si s N > p° 2 , la racine choisie a une partie imaginaire positive. La masse p carac- 
terisant l'echelle d'energie a ete introduite pour raison dimensionnelle. La contribu- 
tion des seuils n'est importante qu'en dega de pL « 1, ce qui pour p « 500 MeV 
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correspond al~ 10 -16 m. Ce terme est done negligeable par rapport a Gi part (p°, L) 
lorsque Ton etudie la propagation sur des distances macroscopiques. 

Retournant a la formule (|2.2|) de l'amplitude de propagation et y inserant l'ex- 
pression fl2.13|) , on conclut que la propagation d'une particule stable sur une distance 
macroscopique L est donnee par 



ou p m = \JE 2 — m 2 et E = E Pi —E Pf > m. 

2.5 Particule instable : representation spectrale 

On a vu plus haut que les poles du propagateur G{p 2 ) correspondent a des etats 
stables a une particule. II est indispensable pour la stabilite de la particule qu'ils 
soient inferieurs au seuil de contribution des etats a plusieurs particules. Pourrait-on 
decrire les particules instables en violant cette condition ? Notons tout de suite que 
ce pole ne peut etre reel sinon la propagation spatiale serait du meme type que celle 
de la particule stable, alors que Ton s'attend a une decroissance de l'amplitude de 
propagation, due a la disintegration possible de la particule. Supposons done que 
le propagateur G(z) ait un pole complexe z dont la partie reelle est superieure a 
M seui i. Pour des valeurs reelles de z = p 2 , le propagateur ne presentera pas de pole 
mais un pic en p 2 = Tie z d'autant plus eleve que la partie imaginaire du pole est 
petite. II est tentant d'interpreter la partie reelle du pole comme la masse de la 
particule instable tandis que la partie imaginaire devrait etre reliee au temps de vie 
(puisque son annulation implique un temps de vie infini). 

II reste un probleme : comme prouve plus haut, la representation spectrale du 
propagateur n'admet pas de poles complexes ! Neanmoins, la coupure sur l'axe des 
reels permet un prolongement analytique de G(z) a travers la coupure sur d'autres 
feuilles de Riemann [§]. II suffit que le pole se trouve sur une de ses feuilles de 
Riemann supplementaires et la representation spectrale n'est plus violeeQ 

2 Toutefois, la correspondence univoque entre les particules instables et les poles complexes 




-iZ 



L 



(2.16) 
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Partons de la deuxieme representation spectrale (equation (PUD) du propagateur 
qui s'ecrit 

G{z) = z-Ml- n(z) 

et rappelons que (equation (|2.8|)) 



Xmli{p ) = lim TmH(z) = —rep a(p ) < 0. 

z—>p 2 +it 

Si la particule n'est pas extremement instable, la partie imaginaire du pole sera 
petite, c'est-a-dire que la partie imaginaire de la solution de z — Mq — Il(z) sera 
petite et done negative en vertu de la propriete Xmli{p 2 ) < 0. On va done effectuer 
le prolongement de G(z) a travers la coupure vers le has ce qui facilite les choses 
puisque le propagateur est connu juste au-dessus de la coupure. 

Le prolongement analytique de G(z) est defini tel que la fonction analytique- 
ment continuee Gn(z) juste en dessous de l'axe reel soit egale a la fonction originale 
G(z) juste au-dessus de l'axe reel. On prolonge done de fagon analogue H(z) en 
Hii(z) et Gu(z) s'ecrit 

GII{Z) = z-MS-IL n (z)- 
Le pole Zq est solution de l'equation 

z - Ml - n n (zo) = . 

Le propagateur complet renormalise peut aussi etre prolonge analytiquement et 
l'equation ( |2.10| ) devient 



^)= z _ M :_ }n{z) . (2.17) 

ou M est la masse renormalisee et fu(z) le prolongement analytique de la fonction 
f(p 2 ) qui est la partie restante de l'energie propre apres renormalisation de la masse 
(voir equation (|2.10|) ). 



de la matrice S (qui apparaissent ici dans le propagateur) n'est pas demontree dans l'approche 
axiomatique de la matrice S. II est possible jlj de construire un modele ou une resonance de type 
Brcit-Wigner ne correspond a aucun pole et vice versa. Par exemple, on pourrait appliquer ce 
modele a la resonance scalaire /o(400 — 1200), si Ton ne veut pas l'inserer dans l'octet scalaire du 
groupe SU(3) /. 
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On aimerait disposer d'une expression equivalente ou le pole complexe du propa- 
gateur apparait explicitement au denominateur. Dans ce but, l'energie propre est 
developpee autour du pole z : 

U n {z) = Ujj(z ) + (z- z ) W n {z ) + h(z) (1 - n'^o)) , 

avec h(zo) = et h'(z ) = 0. Les primes indiquent les derivees par rapport a z. 
Le denominateur du propagateur se reecrit 

z-M^-Un(z) = (l-n^(zo)) (z-z -h(z)) 
= Z^ 1 (z — m 2 + imT — h(z)^j , 

ou 

Z' 1 = 1-Yl h (zq) et ZQ=m 2 -imT. (2.18) 
Avec ces definitions, le propagateur complet peut se reecrire pres du pole comme 



z — m 2 + imT — h(z) 

La constante Z, appelee constante de renormalisation de la fonction d'onde, est le 
residu du propagateur complet au pole. Bien que le residu du propagateur de la 
particule stable ait pu etre interprete comme la probabilite de creer un etat a une 
particule, ce n'est pas le cas pour la particule instable puisque la constante Z est 
complexe. 



2.6 Particule instable : propagation spatiale 

Le modele simplifie ne convient pas tout a fait a la description de la propagation 
d'une particule instable en raison de 1' approximation d'etat stationnaire. II sera 
remedie a cet inconvenient dans le prochain chapitre. Cependant, cette inconsistance 
n'invalide pas les legons que Ton peut tirer de l'analyse de l'amplitude Q2.2p du 
modele simplifie qui, rappelons-le, est proportionnelle a la transformee de Fourier 
du propagateur complet. Bien sur, ce n'est que du point de vue du modele plus 
sophistique que Ton pourra en juger. 



40 



CHAPITRE 2. LE PROPAGATEUR COMPLET 



L'etude de la propagation spatiale ressemble a celle de la particule stable. Avec 
l'objectif d'utiliser le theoreme des residus, on va reecrire le propagateur comme une 
somme d'une fonction singuliere au pole mais analytique ailleurs et une fonction 
reguliere au pole mais non analytique : 

Gn{z) = Gi part {z) + G mu iti{z) 



ou 



Gl p art{ z ) 
G multi (z) 



iZ 



z — m 2 + imT 

i Z i Z 



z — m 2 + imT — h(z) z — m 2 + imT 

G mu iti(z) est une fonction reguliere au pole z = m 2 — imT mais non analytique a 
cause des seuils ou les etats a plusieurs particules commencent a contribuer : pour 
un seuil donne, il existe un J tel que la Jeme derivee de G mu iu(z) est discontinue. 
La transformee de Fourier du propagateur par rapport a sa tri-impulsion est 

En faisant les memes calculs que pour la particule stable (equations (|2.13|) et ( [2.151) ), 
on obtient 

G(p°, L) = ~r—f e^ p ° 2 - m2+imri + V d N (n N L)-l , (2.19) 

AixL N > 2 L 

ou les definitions de s N , p SN et L sont les memes que pour la particule stable et les 
d N sont des coefficients qui ne peuvent etre calcules que dans un modele particulier. 
/ijv correspond a /i pour l'etat N. 

On a utilise l'hypothese que la particule n'est pas extremement instable pour 
negliger les termes en T/m, sauf dans l'exponentielle ou ils sont multiplies par L 
et contribuent notablement a grande distance. Sinon, il n'y a de toute fagon pas de 
propagation macroscopique a observer car la particule se desintegre trop vite. Par 
contre, si la particule a un long temps de vie, le pole se situe pres de l'axe reel et 
T/m < I. 

Le premier terme peut etre identifie comme l'amplitude de propagation de la 
particule instable : on a bien une decroissance exponentielle avec la distance. Sous 
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quelles conditions ce terme est-il dominant par rapport a la contribution des etats 
a plusieurs particules? Sans perte de generalite, posons [i N = \p°\. 

1. Pour m 2 — p° 2 ^> mT, ce terme a une decroissance exponentielle tres rapide 



selon la distance en exp(— y m 2 — p° 2 L) et la propagation spatiale macro- 
scopique de la particule instable est inobservable : G mu m domine G\ part . 

2. Pour \p° 2 — m 2 \ « 0{mT) ou plus petit, la decroissance de ce terme est expo- 
nentielle en exp(-aVmTL), oil a est un nombre d'ordre un. La decroissance 
est plus lente que dans le premier cas. Ce terme domine G mu ui sur des distances 
intermediaires (voir cas suivant) mais Interpretation de V comme l'inverse du 
temps de vie est ici incorrecte. Ce domaine est cependant tres petit dans la 
majorite des cas puisque T/m <C 1. 

3. Pour p° 2 — m 2 ^> mT, on peut developper l'argument de l'exponentielle en 
mr/p m , oil p m = \fp° 2 — m 2 . La decroissance selon la distance est une expo- 
nentielle plus moderee : 

Gi part (p°,L) = — exp [ip m L - —Lj . 

Comme dans le cas d'une particule stable, ce terme sera domine par G mv xu 
pour des distances petites (/j,L < 1), mais il sera aussi domine a grande dis- 
tance en raison de sa decroissance exponentielle, alors que G mu iu ne decroit 
qu'en puissance de L. Pour des distances intermediaires, G\ part est dominant. 
Si Ton utilise la relation L = vT valable pour une trajectoire classique, avec 
v = p m /\p°\, on retrouve la forme habituelle de la decroissance de l'ampli- 
tude de propagation en fonction du temps de propagation T, mais sous forme 
relativiste : 

G lpart ( P ,L) ~ exp \-—L) ~ exp [-^T) - exp ^ , 

oil r est le temps propre de la particule. La largeur T s'interprete dans cas 
comme l'inverse du temps de vie ||. 

La contribution de G mu i t i est-elle observable a grande distance ? Considerons la 
contribution du seuil de production de deux particules. G mu iu commence a dominer 
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G lpart lorsque 

ClGXP ( _ 2^ L ) - C,2(b ° |Lr§ ' 

ou C\ et C 2 regroupent des constantes. 

Prenant le logarithme de cette expression, les facteurs constants C\ et C2 auront 
une contribution negligeable pour L grand et l'on obtient la condition 

^>31n(|p°|L). 

Pm 

Pour des particules se desintegrant faiblement (ou quasi- stables) , on peut se restrein- 
dre au domaine spatial 

r _ r m 

W\ <<TL<< r 

car pour TL m m/F m O(10 lA ) (kaons neutres), G\ vart et G mu m seront tous les deux 
inobservables. La borne inferieure du domaine vient de la condition fiL ^> 1. Des 
lors 

31n(|p°|L) = 31n^r^ =31n^ . 
Gmuiti domine G lpart pour 

rL>3^1n(|p°|/r)- 

m 

Cette borne depend de la valeur de p m : moins la particule va vite, moins elle va loin ! 
II est plus significatif de passer au temps de propagation par L = vT. La condition 
devient 

TT > 3 ^ln(|p°|/r) > 31n(M/r). 
m 

Si M/F ss C(10 14 ), FT > 97, ce qui est inobservable. 

La contribution des etats a plusieurs particules est done inobservable pour les 
particules quasi-stables. Le regime ou la propagation de la particule est bien decrite 
par le propagateur Gi part est fixe par les conditions 



HL > 1, 

FT < 31n(M/r) 
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En conclusion, l'amplitude de propagation ( |2.2| ) d'une particule quasi-stable sur 
une distance macroscopique L est donnee en tres bonne approximation par 



A ~ M P M D 2-k5{E Di -E Df + E Pi -E Pf ) — exp f ip m L - —L \ , 



(2.20) 



on 



p m = VE 2 - m 2 et E = E Pl -E Pp > 



m. 



2.7 Exemple : propagateur complet du kaon 

II sera utile par la suite de disposer de formules explicites dans un cas concret. 
Considerons le kaon neutre, de masse M , en interaction avec des pions charges de 
masse m (M > 2m). Le lagrangien d'interaction s'ecrit 

C = —gKmr* . 

Le calcul a une boucle de l'energie propre du kaon donne, pour s quelconque 

-i u(s) = -i 5M> - ^4 ( i - *< y i n ( i - 4 ^)? +i , 

4vr 2 V s J (l-4m2/s)5 - 1 
ou 5M 2 est une constante divergente qui vaut en regularisation dimensionnelle 

2 m 2 
SM 2 = ln47r + 7 + In -f - 2 , 

Dans cette expression, e = 2 — d/2, ji p est le point de renormalisation et 7 la 
constante d'Euler. SM 2 est incorporee dans la masse renormalisee : 

M 2 = M 2 + 5M 2 . 

Le propagateur renormalise ( [2.10|) s'ecrit 



G(s) = lim 



z-^s+ie Z - M 2 - f(z) 

avec 

4?r 2 \ z ) (l-4mj/ 2 )i-l 
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Pour z reel, la fonction f(z) a une partie imaginaire non nulle si z est superieur a 
Am 2 et on a une coupure sur l'axe reel pour s > Am 2 : 

f{s±ie) = u(s) Tiv(s) , 



avec 



47T 

La condition (12.81) est satisfaite : 



= 4f 1 -i^y ln ( i - 4 </^ +i 

4tt 2 V s / (l-4m|/s)s - 1 



Imn(s) = fl - 2 0(s - Ami) < . (2.21) 

Arc \ s J 

La fonction spectrale peut se calculer explicitement a partir de l'equation ( p.5| ) : 

pis + ie) = —Im(iG(s + ie)) 

— Ttti 

7T s- M 2 - f(s + ie) 
= I V M (2 22) 

7T ( S - M 2 -U(s)f +V 2 ( S ) ' 



Selon l'equation ( p,17|) , le prolongement analytique sur la deuxieme feuille de Rie- 



mann est donne par 

G n (z 



z — M 2 — u(z) + i v(z) 
Le pole zq est solution de 

Zq — M 2 — u(z ) + i v(zq) = . 

Pour un couplage g petit, u(z ) et v(z ) sont petits par rapport a M 2 et z = M 2 . 
Dans ce cas, l'expression du pole au premier ordre en g 2 est donnee par 

z = M 2 + u{M 2 ) - i v{M 2 ) + 0{g 2 ) =m 2 K -im K T . 
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On a done 



m\ = M 2 + u(M 2 ) + 0(g 2 ) , 
m K T = v(M 2 ) + 0(g 2 )=v(m 2 K ) + 0(g 2 ) 



e'est-a-dire 



a 2 / \± 

<m 2 K - Ami) 2 • ( 2 -23) 



4irm 2 K 



On pourrait aussi calculer la constante Z a l'aide de l'equation ( 2.18|) . 



En ne tenant en compte que du seuil a deux pions, on peut ecrire 

I,n n(. +<*) = -v (.) = (^j) 4 S (. - . 

Remarquons que si Ton veut seulement connaitre XmU(s), il est plus facile de 
le calculer a l'aide des regies de Veltman-Cutkosky Je trouve cependant 

interessant de montrer comment tous les parametres de la representation spectrale 
peuvent etre evalues dans un exemple concret. 
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Chapitre 3 

Melanges de propagateurs 



Au chapitre [I], nous avons defini en mecanique quantique la notion de melange 
de particules dans le cas des oscillations par l'impossibilite de faire coincider la base 
d'etats en interaction (etats propres de saveur) et la base d'etats propres de masse. 
Les definitions de melange et d'oscillation en theorie des champs sont analogues, 
excepte le fait que ce ne sont plus les etats physiques qui sont melanges mais les 
champs. Cette distinction permettra d'eviter tous les ennuis lies a la definition de 
bases d'etats propres de saveur et de propagation. 

Le lagrangien total est subdivise en un lagrangien L propag , decrivant la propa- 
gation des particules, et un lagrangien L int , decrivant les interactions produisant ces 
particules. Ces deux parties du lagrangien total peuvent etre distinguees s'il existe 
une transformation dite de saveur, laissant Li nt invariant mais modifiant L propag , a 
laquelle on associe le nombre quantique de saveur. II y a melange de particules, si le 
propagateur construit a partir de L propag et representant la creation d'une particule 
de saveur a au point x et l'annihilation d'une particule de saveur ft au point x' 
est non diagonal, c'est-a-dire non nul pour a ^ ft. Le lagrangien L propag contient 
toujours le terme cinetique et le terme de masse. Dans le cas de particules instables, 
on y incorpore aussi l'interaction provoquant la disintegration. 

Le type de lagrangien utilise depend du systeme de particules etudie. La cons- 
truction du propagateur la plus simple correspond toujours a des particules stables 
dont on connait le lagrangien fondamental. Ce cas est illustre par le systeme compose 
des differents types de neutrinos. Le lagrangien L propag contient la matrice de masse 
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(generee par les interactions de Yukawa) tandis que le lagrangien L int contient les 
interactions faibles (groupe SU(2) L x U(1) Y )- Si Ton se place dans la base de saveur 
des champs, le propagateur du systeme des particules melangees sera non diagonal. 
Comme on va le voir dans ce chapitre, il peut etre relie au propagateur diagonal de 
la base de masse par une transformation unitaire. 

Une construction plus compliquee du propagateur est necessaire si les particules 
auxquelles on s'interesse sont instables. II se peut que l'on connaisse le lagrangien 
fondamental, comme cela pourrait se produire si les neutrinos s'averent instables 
et que l'on connait les interactions a l'origine de la disintegration. Dans ce cas, il 
faut incorporer dans L propag les interactions a l'origine de l'instabilite. Le propaga- 
teur libre est alors remplace par le propagateur complet obtenu par sommation sur 
l'energie propre comme nous l'avons vu au chapitre precedent. Ce propagateur com- 
plet est non diagonal et doit etre diagonalise par une transformation qui cette fois 
n'est pas necessairement unitaire, puisque la matrice de l'energie propre du systeme 
n'est pas hermitienne quand les particules sont instables. 

Une derniere complication surgit si le lagrangien du systeme n'est pas complete- 
ment connu, comme c'est le cas pour le systeme des kaons neutres ou des mesons B. 
Nous sommes en effet incapables de deriver le lagrangien de ces particules a partir du 
lagrangien fondamental de la QCD et des interactions faibles. Neanmoins, il est pos- 
sible de construire un lagrangien effectif decrivant ces particules en employant comme 
contraintes les differentes symetries du lagrangien fondamental de la QCD et des in- 
teractions faibles ainsi que la facon dont elles sont eventuellement brisees. Le terme 
lagrangien effectif signifie que ce lagrangien ne contient que les particules observables 
a basse energie que sont les mesons pseudoscalaires et eventuellement les mesons vec- 
toriels. Le lagrangien effectif contient un nombre infini de termes qui sont classes 
en un developpement perturbatif en l'energie du processus. On parle de theorie chi- 
rale perturbative. Chaque terme de ce lagrangien represente une sommation sur un 
nombre infini de diagrammes de Feynman de la theorie fondamentale sous-jacente, 
chaque diagramme respectant la structure du terme en question. Les differentes 
interactions fondamentales sont done melangees dans chaque terme du lagrangien 
effectif. On definit la partie forte du lagrangien effectif comme l'ensemble des termes 
fournissant les contributions dominantes aux processus respectant les symetries de 
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la QCD. Le terme de masse est le terme non derivatif quadratique en les champs. 
II est genere au niveau fondamental a la fois par les interactions electrofaibles et les 
interactions de la QCD. Dans la meme logique, la partie faible du lagrangien effect if 
est definie comme l'ensemble des termes generant les processus violant les symetries 
de la QCD et de l'electromagnetisme mais respectant les symetries des interactions 
faibles au niveau fondamental. La partie electromagnetique du lagrangien effectif 
est definie de fagon analogue. Les parametres de ce lagrangien effectif sont fixes par 
des resultats experimentaux car leur prediction theorique n'est pas tres satisfaisante 
en raison de la difficulte a evaluer les elements de matrice hadroniques ainsi que les 
corrections a longue distance ||57| . Le traitement du systeme B°B° se fait a une autre 
echelle d'energie ou la theorie chirale perturbative n'est plus valable. Neanmoins, la 
partie a courte distance de la machinerie utilisee pour calculer les parametres du 
systeme K°K° (operateurs effectifs a quatre quarks et coefficients de Wilson) est 
immediatement transposable aux B°B°. 

Dans le but de decrire les oscillations des kaons, le lagrangien effectif peut 
egalement etre separe en un lagrangien Li nt contenant la partie forte du lagrangien 
effectif et un lagrangien L propag contenant le terme cinetique, le terme de masse 
(contenant des masses degenerees puisque le lagrangien decrit une particule et 
son antiparticule) ainsi que la partie faible du lagrangien effectif a l'origine de la 
disintegration des kaons. De nouveau, le propagateur complet est construit par 
sommation sur l'energie propre. Une difference de masse apparait alors entre les 
kaons se propageant et se manifeste par la non degenerescence des parties reelles des 
poles du propagateur complet du systeme des kaons. Cette difference de masse peut 
etre attribute aux interactions faibles (mecanisme GIM inclus!) avec le bon ordre 



de grandeur [58 



Tout d'abord, nous expliquerons comment diagonaliser le propagateur de par- 
ticules stables si le lagrangien est connu. Ensuite nous montrerons comment calculer 
le propagateur complet pour des particules instables en melange et comment le dia- 
gonaliser. Le systeme K°K° sera examine comme exemple. Enfin, la probability 
d'oscillation sera calculee dans le cadre du modele simplifie presente au chapitre 
precedent et les deficiences du modele seront identifiees. 
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3.1 Propagateur de particules stables en melange 

Dans cette section, nous expliquons comment diagonaliser le propagateur de 
particules scalaires stables si leur lagrangien est connu (le cas des fermions est discute 
a la section |B~6|) . Le lagrangien L propag contient uniquement le terme cinetique et le 
terme de masse. La matrice de masse est hermitienne et peut etre diagonalisee par 
une transformation unitaire V sur les champsf] : 



v, 



a 



£vj^, (3-1) 



J 



oil les indices grecs designent toujours la base de saveur et les indices latins la base 
de masse. 

Le propagateur est defini par la fonction a deux points ordonnee dans le temps : 

Gp a {x'-x) =<O\T(p (x') v a {x)) |0> . 

Comme la contraction de champs (theoreme de Wick) ne s'applique que sur des 
champs representant des particules de masses determinees, il faut substituer dans la 
fonction a deux points la relation ( |3.1| ) : 

G Pa {x' -x) = J2 Vl k V ja <0\T (v k {x') v*(x)) | 0> . 
La contraction des champs donne 

G/3a(x' - X) = Vfo G D,jj(x' - X) Vja , (3.2) 

j 

oil GDjj(x f — x) est le propagateur (D pour diagonal) d'une particule scalaire libre 
de masse rrij : 

G D jAx' -x)= I /4 e~^ x '^ . (3.3) 

' nK J J (2vr) 4 p 2 -m 2 j +ie v ; 



1 Remarquons que si l'on definissait des etats comme au chapitre |] (equation fll.lj )), la relation 
ci-dessus donnerait 1^0(0) >= J2j ^aj > > c'est-a-dire que la matrice V diagonalisant le terme 

de masse est la transposed de celle du chapitre |l|. 
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3.2 Propagateur de particules instables en melange 

Le propagateur complet de particules melangees est une matrice non diagonale. 
II s'obtient a partir du propagateur en l'absence de melange par sommation sur 
l'energie propre comme dans le cas des particules non melangees (equation (|2.9|)). 
Cette sommation s'effectue en passant par l'equation de Dyson pDfl , qui s'exprime 
sous forme diagrammatique par 

G(s)= — 





+ 




ou sous forme matricielle 

ou G l p est le propagateur de Feynman des particules libres (done non melangees) 
et —i Ii kl est l'energie propre du systeme des particules melangees dans la base de 
saveur. En sous-entendant la notation matricielle, l'equation precedente s'ecrit 

G = G F + G F {-iU)G 

ce qui donne pour un propagateur inversible 

iG- 1 = iG F 1 - U. 

Posons que le propagateur des particules libres est donne par la matrice 

iG F \s) = sI-Ml, 

ou / est la matrice unite et Mq la matrice de masse non renormalisee. Le propagateur 
complet devient 

iG-\s) = s/-M 2 -n(s). 
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Considerons le melange de deux particules scalaires, par exemple le systeme 
K°K°. Apres avoir effectue le prolongement analytique necessaire dans la deuxieme 
feuille de Riemann, les etats propres de propagation peuvent etre caracterises par 
les poles complexes du propagateur complet, ou encore par les zeros de l'inverse 
du propagateur. Nous allons diagonaliser ce propagateur pour etudier revolution 
spatio-temporelle des particules se propageant sur des distances et intervalles de 
temps macroscopiques. Separons l'analyse en trois etapes. 

1. La diagonalisation est toujours possible si les valeurs propres du propagateur 
complet sont distinctes deux a deux : 

-iG(z) = (zl - Ml -n(z))" 1 (3.4) 
= (zI-M"-f(z)y l 

, x / (z-m\- fi(z))- 1 \ 

= W W\ n , 2 , w ) W 

V (z - m 2 2 - f 2 (z)) J 

Pour z fixe, (z — — fj(z))^ 1 sont les valeurs propres de —iG(z). Chaque 
pole Zj est la solution de l'equation det G^i^Zj) = 0. La matrice de diagonali- 
sation W(z) contient les vecteurs propres correspondants. Si le melange entre 
les particules est faible, les parties reelles des poles du propagateur seront ap- 
proximativement egales aux masses avant calcul de l'energie propre et seront 
positives. Dans l'analyse de la propagation spatiale, la dependance en z des 
matrices W(z) engendrera des contributions a l'amplitude en puissances in- 
verses de la distance macroscopique de propagation L. Elles sont dues aux 
seuils de production d'etats a plusieurs particules et s'ajoutent aux contribu- 
tions de meme type provenant des seuils des energies propres diagonalisees 
fi(z) et f 2 (z). 

2. On aimerait que les matrices de diagonalisation soient independantes de z, 
pour pouvoir les sortir des integrates lors du calcul d'une amplitude. On pose 

-iG^^ui^-f 1 \)v + Q(z). 
V {z-z 2 ) L J 

Les residus Z\ )2 aux poles sont incorpores dans les matrices de diagonalisation. 
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U et V sont des matrices constantes a specifier, z 1)2 sont les poles du propa- 
gateur complet G(z) tandis que Q(z) est une matrice reguliere aux poles. 
Cette expression doit satisfaire a deux conditions : G~ 1 G = I et GG~ 1 = I, ou 
encore 













V + Q(z) iG-\z) 



Pour que ces equations soient satisfaites au voisinage des poles, il faut que les 
matrices U et V soient composees des vecteurs propres de Mq + H(z) evalues 
en z\ et z 2 . 
Plus precisement, si 

U = (ui u 2 ) et V 

alors 

(M 2 + f{z h2 )) u 1)2 =z h2 u 1>2 



vi, 2 (M + /(« li2 )J = z lj2 v 4 li2 

de sorte que Ton ait, pres du pole Zj, 

G- X {z) U ~ z-Zj et V ~ z - Zj , 

sinon le membre de gauche des equations a satisfaire serait divergent et ne 
pourrait etre egal a la matrice unite du membre de droite. Cette diagonalisation 
exacte est a l'origine des methodes utilisant deux angles de melange pour un 
melange a deux particules |5U . 



3. Un cas particulierement interessant est celui de poles quasiment degeneres. 
Des phenomenes d'oscillation de la probabilite de detection d'etats propres 
d'interaction se revelent alors dans revolution spatio-temporelle. Par exemple, 
les poles des particules K°K° sont degeneres si l'interaction faible est negligee. 
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A l'ordre g 2 en l'interaction faible, on peut approximer 

f(z) « f(m 2 ) , (3.5) 

ou m est la masse degeneree (c'est-a-dire la masse si les interactions faibles 
sont absentes). On a alors 

U = V- 1 et Q(z) = 

ou encore 




iG{z) ~ ir-i - - v 

(*-z 2 ) 



(3.6) 



Les residus Z\p sont incorpores a la matrice V. Cette approximation est ana- 
logue a celle qui consiste a prendre un hamiltonien effectif constant dans l'ap- 
proche non relativiste de Wigner-Weisskopf [|13|] : l'energie propre apparaissant 
dans le propagateur correspond a la matrice de masse effective de Wigner- 
Weisskopf si l'on approxime cette energie propre par une constante. 

En conclusion, les approximations ci-dessus menent a analyser revolution d'un 
melange de particules comme la superposition des evolutions des particules non 
melangees. Notons que les matrices de diagonalisation U^ 1 = V ne sont pas unitaires 
si le propagateur n'est pas hermitien, comme dans le cas du systeme K°K° lorsque 
la violation CP est prise en compte. 

Ce traitement peut etre immediatement etendu a des melanges de particules 
vectorielles, en particulier le systeme p°—uj—(f) qui a une grande importance dans la 
description des interactions en-dessous de 1 GeV (modele de dominance vectorielle) . 
Comme il s'agit de particules vectorielles, le propagateur contient un facteur de la 
forme — k ll k u /m 2 1 qui n'influence pas la diagonalisation puisque la contribution 
du terme k^k^/m 2 est toujours nulle en raison du couplage du p° a des courants 
conserves. Le melange des propagateurs a ete etudie dans |JIJ. Par contre, la com- 
posante longitudinale du propagateur doit etre prise en compte dans le melange 
photon-boson Z, qui est d'ailleurs complique par la question de l'invariance de jauge 



et de la brisure spontanee de la symetrie SU(2)l x U(l)y |59]. 
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3.3 Exemple : le systeme K°K° 

Quand il s'agit de melange de particules, le cas d'ecole est le systeme K°K° 
car il presente a la fois toutes les complications possibles et toutes les vertus qui 
les rendent observables : difference de masse et oscillation, violation CP directe et 
indirecte, quasi-stabilite done propagation macroscopique, meme ordre de grandeur 
de la longueur d'oscillation et de la longueur de disintegration. L'analyse qui suit 
s'applique sans modification aux systemes D°D° et B°B°. 

L'energie propre du systeme K°K° ne peut etre calculee que dans une theorie 
effective, dont un exemple est le lagrangien C = —gRim* qui apparait dans le calcul 
de K — > 7T7T vu a la section EO. On se limite ici a parametriser le propagateur. La 



symetrie CPT impose l'egalite des elements diagonaux dans la base de saveur f32|] . 
Une symetrie CP imposerait en plus l'egalite des elements non diagonaux mais nous 
savons qu'elle est violee dans ce systeme. La petitesse de la violation (~ O(10~ 3 )) 
impose neanmoins que les elements non diagonaux soient presque egaux. L'inverse 



du propagateur complet pour des kaons neutres d'impulsion p est parametrise [|| 
par 

iG-V) =( ( - l<5 " |A ' 0> ^^Uf d ° + b ) (3.7) 
\ (K^G-^K ) (K^G-^K ) J \a-b d J 

avec les definitions 

d = p 2 -m 2 - f 0Q (p 2 ) 
a + b = -foo(p 2 ) 
a-b = -foo(p 2 ) 

ou m est la masse renormalisee et les —ifapip 2 ) sont les energies propres complexes 
renormalisees du systeme des kaons neutres. G~ l est l'operateur correspondant au 
propagateur. Notons que les termes non diagonaux dependent de la convention de 
phase choisie pour les kaons. Nous y reviendrons au chapitre |[ 

La base des vecteurs propres de l'operateur CP est reliee a la base de saveur 

par 

' l^i) \ 1 I 1 lW \K°) 



\K 2 ) \/2 I 1 -1 M \K°) 



(3.8) 
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ou CP\K°) = \K°). Le K\ est done pair sous CP tandis que K 2 est impair sous CP. 
Dans cette base, l'inverse du propagateur s'ecrit 

On voit clairement dans cette base que la violation CP se traduit par le parametre 
b. 

Dans le but de diagonaliser ce propagateur, definissons un parametre complexe 
e tel que 





b 



1 + e 2 2a 
L'inverse du propagateur est diagonalise en 



(3.9) 




d + a^§ W 1 
d-a 1 ^ )\-e 



Par consequent, la base physique des kaons neutres consiste en deux etats K L S de 
masses m L S et largeurs T^s deiinies par 

d s = p 2 — mi, + im s Ys = d + a- 



d L = p 2 — m 2 + im L Ts = d — a- 



1-6 2 

1-e 2 



' ' " ° 1 + e 2 

ou l'on a fait l'approximation expliquee dans la section precedente d'evaluer l'energie 
propre a la masse degeneree (f(z) ~ f{m 2 )). Dans cette notation, le propagateur 
s'ecrit 

1 



-iG{p z 




-e 



Les matrices de diagonalisation V et V 1 de la base de saveur K°K° sont donnees 
parQ 



1 I 1 



2(1 - e 2 ) V -e 




r 



2 La normalisation de V est en fait arbitraire car les normalisations de V et V -1 se compensent, 
puisqu'elles apparaissent toujours par paire dans notre approche ou seuls les etats K° et K° sont 
observables. Nous choisissons ici des normalisations egales. 
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Notons que si e est purement imaginaire, la matrice V est unitaire. Dans ce cas, 



il est possible de choisir une convention de phase etrange (voir section |6.4.2| ) dans 
laquelle e = 0. Les etats K\ y2 sont alors etats propres de propagation. On verra aussi 
que dans ce cas il n'y a pas de violation CP dans les oscillations. 

Le propagateur dans la base de saveur est done diagonalise en tres bonne ap- 
proximation par 

-iGV)-v-M d f ° jr. 



d 



-i 

L 



Comme prevu, ces matrices sont non unitaires puisque le propagateur dans la base 
de saveur est non hermitien. Le lien avec le formalisme couramment utilise dans 
la litterature se voit en defmissant une double base physique orthogonale et nor- 
malised (voir equation ( |1.6|) ), e'est-a-dire des etats entrants \ket) et sortants (bra\ 
independants qui sont vecteurs propres a droite et a gauche du propagateur |32|, p3l 





amsi que 




c2 




vT- \ -< 

A moins que e ne soit purement imaginaire (e'est-a-dire qu'il n'y a pas de violation 
CP), ces deux bases ne sont pas hermitiennes conjuguees l'une de l'autre. Si Ton tient 
a travailler uniquement avec la base de vecteurs propres a droite \K StL ), la question 
de la normalisation est incontournable puisque les vecteurs de cette base ne sont pas 
orthogonaux a leurs hermitiens conjugues. On choisit couramment [pll de normaliser 



par 1/yl + |e| 2 au lieu de — e 2 comme ici, de sorte que le produit scalaire de 

\Ks) avec son hermitien conjugue soit normalise a 1 (on procede de meme fagon 
pour \Ki)). Neanmoins, \K$) n'est toujours pas orthogonal a l'hermitien conjugue 
de \Kl). Comme il a deja ete mentionne au chapitre |I], le choix de la normalisation 
influence la forme des expressions du genre M. {Kl — > 7T7t) =<7T7t |T| Kl>, ou T est 
la matrice de transition. On verra au chapitre |6| la maniere correcte de definir de telles 
amplitudes en theorie des champs. L'ambiguite de normalisation en e 2 est beaucoup 
plus petite que la violation CP indirecte (~ e) et n'a done pas pose de probleme 
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pour Interpretation des experiences, en tout cas jusqu'aux premieres mesures de la 
violation CP directe (~ e' ~ e 2 ). 

Les quantites m SL et T s ,l sont mesurees experimentalement tandis que les 
parametres a, b, m et f(m 2 ) peuvent en principe etre calcules theoriquement ||57f . 
Les relations entre les deux ensembles de parametres sont 

a= \ 1 - " 2 ( m i ~ m s - ^( m i r i - m s^s)) , (3.11) 
m 2 - f(m 2 ) = - (m 2 L + m 2 s - i(m L T L + m s T s )) , 
b = ^ - (m 2 L -m 2 - i(m L T L - m s T s )^j . 



3.4 Oscillations dans le modele simplifie 

La propagation de particules en melange va etre decrite par le processus de 



la section |2.1| legerement modifie pour tenir compte de la saveur. Le processus de 
propagation d'une particule u a de saveur a produite a la source en une particule v@ 
de saveur (3 identifiee dans le detecteur est symbolise par : 

Pl{q) ^ P F (k)+u a (p) 

\ 

Mp) + W) ^ D F {k>) 

On suppose que Ton peut identifier la saveur a de la particule intermediaire v pro- 
duite dans la region de (t P ,x P ) au moyen de l'etat sortant Pp(A;) et la saveur j3 
du v detecte dans la region de (t D ,x. D ) au moyen de l'etat sortant D F {k'). S'il est 
impossible d'identifier la saveur sortante (ex : K°,K° — > il suffit de sommer 

sur les differentes saveurs. On garde les hypotheses simplificatrices de representer 
Pi et Di par des etats stationnaires infiniment bien localises dans l'espace et Pp et 
Dp par des ondes planes. 

Soit G/3 a (x'—x) le propagateur complet symbolisant la propagation de particules 
de saveur a produites en x en particules de saveur (3 detectees en x' . Nous venons 
de voir (equations (ft.2|) et (|3.6|) ) que le propagateur peut etre diagonalise en tres 
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bonne approximation par des matrices V constantesj] : 

G Pa (p 2 ) = (V- 1 G D (p 2 ) V)f3 a = £ V p f G DJj V 3 



3 



L'amplitude du processus global est donnee par l'equation fl2.1|) de l'amplitude de 
propagation d'une particule isolee, ou le propagateur G{p 2 ) est remplace par le 
propagateur Gp a (p 2 ). Elle peut done s'ecrire comme une superposition lineaire d'am- 
plitudes de propagation d'etats propres de masse : 

A(a^(3,T,L) = Y,V pj 1 A 1 V 3a , 

3 

ou l'amplitude partielle Aj s'ecrit, apres integration sur x et x' (voir equation (|2.2|)), 
A 3 ~ M P M D 2ir5 (E Di -E Df +E Pi -E Pf ) J -0- e ip " L G Djjj (E, p) 

ou L = x D -x P et E = Epj — Ep F >0. Si Ton integre sur la tri-impulsion, l'amplitude 
partielle devient (voir equation ( |2.20|) ) 

Aj ~ M P M D 2tt 5 (E Di -E Dp +E Pi -E Pf ) ^ exp (ip m .L - ^^l) , 



OU p m] = yJE 2 - m 2 . 

Les oscillations entre les etats propres de masse i et j surgiront des termes 
d'interference AiA* dans la probabilite donnee par la norme au carre de l'amplitude : 

^..exp( i( p m ,- Pm) .)L-g-g). 

Le facteur d'oscillation pour le melange ij est done un cosinus d'argument 

Am?. 

(Pm x -p mj )L = -7r^ L 

^Pm 



ou Am.. = m 2 — m 2 et p m = \fE 2 — m 2 avec m une masse de reference choisie par 
exemple comme m = (m» + mj)/2. On a fait l'approximation |m.j — m| <p m puisque 
dans les cas ou une oscillation macroscopique est observable, \rrii — rrij\ <C p rn . 

3 Rappelons une fois de plus que les matrices V utilisees ici sont les transposees des matrices V 



diagonalisant les kets dans le traitement de mecanique quantique (voir equation (1.1)). 
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Si Ton definit la longueur d' oscillation L° sc pour le melange ij par 

rose = ^Pm 

on retrouve la formule d'oscillation standard (|1.10|) . 

Ce modele repond a plusieurs des objections soulevees au chapitre [l] (section 
( |1 . 1| ) ) , en considerant dans un cadre relativiste (objections 1 et 2) la particule comme 
un etat intermediate non observe directement (objections 3,4,5 et 6) et en engen- 
drant une dependance spatiale (objections 7 et 8). 

Tout n'est pas rose pour autant. Par exemple, il est tentant d'interpreter E et 
p mi comme l'energie et l'impulsion de l'etat propre de masse i. Cette interpretation 
implique curieusement que les different s etats propres de masse de la particule 
oscillante ont la meme energie mais des tri-impulsions differentes alors que dans 
la derivation en mecanique quantique on considere le plus souvent le cas inverse 
(energies differentes mais tri-impulsions egales), sauf dans la prescription de Lipkin 
p8| , f45| . L'origine de notre resultat est evidemment le choix de proceder a l'analyse 
spatiale du propagateur. L'analyse temporelle du propagateur j|] implique de pren- 
dre des energies differentes et des tri-impulsions egales mais n'est pas justifiee par 
un modele de theorie des champs. 

II ne faudrait pas en conclure pour autant que la question de l'egalite ou non 
de l'energie- impulsion des etats oscillants est reglee par notre modele simplifie. 
L 'approximation d'etats stationnaires impose qu'il n'y ait qu'une seule valeur de 
l'energie possible pour les etats oscillants, tandis que l'hypothese de localisation a 
une precision infinie de la source et du detecteur n'impose, par contre, aucune con- 
trainte sur les valeurs des tri-impulsions des etats oscillants, qui sont alors fixees par 
les poles des propagateurs. 

La solution de ce probleme necessite done le recours a un modele plus realiste, 
qui abandonne les hypotheses de stationnarite et de localisation infiniment precises. 
II sera du meme coup possible d'etudier les conditions d'observabilite des oscillations, 
alors que le modele simplifie n'en fournit aucune. Par exemple, les oscillations sont 
interdites si tous les etats entrants et sortants sont des ondes planes, car il n'y 
a alors qu'une seule valeur possible pour l'energie-impulsion des etats oscillants. 
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L'analyse des situations intermediaries devrait done permettre d'etudier 1' influence 
des conditions experimentales. 

Finalement, on voudrait decrire de fagon coherente la propagation d'une parti- 
cule instable. De plus, comme mentionne ci-dessus, les corrections negligees en T/m 
sont du meme ordre de grandeur que le rapport Am/m dans le systeme des kaons 
neutres. II faudrait s'assurer qu'elles n'aient pas d'influence sur le facteur d'oscilla- 
tion. Dans le prochain chapitre, un modele plus sophistique de propagation d'une 
seule particule est developpe dans le but de repondre a ces questions. 
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Chapitre 4 



Propagation : le modele 
sophistique 

Ce chapitre est consacre a l'etude d'un modele de propagation sophistique, te- 
nant compte dans la mesure du possible des conditions de production et de detection. 
Abandonnant la simplification extreme de considerer les particules entrantes et sor- 
tantes soit comme des ondes planes soit comme des etats stationnaires localises a 
une precision infinie, on modelisera ces etats de maniere plus realiste par des pa- 
quets d'ondes |^|, . La particule intermediaire est representee par son propagateur 
relativiste. 

Ce chapitre est plutot technique, debutant par quelques considerations sur les 
paquets d'ondes, suivies de la formulation de l'amplitude. II continue par l'analyse 
temporelle de celle-ci, avec une etude detaillee des contributions non exponentielles 
a la propagation. On verra par exemple que les contributions des etats a plusieurs 
particules vues au chapitre precedent sont souvent absentes pour les processus de 
propagation macroscopiques, car les seuils de production de ces etats se trouvent 
souvent en dehors de la region d'integration. Ces contributions sont remplacees par 
celles des seuils dus a la limitation dans l'espace des impulsions de la grandeur des 
paquets d'ondes entrants et sortants. Le chapitre se termine par l'analyse spatiale 
de l'amplitude et le calcul de la probabilite du processus, en ne conservant que la 
dependance spatiale. 

Les bases du calcul ont ete jetees par Jacob et Sachs ||. La notation a ete 
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definie au chapitre Le processus global de propagation etait symbolise par 



(tp,x P ) 



P F {k) + v(p) 



\ 



v{p) + Dj{q') 



(t D ,X D ) 



D F (k') 



Rappelons la signification des notations. Pj represente l'ensemble des particules ar- 
rivant dans la region de production centree autour du point (t P ,x P ), d'impulsion 
totale q tandis que Pp represente l'ensemble des particules issues de la region de 
production et d'impulsion totale k, a l'exception de la particule intermediaire v 
dont on etudie la propagation. Di, Dp et (t D ,x D ) ont des definitions similaires 
mais concernent le point de detection. Le point d'interaction a la production est 
note par x et le point d'interaction a la detection est note x' . On pose que l'etat 
intermediaire se propageant de x en x' a les nombres quantiques d'une particule, 
pas d'une antiparticule. On ne tiendra pas compte des instabilites eventuelles des 
particules exterieures. Seule l'instabilite de la particule intermediaire sera prise en 
compteQ Si elle se desintegre dans le processus etudie, alors Dj(q') est en realite un 
etat sortant. L'ecriture d'une amplitude tenant compte de l'instabilite des partic- 
ules entrantes et sortantes ne pose en elle-meme aucun probleme. II suffit de con- 
siderer un processus plus global ou toutes les particules instables sont vues comme 
intermediaires. Les etats initiaux et finaux sont stables. Les regies de Feynman 
permettent d'ecrire facilement l'amplitude correspondant a un tel processus en cas- 
cade. Par contre, revaluation des integrales est la plupart du temps un obstacle 
insurmontable. Cest pourquoi on se limite a considerer comme etat intermediaire 
la particule dont on veut etudier les proprietes. 

Toutes les particules sont supposees sans spin. Pour localiser l'interaction de 
production en (t P , x P ), les particules entrantes et sortantes au point x sont modelisees 
par des paquets d'ondes qui ne se recouvrent (dans l'espace de configuration) que 
dans la region de production centree autour de (t P , x P ). En vue d'etudier l'influence 

^ertaines tentatives ont ete entreprises pour tenir compte de l'instabilite de la source dans 
les conditions de production, d'une part dans le cadre de la mecanique quantique non relativiste 
[ pl[ p2[ , et d'autre part en theorie des champs |(|| 
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des facteurs de production, les paquets ont une certaine largeur dans l'espace de con- 
figuration. De plus, comme les experiences mesurent generalement aussi 1'impulsion, 
on peut supposer qu'ils sont bien localises aussi dans l'espace des impulsions au- 
tour de leurs impulsions moyennes. L 'interaction de detection est localisee de meme 
maniere. Notons qu'il n'y aurait aucune difficulte a remplacer certains paquets d'on- 
des par des etats lies, si la modelisation d'un processus particulier l'exige. 
A ce processus correspond le diagramme de Feynman suivant 



D/q')- 



P F (k) 



D 



Ce diagramme represente toute une serie de processus. Le notre se distingue par des 
conditions experiment ales telles que l'interaction en x precede toujours l'interaction 
en x' et telles qu'il y ait un transfert d'energie (positive) de x en x' . Ces situations 
seront selectionnees automatiquement par la forme choisie des paquets d'ondes, avec 
pour consequence que seule la particule v contribuera a la propagation, et non 
l'antiparticule v. 

On ne specifiera pas ici la forme des interactions de production et de detection 
mais il n'y aurait aucun probleme a le faire. Par exemple, si la particule intermediate 
est un neutrino, les methodes de detection sont de trois types J20| : 



vi + X 
vi + X 

Ui + l 



l + Y, 
vi + Y, 
Pi + l, 



ou / symbolise un lepton et X et Y ne contiennent pas de leptons. Pour calculer 
l'amplitude du troisieme processus [pi[ , il faut sommer sur les deux diagrammes 
representant l'interaction par courant neutre et par courant charge. Dans les cas ou 
des neutrinos sont produits dans le detecteur (deuxieme et troisieme processus), il 
faut sommer sur leurs differents etats propres de masse possibles, puisque ce sont les 
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etats physiques asymptotiques. Le neutrino final doit alors etre reexprime dans la 
base de masse : v\ = J2j ^ij v j-. ce Q u i ajoute une matrice de melange supplement aire 
dans l'amplitude. Cette matrice n'a cependant rien a voir avec les oscillations de la 
particule intermediaire. 



4.1 Paquets d'onde 

Un paquet d'ondes representant un etat donne | 0> de masse m s'ecrit dans 
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»= J[dk]<j){k)\k> 



ou 0(k) est la fonction d'onde dans l'espace des impulsions et | k> est un etat a une 
particule d'impulsion k dans la theorie en interaction. On utilise la notation 

[dk]= dk 1 



(2tt) 3 J2E{k) 



ou E(k) = Vk 2 + m 2 . La normalisation des etats libres est donnee par 

<k| p>= 2E{k) (2vr) 3 (5 (3) (p-k) . 

On a 

<6\6>=1 si / -^-\d)(k)\ 2 = 1. 



(x ' t) = /7^ 0(k)e ~ i£(k)t+?kX 



(2tt) 3 

Dans l'espace de configuration, la fonction d'onde s'ecrit 

dk 
(2^ 

de sorte qu'elle satisfasse a l'equation de Klein-Gordon. La relation reciproque est 
donnee par 

0(k) = y"rfx0(x,t) e !£(k,1 - ik ' x 

Si le paquet d'ondes represente une particule dotee d'une impulsion approxima- 
tivement connue, il a un seul maximum global, disons en k = K. Le maximum de 
0(x, t — 0) se situera en x = pour autant que 0(k) ait une symetrie centrale autour 
de K : 

MK+k') = MK-k') 



4.1. PAQUETS D'ONDE 
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Nous ferons l'hypothese que les paquets d'ondes ont une impulsion bien definie : 
leur maximum est bien localise et |0(k)| a approximativement une symetrie centrale. 
On va faire l'hypothese supplementaire que la phase de 0(k) a aussi une symetrie 
centrale, de sorte que le paquet soit centre en x = 0. Sous ces conditions, le paquet 
d'ondes est note 0(k, K). 

Les paquets centres en x ^ sont construits par translation. A l'aide de 
l'operateur de translation exp(iP-x), oui = (T, X), on construit des paquets d'ondes 
centres en X au temps T. Si le paquet d'onde dans l'espace des impulsions est donne 
par 

$(k, K, X, T) = 0(k, K) exp {-i k-X + iE{k)T) , (4.1) 

alors le paquet d'ondes dans l'espace de configuration 
~ r dk 

$(x, t, K, X, T) = J — 0(k, K) exp (i k- (x-X) - iE(k)(t-T)) , 

sera centre en X au temps T. 

Sans perte de generality, on va travailler avec une seule particule dans Pi(q), 
dans Pp(^), dans Dj(q') et dans Dp{k'). L'extension a un nombre plus grand de 
particules est immediate et ne fait que compliquer la notation. Les paquets d'ondes 
correspondant aux differentes particules sont construits de sorte que ceux concernant 
la production du v sont centres en x P au temps t P tandis que ceux concernant la 
detection du v sont centres en x c au temps t D . lis sont notes 

\Pi> = J[dq\$ Pl (q,Q,x P ,t P )|P 7 (q)> 

\P F > = |[rfk]$p F (k,K,x P ,t P )|P F (k)> 

\Di> = yW]$D, (q / ,Q / ,x D ,t D )| J D J (q')> 

\D F > = J[dk']<5> DF (k',K > ,x D ,t D )\D F (k')> . 

On retombe sur le modele simplifie des chapitres precedents en imposant 

1. la condition de stationnarite sur Pj et Di, c'est-a-dire que les energies E P] et 
Edj sont constantes. 

2. la localisation infiniment precise de Pj et Di, c'est-a-dire que 0p 7 (q, Q) et 
^D/Cq'jQ') son t des fonctions constantes. 
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3. que les etats P F et D F soient des ondes planes, c'est-a-dire que les 0p F (k, K) 
et 0D F (k',K') sont des fonctions delta : 

PF (k,K)~5( 3 )(k-K) et DF (k',K')~5 (3) (k'-K'). 

Les energies de ces etats sont fixees par la relation E 2 = p 2 + m 2 si l'on impose 
qu'ils satisfassent a l'equation de Klein-Gordon. 

4.2 Amplitude 

La formule generate de l'amplitude est donnee par 

A =< Pf, D f I T ^exp (-i J d 4 xH^ -1\P I ,D I >, 

oil ?ij est le lagrangien d'interaction de la particule v qui se propage et T l'operateur 
qui ordonne dans le temps. Soit g la constante de couplage du champ v avec les autres 
champs. Developpant cette amplitude a l'ordre g 2 , et y inserant les expressions des 
paquets d'ondes, on obtient 

A = J [dq] $ P/ J [dk] $ Pf J [dq'} $ Dl J [dk'\ 

^ D F Aondes planes 

(q,k,k',q') 

avec 

Aondes P ianes(q, k, q' , k') = J d 4 x M P (q, k) e'^'^ j d 4 x' M D (q', k!) e~ i ( q '~ k '}' x ' 

x G{x' - x) 

ou Mp(q,k) et MD(q',k') sont les amplitudes des processus de production et de 
detection ; la contraction des champs a donne le propagateur 

G(x' - x) = / ^ e"^-*) G(p 2 ) . (4.2) 

ou G(p 2 ) = i(p 2 — Mq + ie)' 1 est le propagateur de la particule libre dans l'espace 
des impulsions. Les particules exterieures sont sur leur couche de masse : 

q° = E Pl (q) = ^q 2 + m 2 Pi , 

et ainsi de suite. 
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On a suppose que la particule (p° > 0) se propage de x en x' et l'antiparticule 
(p° < 0) de x' en x. Si les interactions a la source et au detecteur sont telles que 
la particule se propage de x' en x et l'antiparticule de x en x', la contraction des 
champs v aurait donne le propagateur 

G(x'-x) = J^e ip - [X '~ X) G(P 2 ). (4.3) 

Le signe different de l'exponentielle selectionnerait le pole de l'antiparticule lors de 
l'integrale de contour ulterieure. 

Si l'objectif est de decrire des particules instables, il est necessaire que ce soit 
le propagateur complet qui apparaisse dans l'amplitude et non le propagateur de 
la particule libre. II suffit pour cela de developper la formule generale de l'ampli- 
tude a tous les ordres et de garder tous les diagrammes representant une insertion 
d'energie dans le propagateur de la particule intermediate. La sommation de tous 
ces diagrammes est une serie geometrique en l'energie propre (voir equation (|2.9| )). 
Elle fournit la meme formule que ci-dessus excepte le fait que G(p 2 ) symbolise main- 
tenant le propagateur complet (equation (|2.6|)) : 

G(p 2 



p 2 — ml — Il(p 2 ) + ie 
On procede ensuite a un changement de variable : 

x ^ x + xp et x' — >x' + xd, 

ou xp = (tp, Xp) et xb = (tp), X/)), ainsi qu'a des integrations sur les variables x et 
x' qui donnent des fonctions delta. On arrive a la formule suivante de l'amplitude 



A = J d 4 P ip(p)G(p 2 )e- ipixD - Xp) 



(4.4) 



ou la fonction-poids ip(p) est une integrale de recouvrement des paquets d'ondes 
entrants et sortants. Elle est definie par 

<P(P) = J [dq] 0Px (q, Q) / W <f>*P F (k, K) J [rfq'] <j> Dl (q', Q') f [dk'] ^ (k', K') 

x (2tt) 4 5 (4) (A; + k' -q-q')5 ii \p-q + k)M P {q,k)M D {q',k > ). (4.5) 
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Insistons sur le fait que ce sont bien les paquets d'ondes cj) independants de xp et xo 
qui figurent dans la fonction <p(p). Les fonctions delta y figurant imposent la conser- 
vation de l'energie et de l'impulsion a la production et a la detection de la particule 
intermediaire v. Notons aussi que cette fonction depend de p (et done de la direction 
de l'impulsion totale des etats entrants et sortants) et non de p 2 . La conservation de 
l'energie-impulsion introduit un lien entre p et les valeurs moyennes Q, K, Q', K' 
des impulsions des etats entrants et sortants de sorte que </?(p) depende a la fois de 
la grandeur et de la direction de p. L 'exception est le modele simplifie du chapitre 
puisque la tri-impulsion des etats stationnaires est tout a fait indeterminee. Dans 
ce cas, la fonction-poids vaut (p(p) ~ 5(p° — E) et est independante de p. 

La plupart des calculs en theorie des champs concernent des processus micro- 
scopiques ou les etats entrants et sortants peuvent etre approximes par des ondes 
planes. Dans ce cas, les fonctions 0p 7 , 4>p F , 4>Di et 4>d f sont toutes des fonctions 
delta et la fonction-poids <p(p) est egale a 

<p(p) = {2uf5 {i \K + K'-Q-Q') 5 {4 \p - Q + K) M P (Q,K) M D (Q',K') , 



ou Q = (J Q 2 + m Pi , QJ et ainsi de suite. L'integration dans l'amplitude (|4.4|) est 
alors immediate et Ton obtient 

A = (2n) A 5 {A) (K + K'-Q-Q')G ((Q - Kf) e -*(.Q-K)-(*i>-*p) . 

Ce type d'expression ne peut mener a des oscillations dans une superposition de 
plusieurs amplitudes puisque la phase de l'exponentielle ne depend pas de la masse 
de l'etat intermediaire. Par contre, dans le cas d'une amplitude modelisant des 
processus macroscopiques, la fonction-poids cp(p) n'est pas une fonction delta et 
l'integration dans l'amplitude n'est plus aussi simple. L'integration dans l'amplitude 



(Oj) ne peut alors se faire sans recourir a des methodes d' approximation. D'ailleurs, 
la forme exacte des paquets d'ondes est inconnue. Remarquons enfin que l'amplitude 
( |4.4j ) depend a la fois de la distance et du temps macroscopiques de propagation. 
II faudra integrer sur le temps pour se debarrasser de la dependance temporelle et 
obtenir ainsi une expression applicable aux experiences . 



4.3. ANALYSE TEMPORELLE DE L 'AMPLITUDE 
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4.3 Analyse temporelle de l'amplitude 

On va d'abord etudier la dependance temporelle de l'amplitude en integrant 
sur l'energie p°. Soit l'integrale 



Pour etudier la propagation macroscopique de la particule, il suffit d'examiner le 
comportement asymptotique T — > oo de cette integrale. Cela ne peut cependant se 
faire sans une connaissance minimale de la fonction-poids tp(p). On va supposer que 
le spectre d'energie des particules entrantes et sortantes couvre un domaine limite. 
Par consequent, les fonctions delta figurant dans l'expression de la fonction-poids 
limitent le domaine de p pour lequel la fonction-poids est non nulle : 



oil V est un domaine a support compact. Ceci implique qu'il existe des bornes Mi 
et M 2 pour lesquelles 



On a pose des conditions experiment ales telles que <p(p) ^ pour p° positif puisque 
la propagation macroscopique (T > 0) de x en x' signifie qu'il y a un transfert 
d'energie positive de x en x' . Comme l'amplitude a ete derivee telle que la particule 
transfere une energie p° de x en x' (et l'antiparticule une energie de — p° de x' en x), 
il s'agit bien de poser des conditions experiment ales telles que p°>0. On supposera 
que le domaine (Ei,E 2 ) n'inclut pas zero. 

Si ip{p) etait non nul pour E\ < p° < E 2 < 0, l'integrale de contour ulterieure 
donnerait zero car la condition T > de propagation macroscopique de x en x' 
impose de fermer le contour par le bas. Pour etudier la propagation macroscopique 
d'une antiparticule, il suffit de remplacer G{x' — x) (equation f4.2| ) par G{x' — x) 
(equation [4.3[ ) et de poser des conditions experiment ales telles que ip (p) ^ pour 
E l <p°<E 2 <0. 

Comme pour l'etude de la contribution des seuils au propagateur complet, c'est 
le comportement de l'integrand aux bornes du domaine qui determinera le com- 




(p(p)^0 pour 0<E 1 <p°<E 2 et peV, 



cp(p)^0 pour Ml <p 2 < Ml . 



72 



CHAPITRE 4. PROPAGATION : LE MODELE SOPHISTIQUE 



portement asymptotique de l'integrale. II sera utile par la suite de faire la distinc- 
tion entre les particules se desintegrant faiblement, appelees quasi- stables et celles se 
desintegrant fortement, appelees resonances. Soit zq le pole du propagateur complet. 
En utilisant la methode de Jacob et Sachs ||, on montre que, pour des particules 
stables et quasi-stables, I(T) est donne en tres bonne approximation pour T grand 
par 

I(T) k,Zix (z Q + p 2 ) 2 (p(z , p) exp (-i<Jz + p 2 T 

pour autant que les conditions experiment ales soient telles que 

Ml<Vez <Ml 

sinon la propagation macroscopique serait inobservable. En effet, les conditions expe- 
rimentales ne seraient pas reunies pour qu'une particule intermediaire se propage 
sur sa couche de masse, comme une particule reelle. La derivation detaillee de cette 
formule figure dans l'appendice du chapitre. Notons qu'une integrate de contour 
sur un contour de type demi-cercle ne converge pas pour la plupart des types de 
paquets d'ondes. Par exemple, une fonction-poids gaussienne diverge sur ce contour. 
Le contour propose par Jacob et Sachs permet d'integrer une classe beaucoup plus 
large de fonctions-poids ; entre autres, les gaussiennes y sont incluses. Ce point est 
neglige dans l'etude des oscillations de kaons par Sudarsky et al jpB] . 

Si Tmza <C Tie zq, le coefficient de variation de la phase de l'amplitude en 
fonction du temps est \/Tie z + P 2 - Interpretant cette expression comme l'energie E 
de la particule, on peut ecrire 

E 2 -p 2 = -Rez . (4.7) 

Comme la masse de la particule au carre est egale a la partie reelle du pole du 
propagateur complet, la particule est, dans ce sens-la, sur sa couche de masse. Quant 
a la partie imaginaire du pole (dans le cas ou elle est non nulle), elle fixe la rapidite 
de la decroissance exponentielle de l'amplitude et est proportionnelle a la largeur de 
la particule instable. 

Les corrections a cette formule sont en puissances inverses de T. Elles sont de 
deux types : 



(4.6) 
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1. Des corrections dues au spectre fini de l'energie des particules entrantes et sor- 
tantes. Elles sont done aussi liees a l'incertitude sur la localisation temporelle 
des interactions de production et detection. La grandeur de ces corrections est 
proportionnelle a (AMT)""" 1 , oil AM = Mj — m et n est un nombre positif. 
Ces corrections sont importantes pour des temps tres petits et sont en pratique 
inobservables dans la propagation des les particules stables et quasi-stables. 
Dans le cas des particules quasi-stables, elles dominent aussi l'exponentielle 
decroissante pour des temps tres grands mais a ce moment les deux termes 
sont trop faibles pour etre observables dans la propagation. 

2. Des corrections dues a un seuil de production z = b 2 de plusieurs particules 
reelles, si celui-ci se trouve dans la region ou le spectre de l'energie de la 
particule intermediaire est non nul. Ce sont les corrections que nous avons 
examine qualitativement au chapitre || (equations (|2.15|) et ( |2.19|) ). 



(a) Si la particule est stable, la grandeur de ces corrections est proportion- 
nelle a (mT)~ 3 / 2 g 2 / Q 2 , oil Q = m — b et g est la constante de couplage 
avec les particules produites au seuil. Ces corrections ne sont importantes 
que pour des temps tres petits et sont en pratique inobservables dans la 
propagation. 

(b) Si la particule est instable, la grandeur de ces corrections est propor- 
tionnelle a (QT)~ 3 / 2 T /Q (ou Q = m — b). Pour des temps tres petits, ces 
corrections modifient la decroissance exponentielle de l'amplitude. La for- 
mule habituelle de passage de la largeur au temps de vie est done modifiee 
a petit temps si les deux conditions suivantes ne sont pas satisfaites : 

(m-6)T»l et r/(m-6)<l. 

Ces conditions sont respectees pour les particules quasi-stables mais vio- 
lees pour les resonances. A grand temps, ces corrections sont aussi domi- 
nantes mais si les particules sont quasi-stables, les deux termes sont alors 
trop faibles pour etre observes dans la propagation. 

On constate que ces deux types de correction sont importants a tout temps 
dans le cas des resonances. Ce probleme ne nous concerne cependant pas ici puisque 
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les resonances ne se propagent pas macroscopiquement. Dans l'appendice, figure la 
derivation du comportement non exponentiel a grand temps, qui est en puissance 
inverse de T. Par contre, bien que Ton ait delimite ici l'extension du regime expo- 
nentiel, le comportement non exponentiel a petit temps n'a pas ete calcule. Certains 
calculs effectues en mecanique quantique montrent que, dans ce domaine temporel, 
la probability dependrait quadratiquement du temps [57|, mais cette prediction n'a 
pas encore pu etre testee. 



4.4 Analyse spatiale de l'amplitude 

Rappelons que l'amplitude de propagation (|4.4|) est donnee par 

A = J d 3 P I(T) e ip L , 

oil L e x D — x P , T = t D — t P , et I(T) est donne asymptotiquement par ( J4.6[ ). Si Ton 
y substitue l'expression asymptotique de I(T), l'amplitude devient 

A{T,L) « J d 3 p?p(z ,p) exp (-i ^z + p 2 T + ip ■ Lj , 

oil ip(zo,p) = Zn(z + p 2 )~ 1//2 ip(zo, p) est l'integrale de recouvrement des paquets 
d'ondes entrants et sortants. Comme ceux-ci ont des energies et impulsions bien 
defmies, cette fonction aura un pic prononce. Si Ton se souvient de la definition du 
pole, z = m 2 — imT, la condition r/m < 1 permet de reecrire l'amplitude sous la 
forme 

A(T,L) w J d 3 p^(m 2 ,p) exp (- —J^=^ t\ exp (-i ^m 2 + p 2 T + i p ■ L) 

P (4-8) 
Cette integrate ne peut pas etre calculee exactement et va done etre evaluee par 
approximation autour du point-selle, puisque la fonction ip(zo, p) a par hypothese 
un maximum prononce. 

II est necessaire de verifier prealablement que si r ^ 0, l'exponentielle decrois- 
sante a une influence negligeable sur la position du maximum dans le domaine tem- 
porel observe en pratique. Dans le cas contraire, la position du maximum dependrait 
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du temps T ! Restreignons-nous a une dimension et prenons pour modele 

L 'integrand est maximal est p = P max , qui est la solution de 

Pmax - P = a P max {P max + m 2 )~ 3/2 oil a = a 2 mTT . 
Si T = 0, P max = P. 

Si T > 0, Pmax = (1 + e) -P avec e = a(P max + m 2 ) -3 / 2 , qui est beaucoup plus 
petit que 1 tant que TT <C m 2 /a 2 , c'est-a-dire pour un grand nombre de temps de 
vie. On exige que la position du maximum ne soit pas modifiee par l'exponentielle 
decroissante a une precision de e. On aura e > e lorsque 

e{Pl ax + m 2 fl 2 



FT > 



o 2 m 



Pour un kaon K$, on a par exemple P = m dans l'experience CPLEAR [^] et 
a = 3 x 10 -2 MeV (incertitude sur la masse). Des lors, le niveau de precision voulu 
sera viole lorsque 

TT > 7 x 10 8 e, 

Par exemple, si l'on ne desire pas une precision superieure a e = sur la loca- 
lisation du maximum, sa position ne changera pas sur 70 vies moyennes. Ce chiffre 
est nettement plus grand que les 20 vies moyennes observees dans CPLEAR ||68|| . 

Notons par P la position du maximum de |^(m 2 ,p)|. L'hypothese de symetrie 
centrale de la phase de ip autour du maximum implique les derivees premieres de ip 
sont nulles au maximum P. Dans le domaine oil elle est non nulle, la fonction ip est 
approximee par une gaussienne : 

V>(m 2 , p) w ^(m 2 , P) exp (-(p -P)W(p- P)) . (4.9) 

W est une matrice complexe symetrique contenant les derivees secondes de ln^ 
evaluees en P. Le produit matriciel est desormais implicite dans les expressions du 
type p W p. On developpe aussi les autres termes de l'integrand autour de P : 



m 2 + p 2 = E + v ■ (p - P) + (p - P) R (p - P) , (4.10) 
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oil E = \Jm 2 + P 2 et v = P/E. La matrice R contenant les derivees secondes 
est reelle et symetrique. On neglige la dependance en p du terme en T puisque sa 
contribution a la localisation du maximum a pu etre negligee. II est done a peu pres 
constant dans le domaine oil ip est maximal. 

Apres l'insertion de ces developpements, l'amplitude de propagation ( |4.8|) se 
reecrit 

A(T,L) = V(m 2 ,P) exp(-mTT/2E) exp(-iET) 

x J d 3 p exp [-zv- (p - P) T + ip • L] exp (-(p - P) (W + iRT) (p - P)) 
= ^(m 2 ,P) exp(-mTT/2E) exp(—iET + i P • L) 
x J d 3 p exp[ip- (L-vT)-p(W + iRT)p\ , 

oil l'on a procede a un changement de variable sur p. L 'integration sur p se fait par 
la formule habituelle de l'integrale gaussienne : 

A(T, L) = 7r 3/2 V(m 2 ,P) exp(-mTT/2E) (det(W + iRT)y 1/2 (4.11) 
x exp(-iET + i P • L) exp f-- (L - vT) (W + iRT)* 1 (L - vT)j . 

Si Ton mesure a la fois la distance L et le temps T de propagation, cette expression 
peut etre analysee comme suit : 

1. La premiere exponentielle contient la decroissance temporelle de l'amplitude 
due a la disintegration possible de la particule. 

2. La deuxieme exponentielle contient les termes d'oscillation spatio-temporelle. 

3. La troisieme exponentielle impose la relation approximative L— vT < a x sinon 
l'amplitude est quasiment nulle ; a x est la largeur de la fonction-poids ip(m 2 , x) 
dans l'espace de configuration, e'est-a-dire qu'elle depend de la precision avec 
laquelle on peut localiser la source et le detecteur. 

On va s'interesser a la valeur de l'amplitude dans la direction de v puisque 
e'est la seule direction dans laquelle elle est non negligeable. Le temps T n'etant pas 
mesure dans les experiences, on calcule la probabilite integree sur un grand intervalle 
de temps dont l'etendue peut etre prise de — oo a +oo sans grande erreur grace a la 
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gaussienne en T. La probability integree sur le temps, dans la direction z (L = Le z ) 
fixee par v = ve z = (P/E) e z , se calcule a partir de la formule 

/+oo 
dT \A(T,Le z )\ 2 . 
-oo 

En y inserant l'expression de l'amplitude, on obtient 
V(Le z ) ~ 7T 3 |^(m 2 ,P)| 2 |det(W^ + i J RL/t;)|- 1 



(4.12) 



/+oo 
exp 
-oo 



mr 
IT 



T — 



W 2 



A\W ZZ + iR zz L/v\ 



(L — vTy 



Les indices zz sont des indices matriciels. La dispersion a ete negligee dans les termes 
en RT en y substituant T = L/v. On neglige aussi les termes en T 2 . Les facteurs de 
normalisation sont omis car ils ne jouent pas de role dans notre analyse. L'integration 
sur T donne 



(4.13) 



On reconnait la formule de la probabilite de desintegration d'une particule en fonc- 
tion de la distance, dans sa forme relativiste. 

Remarquons que grace aux fonctions delta figurant dans la fonction-poids (p(p), 
l'energie E et l'impulsion P sont egales a l'energie totale et l'impulsion totale en- 
trantes (et sortantes), a une incertitude a p pres dependant de la taille des paquets 
d'ondes dans l'espace de configuration. On a done 

p 

1 out 



-y" = — — — 

E E in E, 



out 



OU 



Q-K 



out 



K'-Q' 



E Pl (Q) - E Pf (K) - E out = E Df (K') - E Dl (Q') 



La prochaine etape consiste a appliquer cette formule de l'amplitude et de la 
probabilite de propagation a un melange de particules. Ce travail fait l'objet du 
prochain chapitre. 
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4.5 Appendice : integration sur p dans l'ampli- 
tude 

L'objectif est de calculer la valeur asymptotique T^oo de l'integrale flOl ) : 

J(T) = J dp\(p)G(p 2 ) e- ip ° T . 

Le propagateur est donne par l'equation ( |2.(j| ) ou fl2.10| ). La fonction-poids (p(p) n'est 
differente de zero que sur un intervalle fini. Nous prendrons comme modele 

v (p) = (p 2 - Ml) n (p 2 - M\) n n(p 2 , p) pour < M 2 < p 2 < M 2 , 
fip) = sinon. 

p° est positif ainsi que le nombre n. La fonction Q est analytique en p 2 . Le comporte- 
ment symetrique de ip(p) a ses deux bornes a ete choisi pour simplifier le calcul mais 
n'est pas indispensable. L 'exigence d'analyticite de Q n'est pas trop contraignante 
car on peut toujours approximer une fonction reguliere sur un intervalle fini par un 
polynome (qui est analytique). 

Sous le changement de variable z = p 2 = p° 2 — p 2 , l'integrale I(T) devient 

I(T) = - dz(z+p 2 y^(z,p)G(z)e-*V^ T . 

Le propagateur G(z) a des points de branchement aux seuils de contribution des 
etats a plusieurs particules. Soit z = b 2 le premier point de branchement qui cor- 
respond au seuil de production de plusieurs particules reelles. La coupure sur l'axe 
reel commence done en z = b 2 . On fait l'hypothese que les autres points de branche- 
ment sont superieurs a Mf . La continuation analytique de G(z) sur la deuxieme 
feuille de Riemann est notee comme au chapitre precedent Gu(z) et a un pole en 
z = zo = m 2 — imT . On va differencier les cas selon les positions du pole et du seuil 
b, et suivant que la particule est ou non instable. 

1. Particule instable : b < Mi < m < M 2 et T ^ 
On choisit le contour d'integration defini par 
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b 2 M i 



$ M 



i zo 



*i / 



L 'expression analytique des contours V\ et P 2 est 



Vj: z= (-iu + JMf + p* -p 2 (j = 1,2) 



ou uj varie de zero a us^ sur V\ et de us^ a zero sur 7^2- 
L'expression analytique du contour V est 

V : z — (-iu^ + JM^ + p 2 ) 2 - p 2 . 



ou M varie de Mi a M 2 . On prend la limite uj^ — > oo. 
L'integrale recherchee est egale a 

/(T) = J + J x + J 2 + J' . 

J est la contribution du pole z tandis que J x , J 2 et J' sont les contributions 
respectives de Vi, V%, V . Leurs expressions analytiques s'ecrivent 

(a) Contribution du pole : 

Zix (^o + p 2 ) 2 ip(z ,p) exp ^-i^o + p 2 7^ 



J 



Ztt^o + P 2 ) 2 ^(zb,p) exp I -i\J m? + p 2 T 



T 



2 y/mF^pp 1 

On a defini zo = m2 ~ iTnT et fait l'hypothese que T/m <C 1. Comme on 
l'a deja dit, cette relation est tres bien verifiee dans la majorite des cas. 
Par exemple, T/m w e>(l(T 14 ) pour le Posant AM = |m - Afi >2 |, 
l'ordre de grandeur de J est donne par O(J) ~ m 2n ~ l (AM) 2 " a -1 " 7 ^ 2 . 
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(b) Contribution des contours en faucille : 

j j = i(-iy exp(-i v /M2+p 2 T) J™ du<p(z(uj),p) G h {z{uj)) e~ wT . 

Pour T grand , la contribution dominante a l'integrale vient des valeurs 
de oo proches de zero a cause de l'exponentielle decroissante. Comme 
l'integrand tend vers zero lorsque u tend vers zero (ip(M 2 ) = 0)), le 
comportement asymptotique de l'integrale depend du type de conver- 
gence vers zero de l'integrand. On y substitue y = wT et Ton effectue un 
developpement en 1/T : 

Jj = -{\)< 2 -^2 n n\ (M 2 2 -Ml) n (M] + p 2 ) n/2 n(Mlp) G n (Mj) 
x [iTY^ exp (-i^M 2 +p 2 T) . 

Les corrections a cette formule sont en (AAfT)"' n+2 ', ou AM Mj — 
m, c'est-a-dire de l'ordre de grandeur de l'incertitude sur la masse de 
la particule. On fait implicitement l'hypothese que la fonction Q(z, p) 
diverge moins vite que exp(— uT) sur les contours Vj quand ui tend vers 
l'infini, ce qui recouvre une large classe de fonctions. Par exemple, les 
gaussiennes conviennent, alors que leur integrate diverge sur des contours 
de type demi-cercle a de rayon infini. L'ordre de grandeur de Jj est donne 
par Q(Jj) ~ m 2 ™- 1 (AM) n_1 

(c) Contribution du contour a l'infini : 

J' = e~»°° T Q dM VM f +p2 ¥> (z(M)) G n (z(M)) exp (-iy/AP + p*T) . 

Si fl satisfait aux memes conditions a l'infini que ci-dessus, J' ~exp(— uj^T) 
et tend vers zero lorsque ui^, — > oo. 

En conclusion, la contribution du pole est une exponentielle decroissante en 
T tandis que les contributions dues aux bornes du spectre d'energie sont en 
puissances inverses de T. 

Pour des temps petits, on a O(J) ~ m 2n ~ l (AM) 2 " done 0(Jj/J) ~ (AM T)- {n+ ^ 
Comme les interferences entre Jj et J disparaissent par moyenne sur de petits 



APPENDICE : INTEGRATION SUR P° DANS L 'AMPLITUDE 



81 



intervalles de temps, la contribution de Jj par rapport J a la probabilite sera 
plutot donnee par (Jj/J) 2 ~ (AMT)^ ( " +1) , qui est notable pour AMT < 1. 
En deca de cette borne, le calcul asymptotique de Jj n'est plus valable car on 
a neglige lors du calcul des termes en (AMT) -(n+2 ). 

Dans le cas d'une particule quasi-stable, par exemple le K s , la masse est 
mesuree a une precision de 10~ 2 MeV. On a done AM pa 10~ 2 MeV et les 
termes non exponentiels contribueront notablement pour T < 10~ 19 s, ce qui 
est inobservable dans la propagation puisque le temps de vie du K s est de 
0.89 x 1CT 10 s. Pour le B°, AM pa 2 MeV et les termes non exponentiels 
contribueront notablement pour T < 10~ 22 s, ce qui est inobservable dans la 
propagation puisque le temps de vie du B° est de 1.29 x 10~ 13 s. 
Dans le cas d'une resonance, par exemple le A(1232), AM pa 2 MeV et les 
termes non exponentiels contribueront notablement pour T < 10~ 21 s, qui 
est grand par rapport a l'inverse de la largeur egale a 5 x 1CT 24 s. Quel 
que so it le domaine temporel, la propagation des resonances ne peut jamais 
modelisee par la seule contribution du pole : il faut tenir compte des conditions 
de production et de detection. Cependant, les resonances ne se propagent pas 
macroscopiquement done il n'est pas crucial de pouvoir calculer par exemple 
leur temps de vie en fonction de leur largeur. 

Dans le cas de la particule quasi-stable, la contribution du pole ne domine 
pas non plus a grand temps. Les termes non exponentiels prennent le dessus 
lorsque 

TT > 2(n + l) ln(AMT) = 2(n + 1) (ln(TT) + ln(AM/r)) . 

Pour le K° s , soit n = 1/2 et AM/T w C(10 10 ), de sorte que TT > 82, ce qui 
est inobservable. Pour le B°, soit n = 1/2 et AM/T pa 0(1O 9 ), de sorte que 
TT > 75, ce qui est inobservable. 

En fin de compte, I(T) est donnee en tres bonne approximation pour les 
particules quasi-stables par 




,2T 



(4.14) 
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2. Particule instable : M 1 < b < m < M 2 et T ^ 

Dans ce cas-ci, le seuil de production de plusieurs particules reelles se trouve 
dans la domaine d'energies permises. II va generer des contributions en puis- 
sances inverses du temps, qui sont les corrections que nous avons examinees 
qualitativement au chapitre (equations (|2.15 ) et ( |2.19|) ). 
On choisit le contour defini par 



1 ^ 



3. 



M, 



2> 



Mo 



/ 



/ 



/ 



/ / 
/ / 
/ / 
/ / 
m / / 

Jr'j* <P 



/ 



/ 



Les contributions des contours Vj sont les memes que dans le cas precedent, 
sauf que Ton substitue G(z) a Gn(z) sur V\ puisque ce trajet se trouve sur la 
premiere feuille de Riemann. La valeur asymptotique T— » oo de Jj ne changera 
pas. 

La somme des contributions des trajets P&. serait nulle si ces trajets se trou- 
vaient sur la meme feuille de Riemann mais ce n'est pas le cas. Leur contribu- 
tion totale Jb s'ecrit 

■h = If dz(z+p 2 )-^(z,P) (G H (z) - G(z)) e-'V^T 
2 Jv b 

= _ te -*V^T f°° ^(^p) (GjAz) - G(z)) e-" T . (4.15) 
Jo 

L 'expression analytique du contour Vb est donnee par 

V b : z = (-iu; + y/fe 2 + p 2 ) 2 - p 2 , (4.16) 
ou L) varie de zero a et Ton prend la limite Uqq — > oo. L'integrale I(T) est 
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la somme de toutes les contributions des contours : 

I(T) = J+J 1 + J 2 + J b + J'. 

Le comportement asymptotique de J& est etudie de meme maniere que celui des 
Jj. Pour T grand, la contribution dominante a l'integrale vient des valeurs de 
u> proches de zero a cause de l'exponentielle decroissante. Comme l'integrand 
tend vers zero lorsque uj tend vers zero (Gu(z) — > G(z)), le comportement 
asymptotique de l'integrale depend du type de convergence vers zero de l'inte- 
grand. On substitue y = ujT et on developpe l'integrand en 1/T. Si l'on se 
souvient que Gjj(z) est defini par prolongement analytique de G(x + ie) en 
dessous de la coupure (equation (|2.17|) ) et que sa discontinuite a travers la 
coupure est donnee par l'equation (p.5|), 



G(x + ie) — G(x — ie) = 2ti p{x) , 

ou p(x) est la fonction spectrale, on a que 

G n (z) -G{z) =2irp(z). 

Pour fixer les idees, prenons le systeme Kirn dont nous avons explicitement 
calcule la fonction spectrale p(z) (equation ([2.22 )) . Le seuil a deux particules 
est en z = b 2 , avec b = 2m^. 

1 v(z) 
P\ z ) = ~1 o 



7T [z — m 2 — u(z)) + v 2 (z) 
oil m est la masse du kaon et 



L'evaluation de l'integrale Jb donne 

J b = -i (-2^)1 T-*' 2 1- J b2 \ p2)1/ ' 2 V(b\ P) e-tV^T 

Airb (b 2 - m 2 - u(b 2 )) 

Les corrections a cette formule sont de l'ordre de (QT)~ 5 / 2 ou Q = m — b 
est l'energie disponible lors de la disintegration de la particule. Puisque nous 
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avons calcule la largeur de K -^2tt (equation ( |2.23| )), on peut remplacer g 2 par 



son expression en fonction de T, m et Q. L'ordre de grandeur de Jb est donne 
par 

o{j h )^{Q T r^^-- V {b\ V ). 

Q m 

Pour des temps petits, l'ordre de grandeur de J est donne par 0( J) ~ ip(m 2 , p) 
et Ton obtient 

qui est beaucoup plus petit que 1 si QT ^> 1 et T/Q <C 1. C'est le cas pour les 
particules quasi-stables. Par exemple, T/Q pa 10~ 14 pour le Kg. Par contre, ce 
n'est pas vrai pour les resonances. Par exemple, T/Q pa 0.8 pour le A(1232). 
En dega de la borne QT pa 1, le calcul asymptotique de Jb n'est plus valable, 
puisque Ton a neglige des corrections en (QT)~ 5 / 2 . 

Pour des grands temps, Jb prend le dessus pour 

rT-31n(rT) > 51n(Q/r). 

qui se traduit par TgT > 165 pour le K s , TlT > 202 pour le K\ et Tl,hT > 
157 pour le B° L H (pour ce dernier, on a pris Q = 1 MeV). Apres un aussi 
grand nombre de temps de vie, la probabilite de detection d'une particule est 
inobservable et ces corrections non exponentielles sont done indetectables. 

Pour les resonances, la contribution de Jb est comparable ou superieure a celle 
du pole pour tout domaine temporel. 

Dans le cas des particules quasi-stables, l'integrale I(T) est done donnee en 
tres bonne approximation par la meme expression qu'au premier cas (equation 

(EH) = 



I(T)paZ<K (M 2 + p 2 ) ~ 2 <p{zo,p) expf-^o + pST 



Quant aux resonances, elles ne se propagent pas macroscopiquement, leur 
largeur etant du meme ordre de grandeur que les energies typiques des proces- 
sus. 
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Particule stable : M x < m < M 2 < b et T = 

Ce cas se traite similairement au premier cas sauf que le pole est reel. On peut 
reprendre le resultat (^4.14j) . II n'y a pas de decroissance exponentielle de la 
contribution du pole puisque celui-ci est reel. Les corrections a grand temps 
ne dominent jamais la contribution du pole. 

Particule stable : Mi < m < b < M 2 et T = 

Ce cas se traite similairement au deuxieme cas sauf que le pole est reel. On 
peut reprendre le resultat en posant T = 0. Le rapport Jb/J vaut Jb/J ~ 
(mT)~ 3 ^ 2 g 2 /Q 2 , ou g est la constante de couplage avec les particules produites 
au seuil. La contribution de Jb est inobservable en pratique a petit temps et 
sous-dominante a grand temps. 

m < Mi ou m > M 2 

Le pole se trouve en dehors du contour d'integration et ne contribuera pas. 
La propagation macroscopique de la particule sera inobservable. Les energies- 
impulsions des etats initiaux ou finaux ne permettent pas que la particule 
intermediaire soit sur sa couche de masse et cela n'a pas de sens d'etudier sa 
propagation comme celle d'une particule quasi-reelle. 
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Chapitre 5 

Oscillations en theorie des champs 



5.1 Propagation d'un melange de particules 

Maintenant que la propagation d'une seule particule a ete etudiee en detail, nous 
sommes prets a attaquer le probleme de la propagation d'un melange de particules 
dans le cadre du modele sophistique ou les etats entrants et sortants sont modelises 
par des paquets d'ondes. Remarquons que l'etat intermediaire n'est pas un paquet 
d'ondes. Cette methode n'appartient done srorie dite du traitement 

des oscillations par paquet d'ondes, qui consiste a considerer la particule oscillante 
comme un paquet d'ondes dans le cadre de la mecanique quantique |52|, [73|]. Notre 
methode releve bien de la theorie des champs, ou les ondes planes utilisees dans le 
calcul d'amplitudes de Feynman ne sont en fait que des approximations remplagant 
les paquets d'ondes |65| . 

Apres quelques rappels de notation, nous appliquerons la formule de prop- 
agation calculee au chapitre precedent a une superposition de propagations de 
differents etats propres de masse. Les etapes de l'analyse seront identiques. Suite 
a l'etablissement de la formule de l'amplitude, nous en effectuons une analyse tem- 
porelle, lors de laquelle nous verifierons la non observabilite de nouvelles corrections 
non exponentielles dues au melange. L 'etude de la propagation du melange continue 
par une analyse spatiale de l'amplitude et le calcul de la probability de detection 
integree sur le temps. La formule de probability obtenue sera analysee terme par 
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terme. Nous retrouverons les termes de disintegration et d'oscillation, quoique leur 
forme sera differente de celle obtenue dans le modele simplifie. Par contre, les ter- 
mes de decoherence et de localisation des interactions n'apparaissaient pas dans le 
modele precedent. lis imposeront une serie de conditions d'observabilite des oscil- 
lations. Pour terminer, nous passerons en revue les eclaircissement que ce modele 
apporte aux differentes questions posees au chapitre [!] ainsi que les points com- 
muns et differences de notre approche avec les autres traitements existant dans la 
litterature. 

II est peut-etre utile de rappeler les notations utilisees dans la section |3.4| pour 
decrire un processus de propagation avec melange. 



5.2 Amplitude pour un melange 

Le processus de propagation garde la meme forme que celui du modele simplifie 



de la section 3.4 



Pj(q) ^ P F (k) + v a (p) 



\ 



u p (p) + Djffl ^ D F {k') 



La saveur a de la particule intermediaire v a produite dans la region de (t P ,x P ) est 
identified par l'etat sortant Pp(A;) tandis que la saveur f3 de la particule vp detectee 
dans la region de (t D ,x fl ) est identifiee par l'etat sortant Dp{k'). Si ce n'est pas 
possible (ex : K°,K° — > 7r + 7r _ ), on precede a une sommation sur les differentes 
saveurs. 

Les notations ne changent pas : le propagateur complet Gp a (x' — x) symbolise 
la propagation de particules de saveur a produites en x en particules de saveur (3 
detectees en x' . Ce propagateur est diagonalise (equations ( |3.3|) et (|3.6|) ) en tres 
bonne approximation par des matrices V constantesQ : 

G Pa ( P 2 ) = (V- 1 G D (p 2 ) V)j3 a = E V p f Go^p 2 ) V ja , 



1 Rappelons que la matrice V telle qu'elle est definie ici est la transposed de la matrice V 
diagonalisant les kets dans la derivation en mecanique quantique presentee au chapitre |l|. 
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oil Gn^jij) 2 ) = i (p 2 — Zj) 1 . Le nombre complexe Zj est le pole du propagateur 
complet associe a l'etat j, de masse rrij et de largeur Tj, auxquelles il est relie par 

Zj = m) - mtjYj . 

La propagation macroscopique implique qu'il s'agit de particules quasi-stables, c'est- 
a-dire que rrij^Tj. 

L'amplitude du processus global peut s'exprimer comme une superposition 
lineaire d'amplitudes Aj de propagation d'etats propres de masse : 

A(a^/3, T,L) = J2 V p f A, V ja , (5.1) 

j 

oil l'amplitude partielle Aj s'ecrit comme l'amplitude de propagation d'une particule 
isolee (equation ( f4.4|) ) : 

Aj = J d A V <p(p) G DjJ (p 2 ) e -*(*x>-*iO . (5.2) 

La fonction-poids ip(p) symbolise comme auparavant (equation(|4.5|)) l'integrale de 
recouvrement des paquets d'ondes : 

<p(p) = j [dq\ <f> Pl (q, Q) J [dk] <j>* PF (k, K) J [dqfl <j> Dl (q', Q') j [dk'\ ^ (k', K') 

x (2tt) 4 5 (4) (£; + k' -q-q')5 i4 \p-q + k) M P {q, k) M D {q' , k') . (5.3) 

En ce qui concerne la signification des indices dans ce chapitre, les indices inferieurs 
i et j seront des indices de saveur tandis que les indices superieurs seront des indices 
vectoriels ou matriciels (dans l'espace des impulsions ou dans l'espace de configura- 
tion). Notons aussi que le produit matriciel sera implicite dans les expressions du 
genre LVTL. 

5.3 Analyse temporelle de l'amplitude de melange 

Suivant la procedure etablie pour etudier la propagation d'une particule isolee, 
nous effectuerons l'analyse temporelle et spatiale de l'amplitude avant de calculer la 
probability correspondante. 
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II est plus facile pour commencer de ne considerer que l'amplitude partielle Aj 
(equation (|5.2|) ). L'analyse temporelle consiste a integrer sur p°. On recourt a la 
meme methode d'integration de contour qu'au chapitre precedent (voir l'equation 
f)) et Ton obtient 



Aj= I d*pij){m),v) exp I ™ jTj T ) exp(-iJm 2 + p 2 T + ip-L) , (5.4) 

•> ' ' 2Jm 2 + p 2 I 



ou tp(m 2 ,p) = Zn( Zj + p 2 ) l / 2 <f(zj,p). 

A ce stade, on pourrait se demander ce que sont devenues les corrections non 
exponentielles a l'amplitude de propagation, qui ont ete etudiees en detail au chapitre 
precedent dans le cas de la propagation d'une particule isolee. En premier lieu, nous 
avons ete confronted a des corrections dues aux conditions de production et de 
detection, dependant de la forme des paquets d'ondes initiaux et finaux. Dans le cas 
du melange de particules etudie dans ce chapitre, ces corrections peuvent se calculer 
pour chaque amplitude partielle selon la methode expliquee au chapitre precedent 
et restent negligeables. En second lieu, nous avons rencontre des corrections dues 
aux seuils de production des etats a plusieurs particules. Rappelons qu'elles se mani- 
festent par des discontinuites (seuils) de la derivee de l'energie propre en fonction 
de p 2 . Or la methode de diagonalisation par des matrices constantes approxime 



l'energie propre par une constante (equation ( |3.5|) ). Par consequent ces seuils et ces 
corrections disparaissent dans cette approximation. 

II serait neanmoins bon de tester 1' influence des corrections dues aux seuils 
de production sur la propagation de particules melangees, pour verifier qu'elles ne 
perturbent pas les oscillations. Revenons done un pas en arriere, a la diagonalisa- 
tion exacte par des matrices dependant de z = p 2 . En particulier, considerons la 
diagonalisation exacte du propagateur complet du K°K° (equation ( |3.4j) ). 



-iG{z) = W{z) \W (z). 

V (z-m 2 - f 2 (z)) J 

Nous allons negliger la violation CP et approximer les etats de propagation K$ et Kl 
par les etats propres sous CP, K\ et K 2 . L 'incidence sur le melange des corrections 
dues a un seuil en z = b 2 peut etre etudiee en effect uant, par exemple, l'analyse 
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temporelle de l'element non diagonal G q(p 2 ). En effet, si nous definissons un Ip a (T) 
similaire au I(T) apparaissant a la section |4.3| par 

A(a^/3,T,L) = Jd 3 pI Pa (T) e*' L , 

on peut ecrire 



/oo(T)= dp\(p)G - (p 



2^ e -i P °T 



oil ip(p) est la fonction-poids definie en ( |5.3| ) et Ton a extrait de l'expression ma- 
tricielle ci-dessus l'element 00 du propagateur : 

r ^oiCp^foV) | iWm^Wg^) ,,,, 

p 2 - m( - hip 2 ) p 2 -mi- f 2 {p 2 ) 



On reprend la methode du chapitre precedent, section (4J| L 'integration s'effectue 
sur le meme contour dans le plan complexe avec la difference que le contour contient 
maintenant deux poles quasiment degeneres. L 'integrate Iqo(T) se decompose en la 
contribution J des poles zi )2 , la contribution J 12 des seuils de <p(p), et la contribution 
Jb du seuil de production de plusieurs particules (on suppose qu'il n'y en a qu'un 
dans le domaine d'energie permis par les conditions experiment ales) : 

I m (T) = J+J 1 + J 2 + J b . 

J est donne cette fois par les residus de deux poles : 



J = W i(z 1 )W-q 1 (z 1 )tt{z 1 + p 2 ) 3 (p( Zu p) exp ( -iyj Z! + p 2 T 



+ W 02 (z 2 )W- (j 1 (z2)7i(z 2 + p 2 )-^(z 2 ,p) exp (-i^z 2 + p 2 T) . 



Comme la violation CP est negligee, les matrices de diagonalisation evaluees aux 
poles sont donnees par la formule ( ft.8|) de passage a la base propre de CP : 

w( Zl ) - w(z 2 ) - -L 

Ces expressions ne sont pas valables loin des poles. On en tiendra compte dans 
revaluation de Pour evaluer l'ordre de grandeur de J, on developpe l'argument 
des exponentielles en utilisant la definition des poles, Zj = m 2 — irrijTi : 

m 2 - m 2 mTjT 



Zj + p 2 T = —iET - i 



2E 2E 
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ou m est la masse degenee et E = y/rri 2 + p 2 . L'ordre de grandeur de J est alors 
donne par 

O(J) ~ - <^(m 2 , p) e"^ 2 e - iAmT " ArT / 2 - 1 , (5.6) 
m 

ou Am = mi — m-z, AT = Tx — T2 et l'on a approxime E = m. 

Les contributions Ji et J 2 des seuils de (pip) peuvent se calculer separement 
sur les deux termes de la somme ( |5.5|) et ce calcul n'apporte rien de nouveau par 
rapport au chapitre precedent : ces corrections ne dominent l'exponentielle que pour 
des temps tres grands, quand l'amplitude totale est devenue indetectable. 

Examinons la contribution de J^. On a vu que depend de la difference entre 
le propagateur complet et son prolongement analytique (equation ( |4.15p ) : 



Jb = - ie -iV^^ J^du(piziuj),p) (G - 0JI (z)-G - (z)) e 



L'expression Gu{z) — G(z) peut etre calculee dans la deuxieme representation spec- 
trale (equation Q2.6Q ) : 

G o,ii(z) - G o{z) = -i (lIoo,j/(2) - n o(^)) Gqq{z) Gqo , 

ou n o(^) est l'element 00 de la matrice d'energie propre. Lors de revaluation asymp- 
totique de Jb (T grand), seul le domaine proche de u = (c'est-a-dire z = b 2 ) va 
contribuer notablement, en raison de l'exponentielle decroissante. II suffit done de 
connaitre l'expression de l'integrand pres du seuil. Juste en dessous de l'axe reel 
(z = x — ie), la representation spectrale de l'energie propre (equation ( |2.7| )) im- 
plique que 

n o , u(z) ~ n o [z) = n oS (x + ie) - n og ix-ie) 
= 2iXm Il o(a; + ie) . 

II reste a calculer TmU o(x + ie). Comme la renormalisation de la masse ne concerne 
que la partie reelle de l'energie propre, on a l'egalite 

Xm n o (x + ie) = Xm / g (x + ie) . 

Lors de la diagonalisation du propagateur du K°K° a la section [3.3| , la relation entre 
l'element non diagonal de l'energie propre evalue au pole et les quantites mesurees 
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experiment alement a ete donnee dans l'equation : 

2m/ o(x + ze) = —Ima 

= 2 ( m L T L ~ m s T s ) 

- I m (r L -r 5 ) 

^ --mAT, 
2 

puisque suivant les notations de l'equation (|3.11|) , b — 0, e = 0, K\ — K s et 
K2 = Kl en l'absence de violation CP (ce 6-ci n'est pas la valeur du seuil mais un 
des parametres du propagateur). Comme les canaux de disintegration principaux des 
K° et K° sont identiques (disintegrations en deux pions), on modelise Im foo(x + ie) 
par la meme dependance fonctionnelle que celle trouvee pour Im foo(x+ie) (equation 
(|2.21|) ), de sorte que la valeur au pole soit celle indiquee ci-dessus et que la valeur 
au seuil soit nulle : 



I x — b 2 

Im foo (x + ie) = — m AT 



2 V m 2 ~b 2 

Des lors, on peut ecrire pres du seuil 



z-b 2 

Gm,u{z) ~Gw(z) = -mAT J— —G q(z) G Q q n (z) . 

V m z — tr 

Substituons y = uT dans l'equation de Jb et effectuons le developpement en 1/T 
de l'integrand, en utilisant la parametrisation de z en fonction de u donnee par 
l'equation (|4.16|) : 



z = -u 2 + b 2 - 2iu \/b 2 + p 2 

= b 2 - 2*| ^ + p2 + o (^) . 

On obtient 

J b - j T-^ e-Vwr v {b\ p) m AT ^ + ^ 

V m 2 — b 2 
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ou l'on a approxime dans la deuxieme ligne b 2 + p 2 = m 2 et l'on a neglige l'energie 
propre dans les propagateurs evalues en b 2 car elle est negligeable loin du pole. On 
en conclut que l'ordre de grandeur de Jb est donne par 

AF 1 

(D(J b ) ~ (QT)-V 2 — -^(& 2 ,P), (5.7) 
Q m 

ou b est l'energie du seuil et Q = m — b est l'energie cinetique des produits de 
disintegration. 

On peut maintenant comparer les ordres de grandeur de J (equation (|5.6|) ) et 
de Jb (equation (|5.7|)). 
Pour T petit 
0(J b ) ~ 0(J) si 

(QT y 3/2 — ~ s in ~ J !— , 

ou encore si T ~ Q~ l , puisque l'experience dans le cas du K°K° |H| et la theorie 
dans le cas du B°W @ donnent O(AT) < O(Am). Comme Q varie de 220 MeV 
pour les kaons a 1 GeV pour les mesons B, les corrections non exponentielles ne 
dominent a petit temps que pour T < 10~ 24 s. La propagation spatio-temporelle est 
inobservable a des temps aussi petits. Le calcul asymptotique de Jb n'est plus valable 
non plus en dega de cette borne puisque l'integrale Jb a ete developpee en 1/QT. Les 
corrections non exponentielles dues aux seuils de production sont done inobservables 
a petit temps dans la propagation spatio-temporelle des particules quasi- stables. 
Pour T grand 

II faut traiter separement les cas des kaons et des mesons B. 
1. Kaons : Y s ^> T L 



0{J) ~ 0(J b ) si 



ou encore si 



e~ r ^ 2 ~ (QT)- 3 / 2 ^ 

p p3/2 
~ (Tr T)' 3/2 S L 

g 5 /2 ' 

l' L 7 - 3 ln(r L T) ~ 2 In f ~ 174 . 



e'est-a-dire si T^T ~ 190, ce qui est inobservable. 
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2. Mesons B : T L « T H 
O(J) ~ 0(J 6 ) si 



r ff T/2 



(QT)- 3 / 2 



Ar 



e 







Arrt 



ou Ton a approxime le sinus de l'oscillation par sa valeur moyenne. Cette 
condition se reecrit 



ou a q = AT/T H est estime a = 4 x 10 3 pour le et a a s = 10 1 pour le 
£° |70|. On calcule que 



Ces temps sont beaucoup trop grands pour que Ton puisse esperer une cor- 
rection non exponentielle mesurable : la probability de detection aura decru a 
une valeur indetectable. 

Nos estimations des corrections non exponentielles en theorie des champs sont iden- 
tiques aux formules theoriques obtenues par Chiu et Sudarshan 0] (formule (4.50) 
de cet article ; leurs evaluations numeriques sont par contre assez fantaisistes) et 
par Wang et Sanda |72| (formule (59) de cet article, mais leur formule (61) d'or- 
dre de grandeur des corrections est incomprehensible et incorrecte du point de vue 
des unites). Ces auteurs calculent ces corrections dans le cadre de la mecanique 
quantique et proposent des extensions du formalisme de Wigner-Weisskopf. 

En conclusion, les corrections non exponentielles ne sont pas plus visibles dans 
la propagation macroscopique d'un melange de particules que dans la propaga- 
tion d'une seule particule. Les seules contributions importantes pour la propagation 
macroscopique viennent des poles des propagateurs. Par consequent, la formule de 
l'amplitude partielle Aj (|5.4j ) peut servir de base rigoureuse pour l'analyse spatiale 
de l'amplitude. 




B° d : T-3 ln(r H T) ~ 153 T H T = 168 , 

B° s : T-3 ln(r H T) ~ 147 T H T = 162 . 
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5.4 Analyse spatiale de l'amplitude de melange 



L'integration sur p de l'amplitude partielle Q5.4Q demande une certaine pru- 
dence. En effet, nous sommes interesses par les termes d'interference AiA* (i ^ 
j) qui apparaissent dans la probability car ces termes contiennent les oscillations 
dependant de l'espace-temps et de la difference de masse Arriij = — rrij\. II s'agit 
d'evaluer ces termes d'oscillation a une precision superieure a Amy /P. Rappelons 
que P est la position du maximum de ipim'j, p). Cependant, les oscillations ne sont 
observables sur des distances macroscopiques que si Arriij /P est extremement petit. 
Par exemple, Arriij /P ~ O(10~ 14 ) pour les kaons neutres, Amy jP m O(10~ 8 ) pour 
les neutrinos dans l'experience LSND et des valeurs encore plus petites apparaissent 
dans les oscillations des neutrinos atmospheriques et solaires. 

L 'evaluation directe de Aj par une approximation autour du point-selle ne 
fournira en principe pas la precision voulue dans les termes d'oscillation pour deux 
raisons : 

1. Si les positions des maxima de ip(rn^, p) et ^(m^, p) sont notees respect ivement 
Pi et Pj, la difference |Pj — Pj| sera de l'ordre de grandeur de Arriij. Or 
la largeur du pic de la fonction ip est beaucoup plus grande que Arriij. Peu 
importe que Ton utilise Pj ou Pj comme point-selle : la methode d'integration 
n'atteint pas une precision superieure a Arriij/ P. 

2. Pour les particules instables, le terme en T modifie la position du maximum 



de AP = eP, ou e ne peut etre inferieur a 10 (voir section [OP , ce qui est 
beaucoup plus grand que Arriij/ P ~ O(10~ 14 ) pour les kaons neutres. 

Malgre ces obstacles, ce n'est pas la fin des haricots. Nous allons voir que les termes 
d'oscillation sont independants au premier ordre en Amy de la position des points- 
selle Pj et Pj, de sorte que 1' approximation autour du point-selle peut quand meme 
etre utilisee. 

En suivant les memes etapes que pour la particule isolee (voir equation ( [4.1 1|) ), 
on integre par approximation autour du point-selle et l'on obtient 

Aj = ^^(mlPj) exp(-rrijTjT/2Ej) (det(Wj + iRjT)y 1/2 (5.8) 
x exp{-iEjT + iPj- L) exp (-- (L - VjT) (Wj + iRjT)- 1 (L - VjT)) , 
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oil 

Ej = ^m) + P] et vj = ^. 

La probabilite integree sur le temps se calcule a partir de la norme au carre de 
l'amplitude totale ( |5.1| ) et s'ecrit 

V(a - /?, L) ~ Yl V i« V §i V Ja V?- 1 * J dTAi(T, L) A*(T, L) . 

En raison des gaussiennes d'argument (L — VjjT) 2 , la probabilite est tres faible si 
les deux conditions suivantes ne sont pas satisfaites : 

L — VjT = et L — VjT = . 

II est possible de satisfaire ces deux conditions a l'ordre e = Arriij/\Pj\ puisque les 
positions des maxima P» et Pj sont identiques a e pres. On a aussi 

V ; = Vj- + O(e) et Ei = Ej + 0(e) . 

Soit P la position du maximum de i/j(m 2 ,p), pour une masse de reference m qui 
est egale a rrij a e pres. Appelons z la direction de P. Pour les kaons, on peut, par 
exemple, prendre m = (jrii +m,j)/2 et m = pour les neutrinos. La valeur exacte 
de m n'a pas d'importance. Nous allons etudier la probabilite dans cette direction 
z, c'est-a-dire que L sera approximativement parallele aux Vjj . 
On parametrise la deviation des Pjj par rapport a P par 

oil ef y,z = 0(e). Les diverses quantites apparaissant dans l'amplitude peuvent etre 
evaluees avec cette parametrisation et s'ecrivent 

P,-L = (1 + €*■) PL + 0(e 2 ) 

E 3 = ^m) + P] = E + ^ + e)vP + 0(e 2 ) 
v, = W ((l + ^)e 2 + eJe ;c + eJe 2/ )+0(e 2 ) 

oil 

Am 2 = m 2 — m 2 

E = Vm 2 + P 2 
Am? 
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On va maintenant evaluer l'integrale sur T dans la probability. La dependance en T 
de la dispersion est negligee en remplagant T par L/v dans les termes en RT. 

Pour etre tout a fait explicite, nous donnons d'abord l'expression de la proba- 
bility integree sur T pour une direction quelconque L, sans utiliser la parametrisation 
definie ci-dessus : 



x tt^K^P^^K^P,) (detX, + detX;) _1/2 2^ 
x exp(-^ ( V ;XL + v,X*L)) 
x exp(i (p.-p.)-L-i ^~^ (viXiL + Vj^L)) 
x exp(-^L(x + X*)L + ^ (v.^L + v.^l) 



y exp ( -I ( /•;, - /•;, ) 2 + 2i(E t - E 5 ) y ) (5.9) 



ou 



Xj = (^ + ^.|L|/|v J |)- 1 (5.10) 
Y = v, .V, v, • v ; .V : v ; (5.11) 

r« - ^ + ^ (5-12) 

Les termes en T 2 sont bien entendu negliges. 

Nous donnons ensuite la probability dans la direction L = L e z ou la probability 
est maximale, en utilisant la parametrisation definie plus haut des deviations par 
rapport a cette direction : 

V(a^(3, L e 2 ) ~ £ V ia V p f V* a V£ * (5.13) 

x Tr^m 2 ^) ^*(m-,P,) (detX + detX*) -172 2^ (Xf + Xf*)~ 1/2 

/ TijLX (. T ( z z Ei-Ej 
x exp I —J exp U L [p* - p* - 

/ L 2 XrX?**(v?-vff \ ( -{E i -E j f + 2iT {E t -E 3 ) 
x exp — vxx ,' , exp ' 



4 Xf {vff + X?**(v*) 2 J ^ ^ (xf + Xf *) 



t> 2 
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On a effectue un developpement en l'ordre C(e 2 ) dans les facteurs dependant de 
la distance, en ne gardant que les termes en eL ou e 2 L 2 qui peuvent etre non 
negligeables a grande distance. 



5.5 Analyse de la probability 



Dans cette section les differents facteurs apparaissant dans la probabilite d'in- 
terference ( |5. 13|) sont analyses successivement. Posons 



V{a^i3,Le M )~'£v ta V£ 1 v; a V t 



-i * 



3j 



i -J 



x 7T 3 ^(m- , PA ^*(m 2 , P,-) (det X t + det X*) V ' 2^ (x? z + Xf* 



-1/2 



x exp A exp B exp C exp D 



(5.14) 



Les exponentielles e A , e B , e c et e D correspondent dans l'ordre aux exponentielles 
de l'equation ( |5.13|) . 



5.5.1 Disintegration 



Si Tj ^ et/ou Tj ^ 0, alors ^ (equation ( |5.12| )). Dans ce cas, la premiere 
exponentielle e A exprime la decroissance de la probabilite de detection de la particule 
en raison de sa disintegration possible. Notre resultat est relativiste et s'ecrit en 



fonction de T^j comme 




(5.15) 



Comme il s'exprime directement en fonction de la distance, il n'y a pas d'ambiguite 
de passage du temps a la distance. On peut le comparer au traitement non relativiste 
en mecanique quantique qui donnait exp (— (r, + Tj)T/2) dans le repere au repos de 
la particule. L'extension relativiste de cette derniere formule ainsi que sa transfor- 
mation en une expression dependant de la distance peut se faire suivant plusieurs 
prescriptions aboutissant a des resultats differents. 
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5.5.2 Oscillation 



La deuxieme exponentielle e B de l'equation ( p.!4j) fait osciller la probability de 
detection de la particule en fonction de la distance. Au premier ordre en e, l'argument 
de l'exponentielle vaut 



B = iL [pf-p 



Ei — Ej 



iL P (d 



Ami L 

2P 
L 

T OSC 



V 2Ev 



= -2m 



(5.16) 



ou Am?- = mj — m| et Ton a defini la longueur d 'oscillation L° sc pour le melange 
ij paiQ 



(5.17) 



On voit que le facteur d'oscillation ne depend pas au premier ordre en e du choix de 
Pi ou Pj comme point-selle pour evaluer l'integrale ! La raison en est que les etats 




propres de masse sont sur leur couche de masse (comme on l'a vu a la section [O 
lors de l'analyse temporelle de l'amplitude, equation ( |4.7| ) ) puisqu'a grande distance 
seul le pole contribue a la propagation. L'exponentielle e B se reecrit done comme 



(5.1{ 



5.5.3 Decoherence 

La troisieme exponentielle e c de l'equation (|5.14|) montre que l'interference 




disparait a grande distance. Ce phenomene est appele decoherence [74, 53 



Calcul de la longueur de coherence 

Montrons d'abord comment apparait cette troisieme exponentielle. Partons de 
la formule de l'amplitude dans une direction quelconque (equation (|5.9| )) et con- 
2 Si Am?- est negatif, on change le signe de la definition de L°^ c . 
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siderons sa troisieme exponent ielle. Si l'on y introduit la parametrisation autour de 
la direction d'amplitude maximale, l'argument de cette exponentielle devient 

c = -\l(x 1 + x;)l + -^ (v,x 4 l + Vj x;l) 2 (5.19) 

L 2 X zz Xf*{vf-v z ) 2 



4 X? z (v?) 2 + X? z *(v?) 2 

L i ((vf) 2 Xr + {vlfXf + (rffXf- + (rf)-.Vf ') (x zz + Xf * 
~T X? z (vf) 2 + X zz *{vf) 2 

L 2 « - Vj) ((vfX? z + vfxf z )x zz * - (v*X? z * + v)Xf*)X zz ) 
~T X zz (vf) 2 + X zz *(v z ) 2 ' 

ou l'on a conserve les termes en e d'ordre 0(eL) et 0(e 2 L 2 ) car ils peuvent etre non 
negligeables a grande distance. 
Comme cas de figure, posons 



ip(m? Pj) ~ exp 



fr-Pj) 2 > 



Si l'on va rechercher les definitions de Wj, equation fl4.9|) , et de Rj, equation (|4.10|) , 
on voit que Wj est lie a l'incertitude sur la localisation des interactions tandis que 
Rj represente la dispersion en energie : 



(j ab 

W ab = = 5 ab a 2 

3 Aa 2 

xab papb 

Rf = w J ~i E f 6aZ5bz + oie2) 



ce qui donne en utilisant la definition (|5.10|) de Xj : 



En examinant les differents termes de l'argument de l'exponentielle (equation 
(|5.19| )), on remarque que le premier terme (c'est-a-dire la premiere ligne) de la somme 
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domine les deux autres, puisqu'on a toujours la relation L/2Ev ^> a 2 si la distance 
L est macroscopique. L'argument (|5.19|) de l'exponentielle e c se reecrit done 

,z\2 



T 2 



[VJ — V,: 



4 X zz *-\v?) 2 + X t zz -\v?) 2 

Examinons separement le numerateur et le denominateur de cette fraction. 
D'une part, 



v z -v] = v 



Am 



Am?- 



2E 2 



* + (1 - v 2 ) (ef - e z ) = u«=^* 



2E 2 



ou k est un nombre d'ordre 0(1). 
D'autre part, 



X; z *-\v*) 2 + X?'-\vff = 2v z a z x ( 1 + 



.2 2 



2va 2 E* i + 1 J 



ou 77 est un nombre d'ordre 0(1). L'argument de l'exponentielle e° vaut alors 



T 2 



(vf - vj) 2 



4 X zz *- X (vf) 2 



+ X zz -\v z ) 2 



L * 

T coh 



1 + IT] 



L " 

J^coh 



ou la longueur de coherence L1° h pour le melange ij est definie paif] 



AV2 E 2 



■ coh 



0~ r 



(5.20) 



et ou 



2V277 
KEva r 



L'exponentielle e c se reecrit 



(5.21) 




(5.22) 



3 Comme pour la longueur d'oscillation, on change le signe de la definition de L\° h si Am-j est 
negatif. 
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Interpretation de la longueur de decoherence 

Comment interpreter cette longueur de coherence? Montrons d'abord qu'elle 
est beaucoup plus grande que la longueur d'oscillation : 

Llf 1 E 

L°- c \f2i\Kv o v ' 

Meme dans le pire des cas, qui est celui d'une particule relativiste pour laquelle 
y/2-KKV w 0(1), la Ion gueur d'oscillation reste beaucoup plus petite que la longueur 
de coherence : 

L1° h E 
-2- » — > 1. 

TOSC 

u ij a p 

L' interference commence a disparaitre apres un nombre d'oscillations (L^ h / 'L°!j c ) 
egal a E/a p . On peut donner une justification intuitive de ce phenomene par l'ima- 
ge suivante JF5). L'etat intermediaire est represents par la superposition de deux 



paquets d'ondes correspondant a deux etats propres de masse. Si leur vitesse est 
relativiste, la dispersion dans la direction du mouvement est negligeable a cause de 
la contraction de Lorentz. En raison de leurs masses differentes, leurs vitesses de 
group e different de 

Am., 

Av = v\ - v) = vk^- . 

Apres un temps T et une distance L w vT, les paquets d'ondes se sont deplaces l'un 
par rapport a l'autre de 

Am?- 

AL = AvT w k rfL. 

2E 2 

Dans cette image, les oscillations sont les battements des deux paquets d'ondes lors 
de leur deplacement relatif de AL = A, oil A = 2n/E est la longueur d'onde associee 
a l'etat oscillant. La longueur d'oscillation s'obtient en imposant que le deplacement 
relatif des paquets d'ondes soit d'une longueur d'onde : 

AttE 

Pour des etats relativistes, E = P et Ton retrouve la longueur d'oscillation derivee 
ci-dessus. 

Si la taille du paquet est donnee par a x , le nombre maximal d'oscillations est le 
nombre de longueurs d'onde dans le paquet, N max = <j x /\ = E/<Ana p . Si le nombre 
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d'oscillations est superieur a N max , les paquets ne se recouvrent plus et les oscillations 
disparaissent. La longueur de coherence est done donnee par 

E 



T CX>h AT T OSC 

1^*4 — JVrnM l-i 



10 4ira p lJ ' 

On retrouve la formule de la longueur de coherence derivee plus haut a une constante 
pres. 

Par contre, pour des vitesses non relativistes, la dispersion dans la direction du 
mouvement est importante et la decoherence est beaucoup plus lente. 

Revenons a notre facteur de decoherence derive en theorie des champs et exa- 
minons les differentes situations possibles : 

1. Pour les particules instables, la decoherence n'a souvent pas de signification 
car les particules se desintegrent sur une distance beaucoup plus courte. Par 
exemple, pour le systeme Kg — K®, 



coh t/2 r E E 



10 4 



kv m cr„ m, — mi 



'p H"L 

ou l'on a defini L des = v/Y. 

2. Pour les particules stables relativistes, v ^> cr p /E done r] 1 et la decoherence 
est en exp (-(L/ L coh f). 

3. Pour les particules non relativistes telles que v <C cr p /E, le coefficient rj ~ 1/v 
est tres grand et il n'y a pas de decoherence. Le facteur de decoherence de- 
vient une oscillation spatiale de periode beaucoup plus grande que la longueur 
d'oscillation. Nous avons explique dans l'image des paquets d'ondes que cette 
absence de decoherence est due a la dispersion des paquets d'ondes. 

L'approche de la theorie des champs ne permet, cependant, que de fixer une borne 
inferieure pour la longueur de coherence. En effet, elle depend de la position des 
maxima P« et Pj. Par exemple, on peut avoir « = done L coh = oo si |Pi|wj = 
|Pj|mj ! D'autres longueurs de coherence existent mais elles sont superieures a celle 
examinee ici 



51, 75, 62 



Notons que la decoherence fait disparaitre l'oscillation a une certaine distance 
en supprimant le terme d'interference, mais qu'elle ne supprime pas les transitions 
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d'un etat de saveur vers un autre. La probabilite de detection est modifiee par une 
constante au dela de la longueur de coherence et cette modification est observable. 
Par exemple, pour un melange a deux saveurs, d'angle de melange 9 (la matrice de 
melange V est ici tout simplement la matrice de rotation d'angle 8), la probabilite 
de survie de Pet at a devient 

V{a-*a,L) = l ~ \ sin2 20 ■ 

Exemples de decoherence 

Donnons quatre exemples ou la decoherence pourrait en principe jouer un role. 
Le facteur crucial est le rapport entre la longueur de decoherence et la longueur 
d'oscillation, F = L coh /L osc , approximativement egal a E/a p pour des particules 
relativistes. 



1. Dans Pexperience LSND Pb| , Penergie des neutrinos tourne autour de 30 MeV, 
tandis que Pincertitude cr p pourrait etre estimee a 5 MeV Le facteur F 
vaut dans ce cas seulement quelques longueurs d'oscillation et pourrait done 
jouer un role dans les experiences futures. 

2. Pour la detection des neutrinos solaires, Penergie des neutrinos tourne autour 
de 1 MeV et Pincertitude a p peut etre estimee a 0.002 MeV |75| . Le facteur 
F vaut ici 500. On ne peut bien sur faire varier a sa guise la distance Terre- 
Soleil done Poscillation n'est observable que sur la variation saisonniere de 3% 



de la distance [21]. Le modele d'oscillations dans le vide fournit une longueur 



d'oscillation de 3 x 10 7 km done la decoherence n'affecte pas les oscillations 



observables dans ce modele. Par contre, dans le modele MSW [£7]] la difference 
de masse entre les neutrinos est telle que la longueur d'oscillation dans le vide 
tourne autour de 250 km. Le nombre d'oscillations sur la distance Terre-Soleil 
est de 60 000, bien superieur au facteur F. Cette decoherence n'est neanmoins 
pas verifiable, car Pignorance du lieu de production exact et la moyenne sur 
Penergie du neutrino impliquent une moyenne de la probabilite sur un grand 
nombre de longueur d'oscillations. L 'interference est eliminee non seulement 
a cause de la decoherence mais aussi en raison de cette moyenne. Seule la 
decroissance de la probabilite est visible. Bien que cela ne constitue pas un 
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test de la presence de la decoherence, l'observation ou non d'une variation de 
la probability de detection correspondant a la variation annuelle de la distance 
Terre-Soleil eliminerait soit le modele MSW soit le modele des oscillations dans 
le vide. Notons cependant que la moyenne sur les energies des neutrinos peut 
aussi rendre invisible l'influence de la variation de la distance Terre-Soleil dans 
le cas du modele d'oscillations dans le vide J7f], [T7], [78], |75fl . 

3. Si Ton neglige les contraintes des modeles cosmologiques sur la masse des 



neutrinos, la masse du neutrino tau n'est bornee que par 18.2 MeV | 69fl . Si 
la masse du troisieme etat propre de masse est de l'ordre de plusieurs MeV, 
tandis que celles des deux autres etats sont de l'ordre de 1'eV, les longueurs 
de coherence L^ h pour les melanges du troisieme etat peuvent etre inferieures 
a la longueur d'oscillation L°| c des deux autres etats dans les experiences ou 
l'energie des neutrinos est de quelques MeV. Un tel melange de neutrinos 



relativistes et non relativistes a ete etudie par Ahluwalia et Goldman ||79|| , qui 
identifient le troisieme neutrino avec la particule de 33.9 MeV proposee pour 
expliquer une anomalie dans l'experience KARMEN. 

4. Pour terminer, l'emission de neutrinos par des supernovae offre une possibility 
d'etudier soit des differences de masse encore plus faibles (Am 2 ~ 10~ 20 eV 2 ), 
soit d'observer l'effet de la decoherence par l'observation de vagues de neutrinos 
correspondant aux differents etats propres de masse [BO]. 



5.5.4 Localisation des interactions 



La quatrieme exponentielle e D de la formule de la probability (|5.14j) exprime 



une condition sur a x ou une longueur de coherence suivant les cas. En effet, 

Am 2 - . , Am?- 

ou k est un nombre d'ordre 0{1). Supposons que nous ayons affaire a des particules 
telles que le coefficient r\ defini dans le calcul de la longueur de coherence (equation 
( |5.21|) ) est beaucoup plus petit que 1, c'est-a-dire que v^>a v jE. On peut montrer 
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dans ce cas que l'argument de la quatrieme exponentielle devientQ 



II est malheureusement difficile de dire quel terme de la somme domine sans connaitre 
l'ordre de grandeur des masses. 

Si le terme en L est negligeable, on obtient 

D^-2^(^ . (5.23) 

Pour que l'amplitude soit non nulle, il faut done que a x -C L°j c , e'est-a-dire que 
l'incertitude sur la localisation de la source et du detecteur doit etre inferieure a la 
longueur d'oscillation pour que cette derniere soit observable. 
Si le terme en L domine, on obtient 

,2\2 



■13 / 

II s'agit d'un nouveau terme de decoherence qui peut accelerer la decoherence vue 
a la section precedente si les particules sont non relativistes (tout en respectant la 
condition )jCl). 

Pour des particules non relativistes telles que tj^>1, une analyse similaire peut 
etre menee mais n'est pas concluante en raison du grand nombre de parametres 
inconnus. 

Le facteur D depend aussi du choix des points-selle Pj et Pj. Par exemple, si 
Ei = Ej, l'exponentielle vaut 1. En tout etat de cause, la contrainte a x <^L"j c n'est 
jamais perdue car la condition Amfj/E^cTp est necessaire lors de l'integration de 
l'amplitude sur p°, sinon le pole correspondant a une des particules i,j se trouverait 
hors du contour d'integration. On peut done toujours effectuer la substitution 



(5.25) 




4 Le terme en T est omis car il ne donne pas d'information interessante. 
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tout en gardant en memoire les contributions possibles de e D a la decoherence. 

Cette condition resout aussi la question de savoir si une mesure sufnsamment 
precise de l'energie et de 1'impulsion des etats oscillants pourrait determiner quel etat 
propre de masse se propage. Une mesure a la precision de Am| viole la condition 
ci-dessus et l'oscillation ne sera plus observable. 

Elle montre aussi pourquoi on ne pourrait observer l'oscillation de leptons 
charges, mis a part le fait que leurs interactions electromagnetiques les rendent plus 
facilement identifiables et pourraient supprimer l'interference. La difference entre 
leurs masses est trop importante et viole la condition ci-dessus. De plus, la longueur 
d'oscillation est si petite qu'elle serait de toute fagon inobservable. Par exemple, 
la longueur d'oscillation du systeme muon-electron est d'environ 10~ 14 m pour une 
impulsion typique de 100 MeV. 

Finalement, notons que si le neutrino mu ou tau est trop lourd, il ne pourra pas 



osciller avec les neutrinos legers [|6Tf| . Soit une source de neutrinos constitute par un 



noyau radioactif dans un reseau. L'incertitude sur la position du noyau est d'environ 
IQ-io m L 'oscillat ion de neutrinos sera done inobservable si Aro^ est superieur a 



4 x 10 9 eV 2 (E v /lMeV). Si le neutrino le plus le ger a une masse negligeable, la masse 



du neutrino le plus lourd ne peut depasser 60keV J E u /lMeV 



5.5.5 Oscillations sophist iquees et oscillations simplifiees 



Sous quelles conditions retrouve-t-on la formule d'oscillation ( |1.7|) habituelle- 
ment utilisee dans la litterature ? Par hypothese, e = Am^/P est un parametre tres 
petit. On peut done developper tous les facteurs precedant les exponentielles en e 
et ne garder que les termes d'ordre zero. Ces termes sont independants des masses 
rrii, rrij et peuvent done sortir de la somme sur Si l'on insere dans la formule 
(|5.14j) de la probabilite les expressions des exponentielles e A , e B , e c et e D (equations 



(15U5|) , flOg ), (pi) et (|Q5D ), on obtient 



V(a^(3, Le z )~Y: V ia Vtf V* a V^* exp (-T tJ - - 2nr-M 
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Les definitions pour le melange ij de la largeur T^-, de la longueur d'oscillation L°j c , 
et de la longueur de coherence L^° h sont donnees respectivement par les equations 



(15121) , (HID ^ CT) . 

II s'ensuit que si les conditions d'observabilite des oscillations sont satisfaites[], 
c'est-a-dire L<^L coh et AmyCcTp, seuls les facteurs de disintegration et d'oscillation 
ont une contribution fortement dependante en La probability d'oscillation s'ecrit 
alors 



V(a^(3,Le z 



hi 



' L L 
—Tij 2in- 



T OSC 



(5.27) 



Le signe d'egalite provient de la condition de conservation de la probability dans le 



cas d'une evolution unitaire, oil T 



et V~ l = : 



P 



P,Le z 



ce qui est bien le resultat attendu : la somme des probabilites de transition d'une 
particule stable dans les autres particules du melange doit etre egale a l'unite. 

La longueur d'oscillation que nous avons derivee (equation ( |5.17| ) correspond 
done a celle de la formule derivee dans le traitement de mecanique quantique (equa- 
tion (|1.10|) ), mais il n'y a plus de tour de passe-passe ! Les reponses fournies par notre 
formule aux objections soulevees contre la derivation traditionnelle sont passees en 
revue dans la conclusion. 



5.6 Oscillations de fermions 

Quelles sont les modifications a apporter pour les oscillations de fermions ? Les 
paquets d'ondes representant les fermions contiendront des spineurs. Le propagateur 
de la particule intermediaire est le propagateur fermionique. Le propagateur libre 



5 Si les particules sont instables, il faut evidemment aussi que L°j 5C < v/Tij, c'est-a-dire que 
les particules aient le temps d'osciller avant de se desintegrer. Cette condition est incluse dans 



l'equation (5.27) 
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d'un etat propre de masse de type Dirac ou Majorana (si la transition ne viole pas 
le nombre leptonique total) s'ecrit 

d 4 P 1> + m „-ip-(x'-x) 



G ^- X ) =i J (27r) 4 p2 _ m 2 + ^ * 

et decrit la propagation d'un fermion de masse m de x en x', ainsi que d'un anti- 
fermion de meme masse de x' en x. Comme pour les particules scalaires, soit la 
particule, soit l'antiparticule a une contribution negligeable a l'amplitude de propa- 
gation lorsque la source et le detecteur sont separes par une distance macroscopique. 
Si l'on s'interesse a des fermions stables dont le lagrangien est connu, par exem- 



ple des neutrinos, la matrice de masse peut etre diagonalisee J41], ^0] par une trans- 
formation unitaire sur les champs. Dans le cas de termes de masse de Dirac, qui 
melangent les champs de chiralite gauche et droite, la transformation unitaire est 
differente pour les champs gauches et droits {M diag = U^MUr) mais seule la ma- 
trice de melange Ul des champs gauches apparait dans les amplitudes de transition. 
En effet, dans le modele standard, les courants droits n'interagissent qu'en tant que 
courants neutres et le terme les contenant est invariant sous une transformation uni- 
taire, contrairement au terme contenant les courants charges. Si le lagrangien inclut 
des termes de masse de Majorana, on peut considerer les termes de masse de Dirac 
comme une superposition de termes de masse de Majorana et les regrouper avec les 
termes de masse de Majorana deja presents. La matrice de masse correspondante est 
toujours symetrique et peut done etre diagonalisee par une transformation unitaire 
U en une matrice a entrees reelles positives {Mdiag — U l MU). La matrice de melange 
U relie aussi bien les champs de saveur (appartenant a des doublets sous SU(2)l) 
que les champs steriles (singulets sous SU{2)l) aux champs de masse. Comme les 
champs steriles ne subissent pas d'interaction faible, seule la partie de la matrice 
U reliant les champs de saveur aux champs de masse interviendra dans les termes 
d'interaction. Cette sous-matrice de la matrice U apparaitra aussi dans les courants 
neutres car elle n'est pas necessairement unitaire [pifl. 

Les amplitudes de propagation pour des neutrinos de Majorana ou de Dirac 
ne different que par des termes en m u /E |82] qui ne sont pas detectables dans les 
oscillations, pas plus d'ailleurs que les phases supplementaires violant CP dans la 
matrice de melange des neutrinos de Majorana, en raison d'une invariance bien 
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connue de la probability d'oscillation sous une reparametrisation de la matrice de 
melange f83fl . 

Si l'on s'interesse a des fermions instables, il faut d'abord calculer l'energie pro- 
pre puis obtenir le propagateur complet par sommation sur les diagrammes d'energie 
propre, en passant par l'equation de Dyson (comme pour le propagateur matriciel 
d'un melange de particules) : 

G(P 2 ) ' Z 



j) — m — f(]/)) + ie 

Le propagateur complet pour un melange de particules se diagonalise similairement 
a celui des particules scalaires a l'aide de matrices V et V~ x . Si le melange entre les 
particules est faible, les parties reelles des poles du propagateur seront approxima- 
tivement egales aux masses avant le calcul de l'energie propre et seront positives. 

L'evaluation de l'amplitude de propagation ne differe pas notablement de celle 
d'un melange de particules scalaires. Les termes contenant des matrices 7 M , sand- 
wiches entre les spineurs des etats entrants et sortants, sont incorpores dans la 
fonction de recouvrement des paquets d'ondes. 



Contrairement a ce que l'on dit souvent |>3|, |84| , il n'est pas necessaire de prendre 
la limite relativiste pour retrouver la formule d'oscillation classique. II suffit que les 
differences entre les masses des etats intermediaries soient beaucoup plus petites 
que leurs impulsions. De plus, il n'est pas necessaire non plus que l'interaction soit 
chirale, sauf si l'on tient a factoriser explicitement le processus j83] en une amplitude 
de production fois une amplitude de propagation fois une amplitude de detection 
(dans ce cas, il faut aussi prendre la limite relativiste). 



5.7 Conclusion 



L'etude des oscillations de particules en theorie quantique des champs a deja ete 
entreprise dans plusieurs articles, dans les cas specifiques des kaons [^, |3| et des 
511641, 



neutrinos |53|, |6J, |63|, g5], |6|, |84], g7|. Ces articles different entre eux par les modeles 
idealises qu'ils proposent et par les approximations utilisees pour le calcul. Aucun 
modele parfaitement fidele et aucun calcul parfaitement exact ne sont possibles, 
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malgre les affirmations que Fori lit parfois ici et la. Certains auteurs prennent trop 
au serieux les modeles simplifies et les methodes de calcul qu'ils ont choisis de sorte 
que leurs conclusions sont a prendre avec un grain de sel. 

Notre analyse est la premiere a traiter de fagon unifiee les cas stable/instable et 
relativiste/non relativiste. Cette approche conduit a une serie de conclusions, parfois 
nouvelles, que nous enumerons ci-dessous. 

1. Les corrections non exponentielles a l'amplitude de propagation ne sont pas 
plus apparentes pour des particules melangees que pour des particules isolees. 
La demonstration n'en avait pas encore ete faite, a ma connaissance, dans le 
cadre de la theorie des champs, bien que des etudes a ce sujet existent en 
mecanique quantique W% |72[ . 



2. Une deuxieme conclusion, deja connue par ailleurs |53], plj| , est que le traite 



ment du processus de propagation dans sa totalite permet de calculer directe- 
ment une probability dependant de la distance et non du temps de propagation. 
L 'application de la formule aux experiences est ainsi possible sans passage am- 
bigu du temps a l'espace par une formule de physique classique exterieure au 
formalisme. Cette ambiguite a conduit a certaines conclusions fausses dans la 
litterature p7| , j50| , |51| concernant les oscillations de particules. Nous y revien- 
drons dans le prochain chapitre, lors de l'etude du processus 0(1020) — > K°K°. 

3. Un troisieme resultat qui a d'abord ete etabli par des arguments intuitifs en 
mecanique quantique |74|, 0, |61], |7^] , puis en theorie des champs |53| concerne 
les conditions d'observabilite des oscillations. D'une part, la longueur d'oscil- 
lation doit etre superieure a l'incertitude sur les positions de la source et du 
detecteur. D'autre part, les oscillations des particules stables ne peuvent etre 
detectees au dela d'une certaine distance appelee longueur de coherence. Une 
nouveaute de notre analyse est la prise en compte de la dispersion, qui sup- 
prime la decoherence pour des particules stables non relativistes. Malheureuse- 
ment, les oscillations etudiees pour l'instant (kaons, B , neutrinos) ne rePevent 
pas de cette categorie. 

4. Notre etude a eclairci la question de l'egalite ou non des energies et/ou im- 
pulsions des etats oscillants. Cette question est en fait mal posee, car il y a 
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plusieurs fagons d'identifier a posteriori les energies-impulsions de ces etats. II 
serait tentant dans notre approche de les faire correspondre a Ej = yj m| + P| 
et Pj, ou Pj est la position du maximum de la fonction-poids lors de revaluation 
de l'integrale (|5.4|). Ces energies- impulsions sont egales ou differentes selon les 
valeurs de Pi et Pj. Notons que Ei — Ej ~ Pi~Pj ~ O(Amy). Cette difference 
etant tres petite, il est delicat d'imposer une valeur bien precise a Ei, Ej, Pj 
et Pj par des evaluations approximatives d'integrales. 

D'ailleurs, Identification que nous venons de proposer n'est pas univoque. 
Modifions un peu la procedure de calcul en prenant la norme au carre de 
l'amplitude avant d'integrer sur la tri-impulsion. On peut alors choisir d'effec- 
tuer d'abord la moyenne sur le temps T avant d'integrer sur les tri-impulsions 
p et p'. On obtiendra en tres bonne approximation l'egalite 



mf + p 2 = ^m) + p' 2 . 

La precision a laquelle est satisfaite l'egalite depend de l'intervalle d'integration 
de T mais cette precision est meilleure que Arriij pour des temps macro- 
scopiques. Cette egalite imposera que l'integration sur p et p' devra se faire 
avec la contrainte Ei = Ej. De cette fagon, on sera conduit a la conclusion 
que les energies des etats oscillants sont egales. Ce raisonnement s'applique a 
d'autres calculs que le notre, par exemple Q . 

II se fait qu'il n'est pas necessaire de repondre a cette question d'egalite ou non 
des energies-impulsions pour calculer sans ambiguite le facteur d'oscillation, 
puisque celui-ci est independant au premier ordre en e = Am^ / Pj des valeurs 
de E h Ej, Pi et Pj : 



Cette compensation entre energie et impulsion se fait parce que les etats pro- 
pres de masse intermediaires sont sur leur couche de masse (dans le sens de 
l'equation (|4.7|) ). lis ne peuvent en effet franchir une distance macroscopique 
que s'ils sont quasi-reels. 

. Pour le cas de particules intermediaires stables (concretement des neutrinos), 
des etudes assez completes ont deja ete realisees |53|, EM. Notre traitement 



114 



CHAPITRE 5. OSCILLATIONS EN THEORIE DES CHAMPS 



montre de plus que 

- la limite relativiste n'est pas necessaire, seule la condition Am^/P <C 1 
compte. 

- il faut tenir compte de la dispersion pour etudier la longueur de coherence. 

- Pegalite ou non des energies et/ou impulsions associees aux particules in- 
termediaires est un artefact de calcul (voir point precedent). 

6. Notre approche a abouti a un seul type de comportement asymptotique en L de 
l'amplitude, contrairement a Particle [57], qui derive un regime intermediaire 
pour l'amplitude, sous la condition a 2 £ P ^> L (qui est probablement impos- 
sible a satisfaire pour des distances macroscopiques) . Dans ce regime, toute 
la dependance angulaire de l'amplitude est contenue dans le facteur oscil- 
lant exp(iPj ■ L) avec la consequence illogique que l'amplitude a le meme 
ordre de grandeur dans toutes les directions, quelles que soient les impul- 
sions entrantes. En comparaison, notre amplitude (equation ( |5.8|) ) contient, 
en plus du facteur oscillant, une gaussienne avec une dependance angulaire 
en exp (— (L — VjT) 2 /4cr^). Cette gaussienne rend notre amplitude negligeable 
dans toute direction differant notablement du vecteur Vj relie directement aux 
impulsions entrantes. L'absence de ce terme dans l'amplitude figurant dans 
|87| n'est pas etonnante car la methode d'integration utilisee dans cet article 



recourt a une methode incorrecte d'integration sur des angles complexes. 

7. Pour le cas des particules instables, les analyses existantes sont incompletes, 
ne tenant pas compte soit de la possibilite P, ^ Pj |32fl , so it de E,i ^ Ej |66| . 
(les auteurs de ce dernier article se restreignent a une dimension spatiale et 
calculent d'ailleurs des integrates de contour non convergentes) . Ces articles 
n'etudient pas non plus les conditions d'observabilite des oscillations. Dans le 
cadre du traitement unifie stable/instable, nous avons discute de la longueur 
de coherence associee aux particules instables, tout en constatant qu'elle est 
en general bien superieure a la distance typique de disintegration. 

8. Quelques lacunes restent a combler. L'effet du temps de vie fini d'une source de 
neutrinos au repos a ete considere dans [^l], |6^[ , mais seulement en mecanique 
quantique. Une etude des oscillations de neutrinos resultant de la desintegra- 
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tion d'un pion en mouvement a ete tentee en theorie des champs |63| mais 
n'est pas convaincante. Bien qu'il soit aise en theorie des champs d'ecrire 
l'amplitude d'un processus complexe (par exemple en cascade), revaluation 
des integrates constitue souvent une barriere infranchissable. Par exemple, la 
description d'une oscillation resultant d'une source instable au repos ou en 
mouvement implique une integration sur deux propagateurs (celui de la source 
et celui de la particule oscillante) dont il est difficile de tirer des informations 
sans effectuer des simplifications abusives. 

II serait dommage de ne pas appliquer la formule d'oscillation que nous avons 
derivee a quelques cas represent at if s. Cette formule devra etre adaptee a chaque 
cas. Pour ce faire, quelques prescriptions de calcul seraient bienvenues pour eviter 
d'encombrantes manipulations de la formule de la probability. C'est le sujet du 
prochain chapitre. 
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Chapitre 6 
Applications 



6.1 Prescription de calcul et 
matrice de masse effective 

La theorie des champs nous a donne les moyens de deriver la probabilite de 
detection d'une particule en melange se propageant sur une distance macroscopique. 
Cette formule d'oscillation n'est cependant pas tres commode a utiliser, meme si 
Ton suppose que les conditions d'observabilite sont satisfaites. Par exemple, l'ap- 
plication de cette formule aux oscillations simultanees du K° et du K° produits 
par la disintegration du meson 0(1020) necessite un retour a l'amplitude pour y 
inclure le propagateur de la seconde particule oscillante. Sous certaines conditions a 
definir, les calculs ulterieurs ne sont pas compliques mais la formule resultante pour 
la probabilite est extremement peu pratique a manipuler. 

Le phenomene d'oscillation est pourtant deja identifiable dans l'amplitude. 
La somme dans l'amplitude d'exponentielles oscillant selon des energies-impulsions 
legerement differentes est a l'origine des interferences oscillant selon les petites 
differences d'energies-impulsions. Pourquoi devrions-nous calculer la probabilite? 
Le probleme git en la dependance de l'amplitude en le temps de propagation, 
dependance qui ne disparait qu'apres integration de la probabilite sur le temps. II 
serait pratique de travailler sur une amplitude dependant uniquement de la longueur 
L (mesuree dans la direction ou la probabilite est maximale), comme celle derivee 
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dans le modele simplifie, a la section |3.4| . Veriflons done si les approximations ap- 



pliquees a la formule de la probabilite ne sont pas transferases a la formule de 
l'amplitude. 

Tout d'abord, si Ton suppose que les conditions d'observabilite de l'oscillation 
(decoherence et localisation des interactions) sont verifiees, seuls les facteurs d'oscil- 
lation et de desintegration restent a l'interieur de la somme sur les etats propres 
dans l'amplitude, comme on l'a ecrit dans l'equation (|5.27| ). 

Ensuite, la comparaison de la formule de l'amplitude partielle Aj de propaga- 
tion d'un etat j (equation ( |5.8j )), avec la formule de la probabilite integree sur le 
temps, V{a j3, Le z ) (equation ( p.!3| )), montre qu'en ce qui concerne les facteurs 
de desintegration et d'oscillation, l'integration sur le temps revient a remplacer T 
par L/v dans l'amplitude Aj. La substitution T^L/v change l'expression de l'am- 
plitude partielle en 



On peut, de plus, approximer mjTj/Ej par mTj/E car la relation Am/m<Ar/r 
est toujours verifiee, puisque r/mCAr/Am. 

Par consequent, on peut etablir la prescription suivante : si les conditions 
d'observabilite des oscillations (decoherence et localisation des interactions) sont 
verifiees, l'amplitude de propagation, dans la direction z oil elle est maximale, prend 
la forme 






A(a->(3, Le z ) ~ Y, v m ex P ~ { 



j 



( 




(6.1) 



On retrouve facilement la longueur d'oscillation dans la difference de deux facteurs 
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oscillants : 

/ m 2 + m 1 m 2 + m 2 \ L Am?. L L 
i ■ _ _ _ . _ i _ 1 _ - 2tt — 



\ 2E 2E v 2P L°f 

Si Ton definit M comme une matrice diagonale dont les termes diagonaux sont 

Mjj = (m 2 + m 2 ) /2E - imTj/2E , 

on peut reecrire 1' amplitude sous les formes matricielles equivalentes : 

A(a-*f3,Le z ) ~ (V' 1 exp {-iM L/v} v) 

~ (expj-iV^MVL/v}) ^ 

~ (exp{-iM saTCUr L/u}) /3a (6.2) 

ou Ton a defini la matrice M saveur = MV. 

Comme m 2 + m 2 = 2mm j, les elements de M deviennent, pour des particules non 
relativistes, 

Mjj = mj - iVj/2 . 

Ces elements sont identiques a ceux trouves lors de la diagonalisation de la matrice 
de masse effective derivee par la methode non relativiste de Wigner- Weisskopf [T3|] . 
Dans la limite non relativiste, la matrice M saveur est la matrice de masse effective. 
Notre calcul prescrit que l'extension relativiste de cette formule se fait en remplagant 
le temps T par le temps propre commun aux deux particules r = mT/E. Dans 
les sections suivantes, nous appliquerons la prescription de calcul ( |6.1| ) ou ( |6.2| ) aux 
experiences CPLEAR et DA$NE qui ont ete decrites pour la premiere fois en theorie 



des champs, mais de fagon differente, dans [48 . 



6.2 L'experience CPLEAR 



Dans l'experience CPLEAR @, (54], [JJ, [5J, une particule K° ou K° est 
produite au point x P par l'annihilation par interaction forte de pp, et se desintegre 
au point x D par interaction faible. Les mecanismes de production du K° ou K° sont 
pp— » K°K~it + , K°K + 7c~. Le rapport de branchement total est de 0.4% et 10 6 an- 
tiprotons par seconde sont diriges vers la cible d'hydrogene. Comme les interactions 
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fortes conservent l'etrangete, les amplitudes suivantes sont nulles : 

M(pp^K°K + n-) = M(pp^K°K-7r + ) = 0. 

L'etrangete des kaons neutres produits peut done etre connue en identifiant le kaon 
charge associe. Apres leur production, le K° ou K° oscille entre ses deux composantes 
Kl et Ks avant de se desintegrer. Nous voudrions decrire revolution spatiale de 
l'amplitude de desintegration et ses phenomenes d'interference. Notons que malgre 
que les kaons charges et les pions contenus dans l'etat final aient des temps de vie 
comparables au Kl, ils seront traites comme des etats asymptotiques. L'invariance 
CPT donne aussi 

M(pp^K°K-n + ) = M*(pp^WK + ir-) = C . 

L 'impulsion moyenne des kaons neutres est de 550 MeV, de sorte que la distance 
moyenne parcourue par un Ks pendant une vie moyenne est de 3 cm. Le detecteur 
permet d'etudier les disintegrations pendant 20 vies moyennes du Ks- Environ 3% 
des Kl se desintegrent dans cette region. 

6.2.1 Disintegrations semi-leptoniques 

Considerons d'abord des etats finaux contenant un pion et des leptons. II y a 
en principe quatre disintegrations semi-leptoniques pour les kaons : 



K°~ 


-> e + 7r~ 


V 


K*- 


-> e _ 7r + 


V 


K°- 


-> e _ 7r + 


V 


K°- 


-> e + 7r" 


~v 



Les deux premieres sont caracterisees par AS" = AQ et sont permises dans le modele 
standard tandis que les deux dernieres sont caracterisees par AS = —AQ et sont 
negligeables dans ce modele car elles n'apparaissent qu'au second ordre dans l'in- 
teraction faible. Notons cependant que cette regie n'est testee experimentalement 
qu'a 10~ 3 pres La desintegration des kaons neutres n'etant pas immediate, 
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les oscillations introduisent une composante de K° dans le K° et vice versa et 
les disintegrations interdites se produisent. En comptant les electrons produits par 
les K° et les positrons produits par les K° (wrong-sign leptons), on peut calculer 
l'asymetrie At(L) definie par 



\M (k°{x 


= o)- 


-» e + 7r v[x = L)\ 


| 2 - \M (K°(x = 0) - 


-> e 7r + z/(x = L)) 2 


\M (K°(x 


= o)- 


-> e + 7r~z/(x = L)\ 


| 2 + |A4(K°(x = 0) - 


-> e~7r+z/(x = L)) 2 



A T (L) = )== 

' ' " = uj -> e+7r-z/(x = LJJ | z + |A1 (A u (x = (J) -> e-7r+z/(x = 

(6.3) 

Les processus dont on calcule l'asymetrie sont relies par la transformation CP. Par 
consequent, une asymetrie At non nulle signifie une violation CP. Cette asymetrie 
est souvent interpretee [^4], Q comme parametrisant la violation du renversement 



du temps T, ce qui revient au meme dans notre formalisme ou CPT est conserve. 
En effet, les violations de CP et de T ne peuvent etre differenciees que si la violation 
de CPT est aussi possible. En mecanique quantique, on se permet de parametriser 
simultanement les violations de CP, T et CPT. Proceder de meme dans notre modele 
mettrait en jeu sa coherence, puisqu'il est base sur la theorie des champs dont le 
theoreme CPT est un pilier. 

Enfin, notons que l'invariance CPT impose que 

M (K°(x = L) -> e + 7T-u(x = L)) = M* (K°(x = L) -> e -n + u{x = L))=V. 

Calculons d'abord l'amplitude %o(L) du processus complet de la production 
d'un K° se desintegrant apres propagation (et oscillation) sur une distance L en 
e~7r + zA La saveur a a la source est K° (S = +1). La saveur (3 du kaon juste avant sa 
desintegration est identifiee par le signe du lepton charge produit ; dans ce cas-ci, on 
veut j3 = K° (S = —1). Notre prescription de calcul ( |6.2| ) donne, dans la direction 
ou l'amplitude est maximale, 

' C 
/ 

Le vecteur colonne a l'extreme droite selectionne la saveur a = K°, tandis que le 
vecteur ligne a l'extreme gauche selectionne la saveur j3 = K°. La matrice de masse 
effective M est donnee par 

, mm s /E-imT s /2E 
M = | ) (6.4) 

mm L / E — imT l/2E 



Tqo(L) ~ (0 V*) V- 1 exp [~ iM ^] v y 
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ou m = (m s +m L )/2 et E est l'energie du systeme e~TT + u resultant de la desintegration 
du K°. Les matrices V et V' 1 ont ete calculees au chapitre || (equations (|3.10| )) et 
sont donnees par 

V = 

v- 1 = = [ ; ; i [ : ; i (6.5) 

L 'amplitude devient 
CD* 





exp 



Apres injection de la matrice M, l'amplitude devient 




1 





exp ( ( -im s -y)^ 1 ) _6XP {{~ imL ~y)t L 



6.6) 

ou P est l'impulsion totale du systeme e~n + h / resultant de la desintegration du K°. 

Calculons maintenant l'amplitude % (L) du processus complet de la production 
d'un K° se desintegrant apres propagation (et oscillation) sur une distance L en 
e + -n~v. La saveur a a la source est K° (S = —1). La saveur (3 du kaon juste avant 
sa desintegration est identifiee par le signe du lepton charge produit ; dans ce cas-ci 
on veut (3 = K° (S = +1). Un calcul analogue donne 

r-m c * vi + i 



2 1 



oxp ( ( -im s - y) ^ L ) - ex P ((~ im L ~y)j> L 



(6.7) 

On peut enfin calculer l'asymetrie At(L). L'insertion des expressions de %o(L) 
(equation (|6.6|) ) et Tq (L) (equation (|6.7|) ) dans la definition de A T (equation ( |6.3|) ), 
donne une asymetrie independante de la distance : 



1 - \a 



ou 

°~ (6.9) 
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II y a violation CP si le parametre a n'est pas une phase pure : |er| ^ 1. La condition 
correspondante pour e s'obtient a l'aide de la relation 



Wee = 



1 - 


a 


2 


1 + 2Wea + 


a 


2 



L'asymetrie At est non nulle si le parametre e n'est pas purement imaginaire : 
Wee 7^ 0. On parle de violation CP indirecte pour souligner le fait qu'elle est unique- 
ment due au melange et qu'elle n'apparait pas dans les amplitudes de production 
et de disintegration C et T>. L'experience CPLEAR donne Ax = (6.6 ± 1.6) x 10~ 3 
PI- 



6.2.2 Disintegrations pioniques 

Considerons maintenant des etat finaux qui sont des etats propres sous CP, en 
particulier des etats constitues de deux pions charges. Ces etats sont pairs sous CP. 
Calculous d'abord l'amplitude associee au processus complet de la production d'un 
K° se desintegrant ensuite en 7r + 7r~. La saveur a a la source est K° (S = +1). La 
saveur (3 a la detection n'est pas identified. On detecte des etats finaux 7t + tt~ qui 
peuvent etre produits aussi bien par la disintegration de la saveur (3 = K° (S = +1) 
que de la saveur (3 = K° (S = —1). On doit done sommer sur la saveur finale f3. 
En utilisant notre prescription de calcul (|B~2|) , l'amplitude totale dependant de la 
distance, dans la direction ou elle est maximale, s'ecrit 

T + _(L) = $>(cw/3,Le z ) 

P 

~ (M{K g ^u + it-) M(K*^ti + ii-)) V- 1 expj-zM^j V 

L'indice +- de T + _ symbolise la charge des pions resultant de la disintegration du 
K°. Le vecteur colonne a l'extreme droite selectionne la saveur a = K°. Les ma- 
trices V et V~ x figurent a la section precedente (equations |6~5|) ), ainsi que la ma- 
trice M (equation ( |6.4|) ), ou E represente maintenant l'energie des deux pions de la 
disintegration du K°. Le lien entre la base des vecteurs propres de l'opirateur CP 
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et la base de saveur a ete donne chapitre ^] (equation (|3.8|)) et s'ecrit 






L' amplitude devient 

% (L) ~ -j=t^ 

+ V ; v/2(l - e 2 ) 

x («M(#i-»7r + 7r) M(K 2 ^k + tj-)) | \ | Jexpj-iM^Jf 1 "Ml 



Elle peut se mettre sous la forme 
C 

72(1 + 



T + _(L)~ C - MiK^TT+n-) (6.10) 



l+X+-e) exp ((-im s ~ y ) ~p L ) + + exp ff-im £ - 
ou P est la norme de l'impulsion totale des pions finaux et le parametre 

X+ ~ = mIkx^tt+tt-) 

decrit la violation CP dans les disintegrations, puisque CP{K\) = CP (2tt) = +1 et 
CP(K 2 ) = —1. On parle de violation CP directe. 

Considerons maintenant le processus analogue de production d'un K° se de- 
sintegrant ensuite en -k + tt~. En suivant la meme methode que dans le cas de la 
production d'un K°, on obtient l'expression suivante pour revolution spatiale du 
K° : 

~ ^(l-e) ^^ 1 "^ 1 (6,11) 



x 



,1 + x+-e) exp ( ( -im s - — j — Lj - (e + x+-) exp \[-zm L - —J — L 



Notons que si nous etions interesses par le mode de disintegration des kaons neu- 
tres en 7r°7r°, il suffirait de remplacer, dans les equations ci-dessus, A4 (K\ ^7r + 7r~) 
par M.[Ki ^7r°7T ) et x+- P ar Xoo; ce dernier parametre etant defini par 

_ M(K 2 ^7l°7T°) 
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A partir des expressions de T + _(L) et de T + _(L), on extrait les rapports mesurables 
entre l'amplitudes de disintegration du K L violant CP et l'amplitude de disintegra- 
tion du K s conservant CP : 

_ M(K L ^n+n~) _ e + X+ - , , 

% --^(K s ^7T+7r-)-l + X+ _e' [bAZ) 

ainsi que 

_ M{K L ^*W) _ e + Xoo 

Vo °- M(K s ^nW)- l + XQO e- [bA6) 

Les parametres rj + _ et ^ 00 sont couramment utilises pour parametriser la violation 

CP dans les disintegrations a deux pions du K L (voir par exemple p. 107 de p9| ). 

Les relations ci-dessus entre les quantitis mesurables et les parametres quantifiant la 

violation CP directe et indirecte ont iti derivees sans utiliser la symetrie d'isospin. 

En outre, il est possible de montrer explicitement que les parametres r] + _ et r] 00 sont 

independants de la convention de phase etrange choisie pour les kaons neutres, ce 

qui n'est pas le cas pour les parametres e et x+-,oo (voir plus loin). 

L'asymetrie preferee des experimentateurs est definie par 



\r + .(L)\ 2 - 


a 


-2 


r + -(L)\ 2 




a 


-2 


T + -(L)\ 2 



ou a = (1 — e)/(l + e) est le parametre deja rencontre lors de l'analyse des disinte- 



grations semi-leptoniques. A partir des iquations pour T + _(L) et T + _(L), on calcule 
que 

A + _(L) = -2\ V+ _\ ■ - — — — — — , (6.14) 



|2 e (T s -T L )mL/2P 



m, — m 



s- 



ou 7] + _ = \t] + -\ e^ +_ et Am 
Cette iquation se riduit a l'expression devolution temporelle utilisie dans l'analyse 
de la collaboration CPLEAR [p5|, si nous choisissons le repere centre de masse des 
deux pions produits par la disintigration du K° ou K° et que nous substituons le 
temps a l'espace par la relation L = vT. Les quantitis a, \r] + _ \ et + _ s'obtiennent si- 
multaniment par un ajustement de la formule ci-dessus aux donnies expirimentales 



de CPLEAR. La quantiti Am provient d'autres mesures |p9[ . Les quantitis | "?7oo | et 
0oo sont obtenues de fagon analogue |96| . 
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6.2.3 Transposition au systeme B°B° 

Tous ces calculs se transposent sans difficulte aux systemes D°D° et B°B°. 



Regions tout de suite le compte du premier en notant que le modele standard p7 
predit que melange entre D° et D° est extremement faible en raison du petit nombre 
de canaux communs de disintegration qui, en outre, sont fortement supprimes par 
de petits angles de melange. Les D° et D° se desintegreront done avant d'osciller 
notablement. II n'y a meme pas lieu de mentionner la violation CP indirecte. La 
violation CP directe est aussi predite d'un niveau negligeable. 

Passons au B°B°. En premier lieu, en ce qui concerne l'asymetrie semi-leptonique 
A T dans ce systeme, les calculs theoriques montrent qu'elle est proportionnelle 
a m 2 c /m% = O(10~ 4 ), ce qui est bien au dela de la precision experimentale actuelle. 
On s'attend par contre a une violation CP importante dans les disintegrations (vi- 
olation CP directe) |j98| , ce qui est confirme par les premiers indices d'une mesure 
de la violation CP dans ce systeme fTE| . Par consequent, la violation CP dans les 
oscillations (violation CP indirecte) est en general negligee. C'est la situation inverse 
du systeme K°K°. 

En second lieu, en ce qui concerne les asymetries non leptoniques, A + _ et A 00 , 
les formules pour les kaons ne leur sont pas appliquees telles quelles, bien qu'il n'y ait 
pas d'obstacle a le faire. Comme les etats finaux sont de toutes sortes, les parametres 
?7 + _ oo et x+-,oo sont notes rjf et Xf- Lors de l'etude de l'asymetrie At, nous avons vu 
que la violation CP dans les oscillations est estimee theoriquement a O(10~ 4 ), tandis 
que la violation CP dans les disintegrations pourrait etre tres importante puisque les 
trois generations de quarks sont directement impliquees dans des processus comme 
Bd — > J/i/jKs- Des lors, la violation CP indirecte est habituellement negligee et l'on 
pose | cr | = 1 (voir equation (|6.8|) ). Cette condition est facile a appliquer si l'on definit 



un parametre /i/ associe a \f P ar l e meme type de relation qui existe entre e et a : 



Le parametre rjf devient 



_1- Xf M(B° — > /) 

H = TTx/ = m{b° -v 7) • ^ 



/; l + afi f 1 + Zf V ; 
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Le parametre £ j est tres interessant car il peut s'exprimer en bonne approximation, 
dans certains cas (B° d -> J/V>#<?, B° d -> tt+tt-, B° -> D* + D*-, B° s -> etc.), 
uniquement en fonction des elements de la matrice de melange des quarks (matrice 
CKM). Sa mesure permettra de verifier si le modele standard suffit a parametriser 
completement la violation CP. Notons neanmoins que les asymetries CP dependent 
directement des interactions dans l'etat final (qui fournissent les phases conservant 
CP des amplitudes partielles, voir plus loin), dont les contributions a longue distance 
sont difficiles a evaluer. La prediction de ces asymetries varie de fagon significative 
selon le modele choisi pour estimer ces interactions dans l'etat final . 

Comme pour l'experience CPLEAR, on s'interesse a une asymetrie pour des 
etats finaux etats propres de CP, definie par |100|] 

r(B° d (L)->f) + r(B' d (L)->f) 

Les etats propres de propagations ont des largeurs quasiment egales |70| et sont 
differencies par leurs masses. Supposons que l'etat pair sous CP soit plus leger (note 
P>l pour Light) que l'etat impair sous CP (note E>h pour Heavy), comme dans le 
cas des kaons et definissons Am^ = m H — m L . Partant des expressions de T + _(L) 



et T + _(L) (equations Qfi.lUp et ( p.ll] )) et y substituant B pour K et l'expression de 
£f (equation (|6.16|) ), on calcule l'asymetrie A C p avec les approximations \a\ = 1 et 
T H = T L : 

A CP (B° d ^f,L) = Af P ect (B° d ^f) C0S (AmA+ A^B^f) sm(Am d mL 



p j • O-r \— a j i \ " p 



ou la violation CP directe, 



a direct i r>0 . f\ 

A CP K B d , | u 12 = 1 7 |„ 12 ' l 6 " 17 ) 



i-IOI 2 i-W 

a ete separee de la violation CP provenant de l'interference entre le melange et la 
disintegration du B : 

A CP (B d ^f) = - r ^ = - i + M2 . (6.18) 

On voit ici l'importance de tenir compte du melange pour la violation CP meme s'il 
n'y a pas de violation CP indirecte, puisque A™p depend de l'interference entre a 
et /x f . 
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Examinons de plus pres le cas de l'etat final / = J/ipKs dont l'asymetrie vient 
d'etre mesuree. Une etude theorique des differentes contributions au processus (voir 
par exemple fTPf ) montre que le facteur ^j/^k s est en tres bonne approximation egal 



a Cj/4>k s — e 24/3 , oil P est un angle du triangle unitaire (pour la definition de ce 



triangle, voir la reference [69], p. 105). Ces angles sont independants des conventions 
de phase pour les quarks comme il se doit, puisque £f est aussi independant de phase. 
Les asymetries directes (|6.17|) et induite (|6.18| ) deviennent 



A d j P ect {B° d -> J /iff K s ) = et A%* uit (B° d -> J/$ K s ) = - sin 2(3 . 



L 'experience [jnj donne 

A vnduit( B o _^ j/^Ks) = -0.79j££ (STAT+SYST) , 

et constitue le premier indice de l'existence d'une violation CP dans le systeme 
B°B°, ainsi que la premiere etape de la verification de l'unitarite de la matrice de 
Cabibbo-Kobayashi-Maskawa. 



6.3 Production de kaons correles a DA<I>NE 

Cette section est consacree aux oscillations d'une paire de kaons neutres pro- 



duits par l'annihilation electron-positron dans le collisionneur DA$NE a Frascati|88] . 
Les resultats obtenus dans notre formalisme sont immediatement applicables a la 
description de la production de paires de mesons B dans la region de la resonance 

t(4s) n,g. 

Les kaons neutres et charges sont produits en quantite (~ 10 9 paires de K°K° /an) 
dans des collisions electron-positron a une energie de centre de masse dans la region 
de la resonance 0(1020). Le meson produit se desintegre au point x en une paire 
K°K° avec un rapport de branchement de 34%. Les kaons partent en sens op- 
poses avec une faible impulsion et ne parcourent que 6 mm pendant le temps de 
vie moyen du Ks- Cette technique permet de mesurer dans un meme detecteur les 
disintegrations des Ks et des Kl ainsi que leurs correlations : si un Ks est identifie 
par sa disintegration en deux pions, la particule associee est en principe un Kl- 
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Chaque kaon oscille entre ses composantes K L — K S avant de se desintegrer en 
etats finaux f\{k\) et f 2 {k 2 ) aux points y\ et y 2 : 

4>{q) -> -> fi(ki)f2(k 2 ) , 

ou g, fci et &2 sont les energies-impulsions correspondantes. 

Chaque etat final pouvant etre produit soit par K°, soit par K°, il faut sommer 
les deux amplitudes provenant de l'echange de K° et de K° avant de prendre la 
norme au carre. Les nombres quantiques du 0, J PC = 1 , sont conserves par 
les interactions fortes provoquant la disintegration. La paire K°K° se trouve done 
dans un etat antisymetrique sous P, ainsi que sous C. Le signe relatif des deux 



contributions a 0— >fif 2 est done negatif fl7|l . 



Comme deux particules oscillent dans ce processus, notre formule d'oscillation 
n'est pas directement applicable. Considerons momentanement le meme processus, 
mais pour des particules non melangees et non antisymetrisees. L'amplitude du 
processus s'ecrit selon la methode du chapitre 0, section 



A = J [dq\ $ P J [dki] $* Dl J [dk 2 ] $^ A ondes planes 



avec 



A ondes planes (q, h, k 2 ) = J d A x e*** J d A yi M Dl e ik ^ J d A y 2 M Di e lk ^ 

x / 1^ e ~ ipi<yi ~ x) GM) I wt e ~ lpHV2 ~ x) ° M) Mp 

ou Giij>1) symbolise le propagateur de la particule se desintegrant en fi(ki) et les 
particules exterieures sont sur leur couche de masse. Les Mp et Mp. sont les ampli- 
tudes des processus de production et de detection. Les paquets d'ondes $p et 
sont centres dans l'espace des impulsions respectivement autour de Q et Kj tandis 
que leurs transformers de Fourier sont localisees dans l'espace de configuration re- 
spectivement autour de x P et y D . ou x P = (t P ,x P ) et y D . = (t D .,y D .) avec i — 1,2. 
Dans la notation du chapitre ^ (voir equation (|4.1| )), ils s'ecrivent 



$ P (q,Q,x P ,t P ) = P (q,Q)e^ 
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Apres le changement de variables x^x + x P , Vi^Vi + y Di et les integrations sur les 
variables x et y i7 1'amplitude devient 

A = J J d% <p(p U p 2 ) G 1 ( P l) G 2 (pj) e-^ivn^-miyn.-xp) ; 

ou la fonction-poids ip(pi,p 2 ) est une integrate de recouvrement des paquets d'ondes 
entrants et sortants. Elle est definie par 

V(p u p 2 ) = y , [dq]0 P (q,Q)|[dki]^ 1 (k 1 ,K 1 )y , [dk 2 ]^(k 2 ,K 2 ) 

x (2tt) 4 6®(q - Pl - p 2 ) ^ 4 )( Pl - A*) 5^\p 2 - fc 2 ) M P M Dl M Di . 

Les integrations sur p\ et se f° n t P ar integrale de contour, selon la methode etablie 
au chapitre f|, section [0| II s'ensuit que 

A « /" rf 3 pi / rf 3 p 2 V'(-2l,'2 ; 2,Pl,P2) 



x exp ( — iy -2i + Pi Ti + ip! ■ Li - i\Jz 2 + p| T 2 + zp 2 ■ L 

ou Tj ee i D . — i P , Lj ee y D . — Xp et les Zj sont les poles des Gi{pj). 

La norme de la fonction ip est maximale pour pj = Kj et pi + p 2 = Q et 
les Lj alignes avec Kj. L'evaluation des integrales pose des difficultes a cause de la 
correlation entre pi et p 2 , sauf si la contrainte pi + p 2 — Q est imposee avec moins 
de precision que les contraintes pj = Kj. C'est justement le cas pour <p — > K°K° 
puisque la relation Pi + p 2 = Q est exacte a = 4.4 MeV pres, tandis que les 
relations pj = Kj sont exactes a l'incertitude sur la masse des kaons pres, bien 
inferieure a T^. On peut done effectuer independamment les integrales sur pi et p 2 . 
Le reste du calcul s'opere comme si les particules etaient independantes et notre 
prescription de calcul est applicable. 

En fin de compte, si les directions de Li et L 2 sont fixees, 1'amplitude anti- 
symetrisee de detection de fi a une distance L\ et de f 2 a une distance L 2 est 
donnee pour le systeme melange K°K° par 



IfihiLi, L 2 ) 



x 



M(K®^fi) M(K Q —>f{)} V- 1 exp j-iMx^ij V 
M(K°^f 2 ) M(W^f 2 )) V- 1 exp|-zM 2 ^} V I ° 
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meme expression avec 








Ml 



K°K°) 



La notation signifie que la matrice de masse effective M est evaluee a l'energie Ei 
de l'etat final /j. Les Vi sont dermis comme d'habitude par la norme de Vj = Kj/E'j. 
L'amplitude se reecrit en fonction de la base propre de CP comme 



T hh (L)~M(<l>^K»K») 



1 



2(1 -P) 



x 



(M(K t ^f 2 ) M(Kr+h)) 



meme expression avec 




exp < —iMi 



Li 

Vl 



exp <{ —iM 2 — 

V 2 




Si Ton definit Mij = M.(K i —>fj) ) l'amplitude devient 
T hh (L) ~ M^ ^ K°K*) ^ 
x - (Mn + £M 2 i) (eM 12 + M 22 ) exp (- (im s + 
+ (eMn + M 2 i) (M12 + eM 22 ) exp (- (im s + y ) 



r 5 \ mil 



2 / Ki 



im, 



im, 



2 7 
Tj\ mLi 
2 ) Ki 



L'expression obtenue est relativiste : aucun boost n'est necessaire pour passer du 
repere au repos d'une particule oscillante au repere du laboratoire. Observons que 
l'amplitude est nulle |T7J si fi = f 2 et Ki = K 2 et L\ = L 2 . Pour l'application des 
formules theoriques aux experiences, nous referons a la litterature |9l] . 

Une remarque a propos de la frequence d'oscillation semble de mise. Dans le 
repere centre de masse, Ki = K 2 = K et le terme d'interference oscille comme un 
cosinus de 

m 



— (m L - m s ) (Li - L 2 ) . 
K 

Si Ton reecrit ce terme comme lj(Li—L 2 ), la frequence d'oscillation u vaut uj 
Le calcul de cette frequence a suscite une controverse a la suite d'un article [j3?H 



mAm 
K ■ 
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pretendant qu'elle etait deux fois plus grande. Plusieurs articles EH, EM] ont ete 



ecrits pour refuter cette affirmation mais leurs arguments reviennent a contester le 
choix de differents temps de detection et de temps propres pour les differents etats 
propres de masse. Tout cela est assez obscur. Seule une derivation en theorie des 
champs donne sans ambiguite la formule d'oscillation. La valeur calculee est par 
ailleurs confirmee KJ par les experiences ARGUS et CLEO dans le cas du 



processus analogue Y— >£>°£> . 



Un autre article controverse pi soutient que dans un processus tel que 



la particule A associee a la production de la particule oscillante K° oscillera aussi : 
elle serait une superposition de deux etats d'energie-impulsion differentes par conser- 
vation de l'energie-impulsion a la production. Notre approche montre qu'il n'en est 
rien. Le processus peut etre decrit completement comme ci-dessus, avec un propa- 
gateur pour le kaon et un propagateur pour le A. Le propagateur du kaon est une 
superposition de deux propagateurs correspondant aux etats propres de masse : a 
grande distance, deux poles differents contribueront et l'amplitude sera la somme 
de deux amplitudes correspondant a deux etats propres de masse differents sur leur 
couche de masse. Ce n'est pas le cas du A. On voit clairement dans cet exemple 
que la question ne peut etre tranchee sans le formalisme de theorie des champsQ. 
Comme nous l'avons vu dans l'equation ( |5.16| ), les differences d'energie et d'impul- 



sion se compensent dans le facteur d'oscillation saufla difference de masse entre les 



etats propres. La meme question s'est posee |51|, (44|, |35|, |86|, |92[ pour le processus 



7i— >fiv. Le muon oscillerait-il ? Par le meme raisonnement, on conclut que non. 



^es refutations recourant a la mecanique quantique ||39(1 font de nouveau appel a la notion de 
differents temps propres. 
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6.4 La parametrisation correcte de la violation 
CP 

Le developpement du formalisme theorique adapte a la description de l'expe- 
rience CPLEAR nous a montre que la violation CP peut etre parametrisee dans 
le systeme K°K° ou B°B° par des parametres complexes e et \ s . Cette simplicite 
cache certaines subtilites dans l'extraction des parametres a partir des donnees. 
Tout d'abord, il faut mentionner que ces parametres ne sont pas des observables 
physiques car ils dependent de certaines conventions de phase. Ensuite, des approxi- 
mations assez violentes sont souvent appliquees dans le systeme K°K°, par exemple 
negliger les termes en e 2 , alors qu'ils sont du meme ordre que les parametres Xs- 
Les ambiguites d'ordre 0(e 2 ) dues au probleme de la normalisation des etats K$,l 
sont aussi oubliees. La plupart des analyses ne tiennent pas compte de ces points 
delicats. II est vrai que les experiences ne sont pas encore assez precises pour dis- 
criminer les demarches correctes des mauvaises mais il n'est plus l'heure de passer 
a cote de telles incoherences quand la violation CP directe est confirmee par deux 
experiences independantes (T^|. Par ailleurs, les conventions de phase choisies pour 



le K°K° ne sont pas necessairement pratiques pour le B°B° ; il est des lors utile 
de garder la liberte de phase dans le formalisme. Etant donne que les definitions 
de nos parametres ne souffrent pas de probleme de normalisation, il est interessant 
de montrer que l'analyse peut etre continuee rigoureusement jusqu'au bout dans 
l'extraction, a partir des resultats experimentaux, des parametres violant CP. Nos 
resultats finaux ressemblent a ceux de la reference [ 101|| mais celle-ci utilise des 



etats intermediaires de kaons et ne discute pas explicitement de la transformation 
de phase. Avant de determiner les valeurs des parametres, on etudiera la question 
de la convention de phase. 



6.4.1 Phase etrange et parametres observables 

Une premiere manifestation de la dependance de phase de e est apparue lors 
du calcul de l'asymetrie semi-leptonique At- Une asymetrie non nulle est une mani- 
festation de la violation CP et implique que Ve e ^ 0. Par contre, si Tie e = mais 
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Xme 7^ 0, cette asymetrie s'annule et n'implique pas de violation CP. La partie 
imaginaire de e n'a done aucune importance dans ce cas-ci. Cet arbitraire dans la 
valeur de Xm e est en fait la manifestation d'un phenomene plus general. La raison 
provient de la symetrie U(l) associee a chaque saveur, et qui laisse le lagrangien 
total invariant, excepte le terme de masse (ou la matrice de melange CKM apres 
diagonalisation de la matrice de masse, ce qui revient au meme) . On s'interesse ici a 
la transformation S de la phase etrange ||103j , |102|| , puisque les kaons possedent une 



etrangete non nulle. De nouveaux etats peuvent etre definis a l'aide de S : 

\K° a ) = e~ iaS \K°) = e~ ia \K°) , 
\K°) = e~ ta§ \K°) = e +ia \K*) , 

ou 5* est l'operateur d'etrangete. Cette redefinition est possible parce que K° et 
K° sont produits par les interactions fortes qui conservent l'etrangete. Elle ne peut 
evidemment avoir d'effet observable : toute quantite observable doit etre indepen- 
dante de cette phase. Dans le but de conserver la relation \K°) = CP \K°), la trans- 
formation CP est redefinie dans la nouvelle base : 

(CP) Q = e~ ia§ CPe ia§ , 

de sorte que 

\Kl) = (CP) a \K°) . 
La base propre sous (CP) a est definieQ par : 



V2 

Comme la transformation ne laisse pas invariants les termes melangeant des par- 
ticules d'etrangetes differentes, les elements non diagonaux du propagateur complet 
du systeme K°K° sont modifies. Se rappelant la parametrisation du propagateur 
(equation 



, (K°\G-'\K°) (K°\G-i\K°) 
' 1 (K°\G- l \K°) (K Q \G- l \K°) 




2 Une autre possibilite consiste a ne pas redefinir l'operateur CP. On a alors CP\K^} 
Va F \Ka)> ou r fa > = e _2lQ . II faut definir difFeremment les etats propres sous CP 
\Kl,2a) = js (K)± CP \K°)) 
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on voit que les termes non diagonaux sont transformes sous S en 

a a + b a = (a + b)e 2ia , 
a a -b a = (a- b) e~ 2ia . 

Quel est l'effet de S sur le parametre e ? Ce parametre a ete defini (equation 
)) comme la solution de 



1 + e 2 2a ' 

Cette equation possede en fait deux solutions e\ et €2 inverses Tune de 1' autre : 



a 



e\ 2 = - \1±\1 



fe2 



ou le signe de la racine est choisi de sorte que sa partie reelle soit positive. Comme 
ii = l/e.2, une des deux solutions a necessairement une norme inferieure a 1. On 
choisit de travailler avec cette solution, dans l'esprit de representer la petite violation 
CP observee par des parametres theoriques petits. Nous choisissons done la solution 
avec le signe — : 

Lors du calcul de l'asymetrie At, un parametre a a ete associe a e (equation 
par 

1 - e A I- a 

= e= . 

1+e 1+a 

Son expression en fonction de a et b s'ecrit 



Va 2 - b 2 

a — ; — • 

a + b 

Notons que 

1 — lei 2 

'Rea = , ,0 > 

1 + 2-Ree+ \e\ 2 ~ 

puisque |e| < 1. 

Sous la tranformation S, le parametre a devient 
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Partant de l'equation liant e a et a a , 

1 - rr. 



1 + cr a 

on etablit la loi de transformation de e : 

e + i tan a 



a 

i 



1 + ie tan a 

1 — I £ 1 «in 9w 4- 9 7rn 2 pns 9n/l 

(6.19) 



2 73ee + i[(l - 


e 


2 ) sin 2a + 2 Tm e cos 2a] 


1 + cos 2a — 2 Xm e sin 2a + 




2 (1 - cos 2a) 



Remarquons d'abord que Tie e a oc Tie e. La condition 73e e 7^ est done bien invari- 
ante sous S comme il se doit. Cependant, la transformation S ne peut etre tout a 
fait arbitraire. En effet, le choix de la valeur a = ±7r/2 resulte en la transformation 
e a = 1/e. On doit done restreindre le domaine de a si Ton veut respecter la condition 
|e| < 1. Comme les expressions ne dependent que de 2a, il suffit de travailler sur 
l'intervalle — n/2<a<ir/2. On calcule que |e Q | < 1 pour 

If 1 - le| 2 \ 1 1 - |e|~ 

- £ ' 



arctan — — — - — ir < a < - arctan ■ 



2 V 2Xme J ~ ~ 2 21me 

ou F arctan est pris dans le premier ou deuxieme quadrant. 



6.4.2 Determination des parametres 

Examinons maintenant la premiere donnee experimentale dont nous disposons 
pour les kaons. La petitesse de l'asymetrie semi-leptonique, Aj- = (6.6 ± 1.6) x 10 -3 
P^| , nous dit a travers l'equation ( |6.8| ) que \a\ — 1 + O(10~ 3 ). Comme 



Ttee = 



1 - 


a\ 


2 


1 + 2 Vea + 


a\ 


2 



avec Tie a > 0, on conclut que T&e ~ (10 -3 ). La valeur de TZee depend bien sur du 
choix de la phase etrange mais on peut deja donner son ordre de grandeur qui sera 
correct tant que la condition |e| <1 est maintenue. Par contre, on ne peut rien dire 
sur Tm e sans faire un choix de phase etrange, puisque 

—2Xma 
l7Ui= l + 2-Rea+\a\i- 
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Peut-on choisir la phase etrange telle que \Tme\ <C 1 ? En utilisant l'equation 
(|6.19|) , on obtient 

2 In 6 

Xm e a = si tan 2a = — — . 

1 — \e\ z 

Sif] Tm i > 0, — tt/4 < a < 0. Si Ime < 0, < a < tt/4. Par consequent, il est 
toujours possible de choisir la phase etrange telle que Tme a = 0. Toute une serie 
d'autres choix sont possibles avec \Tme a \ <C 1. 

Poursuivons notre examen des donnees experiment ales. Dans les disintegrations 
des kaons neutres en pions, les observables lies a la violation CP ont ete definis 
(equations ( |6.12| ) et ( |6.13j )) par 

M.(K L — > f) _ e + x f 
Vf M{K s ^f) l+x f i' 
ou Xf = M.{K<i — > f)/M.(K\ — > /) et / est un etat pair sous CP. Comme les 
observables doivent etre invariants sous la transformation de phase S, les quantites 
Xf compensent la variation de e sous S. On le verifie le plus aisement a l'aide du 
parametre \i s associe a Xs P ar l'equation ( |6.15|) . La transformation de fi f sous S est 
tout simplement la transformation inverse de a : 

La quantite £ f = a\i s est done invariante sous S et des lors 77/ aussi (equation ( |6.16| )). 
On a verifie du meme coup que le parametre £ /; decrivant la violation CP dans le 
systeme B°B°, est invariant de phase, comme il se doit puisqu'il est mesurable. Pour 
un usage futur, je donne la transformation explicite de Xs '■ 

Xt — itanct 

X f ,a= , Xf . I • 6.20 

1 — %Xj tan a 

Les parametres r] f = \r) f \ e^f ont ete mesures par CPLEAR |^5|, : 

\r] + _ I = (2.254 ± 0.024 STAT ± 0.026 5yST ) x 10" 3 , 

+ _ = (43.63 ± 0.54 5TAT ± 0.48 syST )° , 

|^ 00 | = (2.47±0.31 STAT ±0.24 SYST ) x 10- 3 , 

0oo = (42.0 ± 5.6 STAT ± 1.9 5ysr )° . (6.21) 

3 On a choisi des valeurs de a appartenant ail domaine tel que |e a | < 1. 
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Des lors, si la convention de phase est telle que |e| ~ 0(7Zei) ~ 0(1CT 3 ), les 
parametres Xs doivent etre au plus du meme ordre : \ s < (9(10 -3 ). La violation 
CP indirecte pour les kaons est done au maximum du meme ordre que la violation 
CP directe. Une difference entre r/ + _ et r] 00 etablit l'existence d'une violation CP 
directe. C'est ce qu'ont fait les experiences NA31, E731, NA48 et KTeV |12| en 
mesurant 

| r?00 |2/| r/+ _| 2 = 1 - 6 x (21.2 ±4.7) x 1(T 4 , (6.22) 

ou l'on a fait la moyenne des resultats ||68|| . 

Dans le cas ou I'etat / est impair sous CP, on definit le parametre 

_ M(K s -+f) _ e + Xf 
Vf Ji4(K L — > /) 1 + X/e' 

ou Xf = M{K\ — > f)/M(K 2 — > f). Si la violation CP directe est petite par rapport 
a la violation CP indirecte, il existe une convention de phase telle que rij = e = rjf. 
On predit par exemple 77000 = ?7oo + O(10~ 6 ), ou l'indice 000 symbolise I'etat a trois 
pions neutres. 

Pour aller plus loin et extraire les quantites e, x+- e ^ Xoo (dependant de 6 
nombres reels) de At, et 77 00 (dependant de 5 nombres reels), il faut disposer 
d'une relation supplement aire fournie par la symetrie d'isospin des pions. 

Cette technique bien connue (voir par exemple ||104|| ), permet de decomposer 
les amplitudes de transition des kaons en pions en amplitudes correspondant a des 
isospins dermis : 

M(K° -tt+O = ^e^Ao + ^e^A,, 
M(K* ^tt+tt-) = fit&AZ + Ji&AZ, 



Mi K [) — Tr'V i - -^e iS °A + J^A 2 , 



M(if°-A°) = -Jl e iS °A* + J^e iS >A*. 
L'etat 7r + 7r~ a ete symetrise et l'on a introduit un facteur l/y/2 pour I'etat final 



0^0 



7T 7T 



pour fournir le facteur 1/2 dans la probability de detection de deux parti- 
cules identiques. La brisure d'isospin est prise en compte dans la parametrisation en 
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redefinissant les parametres Aj par combinaison lineaire. La decomposition garde la 
meme forme bien que la symetrie d'isospin ne soit plus respectee. Soit Aj = \Aj \ e 1 ^ 1 . 
Sous une transformation de phase S, Ai a = e~ ta Aj et il est par exemple possible 
de rendre A reel (convention de Wu-Yang). La difference de phase (Co — C2) reste 
invariante. L'asymetrie A^ definie par 

_ \M(K° -> 7T + 7T-)| 2 - \M(W -> 7T + 7T-)| 2 
+ ~ = |A^(K° -> 7T + 7T-)| 2 + |7W(^0 _, Tr + T,-)^ 

est donnee en terme des amplitudes d'isospin par 









2^2 sm(S - 


-<y 


sin 


(Co - Ci) 1 Ail 


\A 2 \ 




2|A)! 


2 + 


A 2 \ 


2 + 2^^0! 


|^| 


cos(<5 - 


■5 2 ) 


cos 


(Co - C 2 ) 



A+ _ = -v- ~~ v ~ u ^ ^i^u, _ (6>23) 



Elle est non nulle s'il y a une violation CP directe. C'est le cas s'il y a au moins deux 
amplitudes partielles non nulles (A et A 2 ) avec une phase relative (Co — C2) violant 
CP et une phase relative (8 — 5 2 ) conservant CP. L'asymetrie A 00 correspondant 
aux pions neutres a une forme similaire. 

Ces asymetries ne sont cependant pas observables telles quelles puisque les kaons 
se propagent et oscillent avant de se desintegrer, donnant lieu a une violation CP 
indirecte en plus de la violation CP directe. Comme on l'a vu ci-dessus, les quantites 
mesurees sont 77+ _ et i] 00 . On voudrait parametriser independamment les violations 
CP directe et indirecte. Ce n'est pas possible avec la definition de 77/ en fonction de 
e et Xf puisqu'ils ne sont pas separement invariants sous la transformation de phase 
S : meme en absence de violation CP directe, le parametre Xs ne sera nm Q ue pour 
un choix particulier de phase. Pour determiner l'expression de r] + _ en l'absence de 
violation CP directe, calculons-le d'abord en fonction des amplitudes d'isospin : 

V2KeA (e + i«o) + e iS T?eA 2 (e + in 2 ) 
V+ ~ ' \/2 / ReA (l+ieK^+e i&/ ReA 2 (l + ieK 2 ) ' 

ou Kj = Tm AijTZe Ai = tan(j et 5 = 5 2 — 5 . La violation CP directe est par exemple 
absente si A 2 = 0. Dans ce cas, 

e + iK 

V+- = ; , ~ = e , 
1 + te k 

que l'on prend comme definition d'un nouveau parametre, e, invariant sous la trans- 
formation de phase S. On peut verifier explicitement que ik se transforme comme 
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X f (equation ( |6.20| )) et joue done le role de Xs dans la compensation de la variation 
de e sous S : 

iK — itana 
1 — i(iK, ) tana 

D'une part, la condition Kee ^ signale une violation CP dans l'asymetrie A + _ 
car Tie e oc Tie e. Comme elle rend l'asymetrie At non nulle, on l'interprete comme 
une violation CP dans les oscillations. D'autre part, Ime 7^ signale aussi une 
violation CP dans l'asymetrie A + _(L), qui resulte de l'interference entre l'oscillation 
et la disintegration. On pourrait l'appeler violation CP induite, en analogie avec la 
violation CP induite definie dans le systeme B°B° (equation ( |6.18| )). En utilisant la 
definition de e, on peut maintenant ecrire l'expression de rj + _ en presence de violation 
CP directe, e'est-a-dire pour A 2 7^ : 

"^ £ + I7f77r (6 - 24) 

ou e garde sa definition et ou u et e' sont definis par 

1 + ie k 2 TieA 2 iS 

00 = — e , 

1 + ie k Tie A 

£ , = iu(l- e 2 ) (k 2 - k ) 



y/2 (1 + ie k 2 ) (1 + ie k ) 

Les parametres e, uj et e' sont invariants sous la transformation de phase S. 

On precede de meme pour la disintegration des kaons en pions et Ton obtient 

"» = £ -T37S- (6 - 25) 

Ces equations sont exactes, contrairement a la plupart des expressions des r\ s en 
fonction de e et e' trouvees dans la litterature, ou e et e' ne sont meme pas invariants 
de phase. Cette parametrisation exacte nous evite de devoir negliger des termes en 
e 2 qui sont du meme ordre que e' . Dans la litterature (par exemple dans ||102|| ), 
les definitions de uj et e' different des notres et ne sont pas invariantes de phase, 
avec pour consequence que les expressions pour r] + _ et r] 00 , bien qu'identiques a nos 
formules (|6.24|) et ( |6.25|) , ne peuvent etre obtenues qu'apres des approximations 
negligeant e 2 par rapport a e'. 
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Que deviennent les different s cas (voir l'equation ( |b.23p ) pour lesquels l'asy- 
metrie A + _ s'annulait ? Entrainent-ils l'annulation de e', le nouveau parametre 
decrivant la violation CP directe? 

1. Si k 2 = k (c'est-a-dire ( 2 = ( ), alors e' = 0. 

2. Si A 2 = 0, alors e' = 0. 

3. Si A = 0, la decomposition ci-dessus est mal definie. On a plutot 

e + ik 2 2g' 
Q — ~ — et rj + _ = g 



1 + ien 2 1 + ^/2w ' 

ou q' et w sont obtenus a partir de e et u en echangeant les indices et 2 
(done q' = —e/uj 2 et vj — uj^ 1 ). Si A = 0, le parametre g' est nul. 

4. Si 5 = et TZee = 0, alors 7^ee' = mais il est toujours possible que 
2m e' 7^ 0. Le choix de la phase tel que e = ne change rien. On pourrait 
se demander comment il peut y avoir une violation CP alors qu'une des condi- 
tions de son apparition n'est pas satisfaite, c'est-a-dire l'existence d'une phase 
relative conservant CP entre les amplitudes partielles. Contrairement aux ap- 
parences, il n'en est rien ! En effet, les conditions 5 = et e = impliquent que 
7ter] + _ ~ TZee' = 0, ou encore + _ = ±7r/2 (rappelons que r/ + _ = \r] + _ \ e^ +_ ). 

L'examen de l'asymetrie A^ (equation ( |6.14| )) montre que cette asymetrie 

serait nulle sous ces conditions, si ce n'etait la presence de Am ^ dans 1' ar- 
gument du cosinus. La phase relative conservant CP est done fournie par la 
difference de masse entre les etats se propageant. 

Observons aussi que meme si l'asymetrie semi-leptonique At est nulle (TZee = 0), le 

parametre e joue toujours un role dans l'asymetrie A^ a travers sa partie imaginaire 

Xm e. Non seulement Tm e peut se combiner avec k pour donner Xm e ^ comme 
on l'a vu ci-dessus, mais il interfere aussi avec les parametres Aj, £ r et 5 a l'interieur 
de e'. En bref, le parametre de melange e intervient dans tous les types de violations 
CP. Les expressions des r\ f se simplifient si l'on utilise l'information experimentale 
(regie AI = 1/2) que |A 2 |/|A | = 1/22 ||105|| . Les termes en ui peuvent alors etre 
negliges et les equations ( |6.24| ) et ( |6.25| ) prennent la forme 
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Les valeurs trouvees par CPLEAR (voir equations fl6.21| )) etant egales aux erreurs 
experiment ales pres, on sait que e' est beaucoup plus petit que e. Les experiences 
dediees a la mesure de la violation CP directe mesurent plutot Tie (e'/e) : 



\ Voo \ 2 /\ V+ _\ 2 = l-6'Re- + Ol- 



Cette equation est identique a celle trouvee dans la litterature (voir par exemple 
p8|), mais a ete obtenue sans approximations douteuses. Nous avons en effet resolu, 
d'une part, le probleme de la normalisation a l'ordre e 2 des etats oscillants et, d'autre 
part, nous avons defini les parametres uj et e' de sorte qu'il ne nous est pas necessaire 
de recourir a des approximations incorrectes du type e 2 <C e'. La compilation des 
resultats experimentaux (voir equation fl6.22p ) fournit la valeur 

Tie - = (21.2 ±4.7) x 10" 4 . 

Avec une convention de phase telle que |e| <C 1, les expressions de e, u et e' prennent 
la forme approximative (et non invariante de phase) : 

e = e + iK , 
TteA 2 ig 
Tie A 

e' = -^=u(k 2 -k ), 

Dans ce cas, Tiee = TZee et la mesure des r\ s fournit la valeur de TZee. Dans la 
meme convention de phase, l'asymetrie semi-leptonique At devient At = ATiee. 
On peut done comparer les resultats experimentaux provenant de At et de A + _ et 
A 00 . lis coincident aux erreurs experiment ales pres. Dans la convention de phase de 
Kobayashi-Maskawa, Xm A 2 = si la symetrie d'isospin est respectee. Le terme k 2 
dans l'expression de e' parametrise done la violation de l'isospin dans cette conven- 
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Sous les memes conditions, les parametres \s de violation CP directe se reecrivent : 

X+- = e' ± in , 
Xoo = -2e' + iK , 

Dans la convention de Wu-Yang (qui est une des nombreuses conventions telles que 
|e| < 1), k q = done x+- - e' et Xoo = -2e'. 



Conclusion 



Les oscillations spatio-temporelles de la probability de detection dans la propa- 
gation de particules en superposition quantique ont une importance capitale aussi 
bien du point de vue fondamental, dans l'etude de la violation CP dans les systemes 
de kaons neutres et de mesons B et dans la determination du spectre de masse des 
neutrinos, que du point de vue phenomenologique dans l'elucidation des anomalies 
observees dans le flux des neutrinos atmospheriques et solaires. Ces phenomenes en 
essence identiques n'ont pas regu jusqu'a present de description theorique unifiee, 
couvrant a la fois le cas des particules stables et instables, ainsi que les domaines 
relativistes et non relativistes. 

Dans le calcul traditionnel de la probabilite d'oscillation en mecanique quan- 
tique, les etats oscillants sont considered comme une superposition d'etats propres 
de masse et sont supposes etre dotes d'une energie-impulsion bien definie. Ce traite- 
ment simpliste mene directement a un paradoxe, puisque la connaissance exacte de 
l'energie-impulsion force l'etat a se trouver dans un etat propre de masse et supprime 
du meme coup les oscillations. De plus, comment une oscillation de la probabilite 
de detection de l'etat pourrait-elle etre observable si la connaissance precise de son 
energie-impulsion implique une incertitude infinie sur sa position ? 

II parait logique de resoudre ce paradoxe en postulant l'existence d'une incer- 
titude sur l'energie-impulsion, c'est-a-dire en traitant l'etat oscillant comme une 
superposition de paquets d'ondes. Cependant, des questions de principe subsistent, 
que le formalisme soit relativiste ou non : d'une part, l'interference entre etats est 
interdite en mecanique quantique non relativiste par la regie de superposition de 
Bargmann, d'autre part, il parait impossible de construire un espace de Fock pour 
une particule qui n'est pas un etat propre de masse. Une formule relativiste ne peut 
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etre derivee rigoureusement dans ces conditions. Par ailleurs, la taille et la forme du 
paquet d'ondes sont indeterminees, ce qui est insatisfaisant vu l'influence qu'elles 
ont sur l'observabilite des oscillations. L'effet des conditions de production et de 
detection n'est pas non plus pris en compte. Pour donner un exemple, une mesure 
precise de l'energie-impulsion des etats issus du lieu de detection de la particule 
oscillante, identifie l'etat propre de masse se propageant et supprime les oscillations. 
Enfin, les particules oscillantes instables ne peuvent pas etre decrites par un paquet 
d'ondes. 

II ne reste plus qu'a se tourner vers la theorie des champs, ou les particules 
oscillantes sont decrites comme des etats intermediaires virtuels, non observes di- 
rectement et se propageant entre une source et un detecteur. Dans le meme esprit, 
la question de savoir comment un etat quantique peut se signaler par une trajec- 
toire bien definie dans un detecteur, a deja ete resolue par Mott [ |107| ] en 1929, 
en ne considerant comme observables que les etats excites le long de la trajectoire 
par la particule se propageant dans le detecteur (ces etats excites se manifestent 
par exemple par des bulles dans une chambre a bulles). Pris deux a deux, les etats 
excites forment une suite de systemes source-detecteur. Contrairement a Mott qui 
utilise le theorie des perturbations non relativiste de la mecanique quantique, nous 
representons les etats oscillants par un propagateur covariant sous les transforma- 
tions de Lorentz et qui, sous sa forme complete, contient la description de l'instabilite 
des particules a travers la localisation de ses poles complexes. Un systeme de par- 
ticules en melange est represents par un propagateur matriciel non diagonal, qui 
rend possible la propagation d'une particule entre deux points avec changement de 
saveur. Ce formalisme permet d'eviter la definition d'etats de saveur, c'est-a-dire 
d'etats de masse indefmie. La matrice melangeant les etats en mecanique quantique 
est remplacee par la matrice diagonalisant le propagateur. Les corrections non ex- 
ponentielles a la propagation d'un melange peuvent etre facilement analysees dans 
ce contexte et sont negligeables. 

Dans ce formalisme, les oscillations sont mesurees indirectement, comme dans 
les experiences, par la detection des particules issues de la source et du detecteur. 
Les notions de source et de detecteur sont en fait fictives, et symbolisent les processus 
de production et de detection de la particule oscillante. Les particules entrantes et 
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sortantes sont modelisees de maniere realiste par des paquets d'ondes. Par ce biais, 
en jouant sur la largeur des paquets d'ondes, il est possible d'etudier l'influence des 
conditions de production et de detection sur l'observabilite des oscillations, puisque 
celles-ci disparaissent dans le cas ou les etats asymptotiques sont des ondes planes. 

Les problemes resultant du traitement traditionnel des oscillations n'appa- 
raissent plus dans notre calcul. lis etaient en effet lies d'une part a l'attribution 
d'une serie de proprietes (saveur, energie-impulsion, etc.) a la particule oscillante, 
et d'autre part a l'introduction de concepts classiques (temps moyen de propaga- 
tion, vitesse, etc.) par des hypotheses exterieures au formalisme. Notre methode 
met du meme coup un point final aux controverses sur la longueur d'oscillation ainsi 
que sur l'oscillation des particules associees a la production de l'etat oscillant. La 
question de l'energie-impulsion de la particule oscillante est resolue en demontrant 
en meme temps l'influence des conditions experimentales sur l'energie-impulsion et 
l'independance du resultat final (en tout cas au premier ordre en la difference de 
masse) par rapport a elle. La formule d'oscillation obtenue contient differents fac- 
teurs que l'on peut identifier comme la decroissance exponentielle en fonction de la 
distance pour une particule instable, l'oscillation dependant de la difference de masse 
et de la distance, la decoherence supprimant les oscillation a grande distance pour les 
particules relativistes et finalement l'influence des conditions experimentales sur l'ob- 
servabilite des oscillations. Les conditions sous lesquelles la formule obtenue coincide 
avec la formule classique apparaissent done explicitement dans notre resultat. 

Les caracteristiques principales de la probability d'oscillation sont ensuite repri- 
ses dans une prescription de calcul s'appliquant directement a l'amplitude d'un 
processus. Les cas des experiences CPLEAR (oscillation d'une seule particule) et 
de DA$NE sont examines. Un traitement coherent de la violation CP dans le 
systeme des kaons neutres s'ensuit, contrairement au traitement traditionnel oil 
des problemes de normalisation des etats entachent la precision des predictions 
theoriques. 

Bien que ayons considere un cas d'oscillation double (DA$NE), les hypotheses 
que nous avons posees ne permettent pas une extension immediate a d'autres situa- 
tions. II serait par exemple interessant d'etudier les oscillations en cascade B°(B°) — > 
K°(K°) ^7T7r, \11xv qui pourront etre etudiees au LHC ||108| , |109|| . Une autre applica- 
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tion interessante est l'etude de l'influence de l'instabilite de la source de la particule 
oscillante, qui peut en principe se faire en considerant aussi cette source comme un 
etat intermediaire non observe directement. Nos techniques d'integration ne con- 
viennent pas a cette analyse. Enfin, notre formule s'applique a la propagation de 
neutrinos instables, qu'ils soient relativistes ou non. Cette instabilite couplee a un 
melange est un modele explicatif des anomalies observees dans la mesure des neutri- 
nos atmospheriques et solaires. Pour l'instant, une modification simple de la formule 
de mecanique quantique par une exponentielle decroissante suffft a parametriser les 
mesures en raison du peu de donnees disponibles, mais l'astronomie au moyen des 
neutrinos ne fait que commencer ! 



Bibliographie 



[1] L. Fonda, G. C. Ghirardi et A. Rimini, Decay theory of unstable quantum systems, 
Rep. Prog. Phys. 41, 587 (1978). 

[2] L. A. Khalfin, Contribution to the decay theory of a quasi-stationary state, Sov. 
Phys. JETP 6, 1053 (1958) ; une liste des travaux pionniers concernant les corrections 
non exponentielles se trouve dans : R. G. Newton, Scattering theory of waves and 
particles (McGraw-Hill, New York, 1966), p. 608. 

[3] Voir par exemple : J. J. Sakurai, Modern quantum mechanics, rev. ed. (Addison- 
Wesley, Reading, 1994), p. 481. 

[4] J. Schwinger, Field theory of unstable particles, Ann. Phys. 9, 169 (1960). 

[5] R. Jacob et R. G. Sachs, Mass and lifetime of unstable particles, Phys. Rev. 121, 
350 (1961). 

[6] L. S. Brown, Quantum Field Theory (Cambridge University Press, Cambridge, 1992) 
p. 293 et p. 339. 

[7] P. T. Greenland, Seeking non-exponential decay, Nature 335, 298 (1988). 

[8] R. E. Peierls, Proc. of Glasgow Conference on Nuclear and Meson Physics (Pergamon 
Press, 1954) p. 296. 

[9] M. Gell-Mann et A. Pais, Behavior of neutral particles under charge conjugation, 
Phys. Rev. 97, 1387 (1955). 

[10] K. Lande, E. T. Booth, J. Impeduglia, L. M. Lederman et W. Chinowsky, Observa- 
tion of long-lived neutral V articles, Phys. Rev. 103, 1901(L) (1956). 

[11] J. H. Christenson, J. W. Cronin, V. L. Fitch et R. Turlay, Evidence for the 2n decay 
of the K% meson, Phys. Rev. Lett. 13, 138 (1964). 

[12] H. Burkhardt et al, NA31 Collaboration, First evidence for direct CP violation, Phys. 
Lett. B206, 169 (1988) ; G. D. Barr et al, NA31 Collaboration, A new measurement 
of direct CP violation in the neutral kaon system, Phys. Lett. B317, 233 (1993) ; L. 
K. Gibbons et al, E731 Collaboration, Measurement of the CP-violation parameter 
Re{e'/e), Phys. Rev. Lett. 70, 1203 (1993) ; V. Fanti et al, NA48 Collaboration, A 
new measurement of direct CP violation in two pion decays of the neutral kaon, Phys. 
Lett. B465, 335 (1999) ; KTeV Collaboration, Observation of Direct CP Violation 
in K s ,l 7T7T Decays, Phys. Rev. Lett. 83, 22 (1999). 



147 



148 



BIBLIO G RAP HIE 



[13] E. Wigner et V. F. Weisskopf, Berechnung der natiirlichen Linienbreite auf Grund 
dei Diracschen Lichttheorie, Z. Phys. 63, 54 (1930) ; ibid., Uber die natiirliche Lin- 
ienbreite in der Strahlung des harmonischen Oszillators, Z. Phys. 65, 18 (1930) ; 
pour une presentation moderne, voir par exemple la ref. p. 509. 

[14] N. R. Lipshutz, Invariance principles and the K°—K° propagator matrix, Phys. Rev. 
144, 1300 (1966). 

[15] H. Albrecht et al, ARGUS Collaboration, Observation of B°-W mixing, Phys. Lett. 
B192, 245 (1987) ; T. Affolder et al, CDF Collaboration, Measurement of the B°W 
oscillation frequency using l~D* + pairs and lepton flavor tags, Phys. Rev. D60, 
112004 (1999). 

[16] T. Affolder et al, CDF Collaboration, A measurement of sin 2/? from B — ► J/ipK^ 
with the CDF detector, Phys. Rev. D61, 072005 (2000). 

[17] T. B. Day, Demonstration of quantum mechanics in the large, Phys. Rev. 121, 1204 
(1961) ; H. Lipkin, CP violation and coherent decays of kaon pairs, Phys. Rev. 176, 
1715 (1968). 

[18] R. Davis, D. S. Harmer et K. C. Hoffman, Search for neutrinos from the sun, Phys. 
Rev. Lett. 20, 1205 (1968). 

[19] J. N. Bahcall, Neutrino astrophysics (Cambridge University Press, Cambridge, 
1989). 

[20] S. M. Bilenky, C. Giunti et W. Grimus, Phenomenology of neutrino oscillations, 
Prog. Part. Nucl. Phys. 43, 1 (1999) et references incluses. 

[21] B. Pontecorvo, Neutrino experiments and the problem of conservation of leptonic 
charge, Sov. Phys. JETP 26, 984 (1968) ; V. Gribov et B. Pontecorvo, Neutrino 
astronomy and lepton charge, Phys. Lett. 28B, 493 (1969). 

[22] J. N. Bahcall, N. Cabibbo et A. Yahil, Are neutrinos stable particles ?, Phys. Rev. 
Lett. 28, 316 (1972) ; J. N. Bahcall, S. T. Petcov, S. Toshev et J. W. F. Valle, Tests 
of neutrino stability, Phys. Lett. B181, 369 (1986). 

[23] Y. Fukuda et al, Super-Kamiokande Collaboration, Evidence for oscillation of atmo- 
spheric neutrinos, Phys. Rev. Lett. 81, 1562 (1998) ; K. S. Hirata et al, Kamiokande- 
II Collaboration, Experimental study of the atmospheric neutrino flux, Phys. Lett. 
B205, 416 (1988). 

[24] J. G. Learned, S. Pakvasa et T. J. Weiler, Neutrino mass and mixing implied by un- 
derground dehcit of low energy muon-neutrino events, Phys. Lett. B207, 79 (1988) ; 
V. Barger et K. Whisnant, The effects of neutrino oscillations with one mass scale on 
the atmospheric neutrino flux, Phys. Lett. B209, 365 (1988) ; K. Hidaka, M. Honda 
et S. Midorikawa, Neutrino oscillations and the anomalous atmospheric neutrino 
flux, Phys. Rev. Lett. 61, 1537 (1988). 

[25] V. Barger, J. G. Learned, S. Pakvasa et T. J. Weiler, Neutrino Decay as an Expla- 
nation of Atmospheric Neutrino Observations, Phys. Rev. Lett. 82, 2640 (1999) ; V. 
Barger, J. G. Learned, P. Lipari, M. Lusignoli, S. Pakvasa et T. J. Weiler, Neutrino 
Decay and Atmospheric Neutrinos, Phys. Lett. B462, 109 (1999). 



BIBLIO G RAP HIE 



149 



C. Athanassopoulos et al, LSND Collaboration, Evidence for — ► v e oscillations 
from the LSND experiment at LAMPF, Phys. Rev. Lett. 77, 3082 (1996) ; C. 
Athanassopoulos et al, LSND Collaboration, Results on — ► u e neutrino oscilla- 
tions from the LSND experiment, Phys. Rev. Lett. 81, 1774 (1998). 

L. Wolfenstein, Neutrino oscillations in matter, Phys. Rev. D17, 2369 (1978) ; S. 
P. Mikheyev et A. Yu. Smirnov, Resonance enhancement of oscillations in matter 
and solar neutrino spectroscopy, Sov. J. Nucl. Phys. 42, 913 (1985) ; H. A. Bethe, 
Possible explanation of the solar-neutrino puzzle, Phys. Rev. Lett. 56, 1305 (1986). 

S. M. Bilenky et C. Giunti, See-saw type mixing and v„ —* v T oscillations, Phys. Lett. 
B300, 137 (1993). 

Z. Maki, M. Nakagawa et S. Sakata, Remarks on the unified model of elementary 
particles, Prog. Theor. Phys. 28, 870 (1962). 

C. Giunti, C. W. Kim, U. W. Lee, Remarks on the weak states of neutrinos, Phys. 
Rev. D45, 2414 (1992). 

O. Nachtmann, Elementary particle physics (Springer- Verlag, Berlin, 1990) p. 443. 

R. G. Sachs, Interference phenomena of neutral K mesons, Ann. Phys. 22, 239 
(1963). 

C. P. Enz et R. R. Lewis, On the phenomenological description of CP violation for 
K-mesons and its consequences, Helv. Phys. Acta 38, 860 (1965). 

L. Alvarez-Gaume, C. Kounnas, S. Lola et P. Pavlopoulos, Violation of time-reversal 
invariance and CPLEAR measurements, Phys. Lett. B458, 347 (1999). 

C. Jarlskog, Commutators of the quark mass matrices in the standard electroweak 
model and a measure of maximal CP violation, Phys. Rev. Lett. 55, 1039 (1985). 

K. Gottfried et V. F. Weisskopf, Concepts of Particle Physics (Clarendon Press, 
Oxford, 1984) vol. I p. 151. 

Y. Srivastava, A. Widom et E. Sassaroli, Spatial correlations in two neutral kaon 
decays, Z. Phys. C66, 601 (1995). 

H. J. Lipkin, Theories of non-experiments in coherent decays of neutral mesons, 
Phys. Lett. B348, 604 (1995). 

J. Lowe, B. Bassaleck, H. Burkhardt, A. Rusek, G. J. Stephenson Jr. et T. Goldman, 
No A oscillations, Phys. Lett. B384, 288 (1996) ; H. Burkhardt, J. Lowe, G. J. 
Stephenson Jr. et T. Goldman, Oscillations of recoil particles against mixed states, 
Phys. Rev. D59, 054018 (1999). 

B. Kayser, The frequency of neutral meson and neutrino oscillation, SLAC-PUB- 
7123 ; B. Kayser, CP violation, mixing, and quantum mechanics, Proc. of 28th HEP 



Conf., Varsovie 1996, p. 1135 ( |hep-ph/9702327|) 



Voir par exemple : R. N. Mohapatra et P. B. Pal, Massive neutrinos in physics and 
astrophysics, (World Scientific, Singapore, 1991), p. 156. 

R. G. Winter, Neutrino oscillation kinematics, Lett. Nuovo Cimento 30, 101 (1981). 



T. Goldman, Source dependence of neutrino oscillations, hep-ph / 9604357 . 



150 



BIBLIO G RAP HIE 



[44] A. D. Dolgov, A. Yu. Morozov, L. B. Okun et M. G. Schepkin, Do muons oscillate ?, 
Nucl. Phys. B502, 3 (1997). 

[45] H. J. Lipkin, Quantum mechanics of neutrino oscillations - Hand waving for pedes- 
trians, |hep-ph/990i399| . 

[46] V. Bargmann, On unitary ray representations of continuous groups, Ann. Math. 59, 
1 (1954) ; A. Galindo et P. Pascual, Quantum Mechanics, Springer- Verlag, 1990, p. 
292; F. A. Kaempffer, Concepts in quantum mechanics (Academic Press, 1965), p. 
341 ; cette derniere reference est aussi interessante pour sa presentation des theories 
de Yang-Mills et de la gravitation comme theorie de jauge. 

[47] M. Blasone, P. A. Henning et G. Vitiello, The exact formula for neutrino oscillations, 
Phys. Lett. B451, 140 (1999) ; M. Blasone et G. Vitiello, Remarks on the neutrino 
oscillation formula, Phys. Rev. D60, 111302 (1999). 

[48] M. Beuthe, G. Lopez Castro et J. Pestieau, Field theory approach to K° — K° and 
B°-W systems, Int. J. Mod. Phys. A13, 3587 (1998). 

[49] M. Veltman, Unitarity and causality in a renormalizable field theory with unsta- 
ble particles, Physica 29, 186 (1963) ; ibid., Diagrammatica (Cambridge University 
Press, Cambridge, 1994). 

[50] Y. Srivastava, A. Widom et E. Sassaroli, A oscillations, Phys. Lett. B344, 436 
(1995). 

[51] Y. Srivastava, A. Widom et E. Sassaroli, Charged lepton oscillations, |hep-| 
ph/9509261 ; ibid., Charged leptons and neutrino oscillations, Eur. Phys. J. C2, 769 



(1998) ; Y. Srivastava et A. Widom, Of course muons can oscillate, |hep-ph/9707268 



[52] B. Kayser, On the quantum mechanics of neutrino oscillation, Phys. Rev. D24, 110 
(1981). 

[53] C. Giunti, C. W. Kim, J. A. Lee et U. W. Lee, Treatment of neutrino oscillations 
without resort to weak eigenstates, Phys. Rev. D48, 4310 (1993) ; C. Giunti, C. W. 
Kim et U. W. Lee, When do neutrinos cease to oscillate ?, Phys. Lett. B421, 237 
(1998). 

[54] R. G. Stuart, Gauge invariance, analyticity and physical observables at the Z° reso- 
nance, Phys. Lett. B262, 113 (1991) ; U. Baur et D. Zeppenfeld, Finite width effects 
and gauge invariance in radiative W production and decay, Phys. Rev. Lett. 75, 
1002 (1995) ; J. Papavassiliou et A. Pilaftsis, Gauge-independent approach to reso- 
nant transition amplitudes, Phys. Rev. D53, 2128 (1996) ; M. Beuthe, R. Gonzalez 
Felipe, G. Lopez Castro et J. Pestieau, Behaviour of the absorptive part of the 
electromagnetic vertex, Nucl. Phys. B498, 55 (1997) et references incluses. 

[55] G. Kallen, On the definitions of the renormalization constants, Helv. Phys. Acta 
25, 417 (1957) ; H. Lehmann, On the properties of propagation functions and renor- 
malization constants of quantized fields, Nuovo Cimento 11, 342 (1954) ; pour une 
presentation moderne, voir par exemple H, p. 282. 

[56] R. E. Cutkosky, Singularities and discontinuities of Feynman amplitudes, J. Math. 
Phys. 1, 429 (1960) ; voir aussi ||, p. 235. 



BIBLIO G RAP HIE 



151 



[57] G. Buchalla, A. J. Buras et M. E. Lautenbacher, Weak decays beyond leading loga- 
rithms, Rev. Mod. Phys. 68, 1125 (1996). 

L. Okun, Leptons and quarks (North-Holland, Amsterdam, 1982), p. 85. 

L. Baulieu et R. Coquereaux, Photon-Z mixing in the Weinberg-Salam model, Ann. 
Phys. 140, 163 (1982). 

S. Coleman et H. J. Schnitzer, Mixing of elementary particles, Phys. Rev. 134, B863 
(1964) ; H. B. O'Connell, Recent developments in rho-omega mixing, Aust. J. Phys. 
50, 255 (1997), [hep-ph/ 9604375) . 

J. Rich, Quantum mechanics of neutrino oscillations, Phys. Rev. D48, 4318 (1993). 

W. Grimus, P. Stockinger et S. Mohanty, The field-theoretical approach to coherence 
in neutrino oscillations, Phys. Rev. D59, 013011 (1999). 

J. E. Campagne, Neutrino oscillations from apion decay in flight, Phys. Lett. B400, 
135 (1997). 

W. Grimus et P. Stockinger, Real oscillations of virtual neutrinos, Phys. Rev. D54, 
3414 (1996). 

M. E. Peskin et D. V. Schroeder, An introduction to quantum field theory (Addison- 
Wesley, Reading, 1995) p. 102. 

D. Sudarsky, E. Fischbach, C. Talmadge, S. H. Aronson et H.-Y. Cheng, Effects of 
external fields on the neutral kaon system, Ann. Phys. 207, 103 (1991). 

A. Peres, Nonexponential decay law, Ann. Phys. 129, 33 (1980). 

J. R. Fry, CP violation and the standard model, Rep. Prog. Phys. 63, 117 (2000). 

C. Caso et al, Particle Data Group, .Review of Particle Physics, Eur. Phys. J. C3, 
1 (1998) ; la nouvelle edition vient de paraitre : D. E. Groom et al, Particle Data 
Group, Review of Particle Physics, Eur. Phys. J. C15, 1 (2000) ; on peut aussi 



consulter le site http ://pdg.lbl.gov/ 



R. Fleischer, CP Violation and the role of electroweak penguins in nonleptonic B 
decays, Int. J. Mod. Phys. A12, 2459 (1997), et references incluses. 

C. B. Chiu et E. C. G. Sudarshan, Decay and evolution of the neutral kaon, Phys. 
Rev. D42, 3712 (1990). 

Q. Wang et A. I. Sanda, Neutral kaon system reinvestigated, Phys. Rev. D55, 3131 
(1997). 

C. Giunti, C. W. Kim et U. W. Lee, When do neutrinos really oscillate ? Quantum 
mechanics of neutrino oscillations, Phys. Rev. D44, 3635 (1991) ; C. W. Kim et 
A. Pevsner, Neutrinos in physics and astrophysics (Harwood Academic Publishers, 
Chur, 1993), chap. 9 ; C. Giunti et C. W. Kim, Coherence of neutrino oscillations in 
the wave packet approach, Phys. Rev. D58, 017301 (1998). 

[74] S. Nussinov, Solar neutrinos and neutrino mixing, Phys. Rev. 63B, 201 (1976) ; S. 
M. Bilenky et B. Pontecorvo, Lepton mixing and neutrino oscillation, Phys. Rep. 
41C, 225 (1978). 



152 BIBLIO G RAP HIE 



[75] K. Kiers, S. Nussinov et N. Weiss, Coherence effects in neutrino oscillations, Phys. 
Rev. D53, 537 (1996). 

[76] J. N. Bahcall et S. C. Frautschi, Lepton non-conservation and solar neutrinos, Phys. 
Lett. 29B, 623 (1969). 

[77] S. L. Glashow et L. M. Krauss, "Just so" neutrino oscillations, Phys. Lett. B190, 
199 (1987). 

[78] L. Krauss et F. Wilczek, Solar-neutrino oscillations, Phys. Rev. Lett. 55, 122 (1985). 

[79] D. V. Ahluwalia et T. Goldman, Interplay of non-relativistic and relativistic effects 
in neutrinos oscillations, Phys. Rev. D56, 1698 (1997). 

[80] P. Reinartz et L. Stodolsky, Neutrino masses and mixings in supernova bursts, Z. 
Phys. C27, 507 (1985). 

[81] J. Schechter et J. W. F. Valle, Neutrinos masses in SU(2) x 17(1) theories, Phys. 
Rev. D22, 2227 (1980). 

[82] B. Kayser, F. Gibrat-Debu et F. Perrier, The physics of massive neutrinos (World 
Scientific, Singapore, 1989), p. 59. 

[83] I. Yu. Kobzarev, B. V. Martem'yanov, L. B. Okun' et M. G. Shchepkin, Phenomenol- 
ogy of neutrino oscillations, Sov. J. Nucl. Phys. 32, 823 (1980) ; S. M. Bilenky, J. 
Hosek et S. T. Petcov, On the oscillations of neutrinos with Dirac and Majorana 
masses, Phys. Lett. 94B, 495 (1980). 

[84] C. Y. Cardall, Coherence of neutrino flavor mixing in quantum field theory, Phys. 
Rev. D61, 073006 (2000) ; C. Y. Cardall et D. J. H. Chung, The MSW effect in 
quantum field theory, Phys. Rev. D60, 073012 (1999). 

[85] Yu. V. Shtanov, Space-time description of neutrino flavour oscillations, Phys. Rev. 
D57, 4418 (1998). 

[86] K. Kiers et N. Weiss, Neutrino oscillations in a model with a source and detector, 
Phys. Rev. D57, 3091 (1998). 

[87] A. Ioannisian et A. Pilaftsis, Neutrino oscillations in space within a solvable model, 
Phys. Rev. D59, 053003 (1999). 

[88] The DA&NE Physics Handbook Vol. I, eds. L. Maiani, G. Pancheri et N. Paver 
(INFN-LNF Publications, 1992) ; P. H. Eberhard, in : The Second DA&NE Physics 
Handbook, Vol. I, p. 99, eds. L. Maiani, G. Pancheri et N. Paver (INFN-LNF Pub- 
lications, 1995). 

[89] H. Albrecht et al, Collaboration ARGUS, A study of 5° -> D*+l~i> and B°W 
mixing using partial D* + reconstruction, Phys. Lett. B324, 249 (1994) ; J. Bartelt 
et al, Collaboration CLEO, Two measurements of B°B° mixing, Phys. Rev. Lett. 
71, 1680 (1993) ; H. Schroder, BB mixing, in : B decays, Rev. 2nd ed., ed. S. Stone 
(World Scientific, Singapore, 1994) p. 449. 

[90] KEK-B 1995 B factory Design Report KEK; BABAR 1995 Technical Design Report 
SLAC-R-95-0457. 



BIBLIO G RAP HIE 



153 



[91] I. Dunietz, J. Hauser et J. Rosner, Proposed experiment addressing CP and CPT 
violation in the K -K® system, Phys. Rev. D35, 2166 (1987) ; C. D. Buchanan, R. 
Cousins, D. Dib, R. D. Peccei et J. Quackenbush, Testing CP and CPT violation in 
the neutral kaon system at a factory, Phys. Rev. D45, 4088 (1992) ; M. Hayakawa 
et A. I. Sanda, Searching for T, CP, CPT, and AS = AQ rule violations in the 
neutral K system : A guide, Phys. Rev. D48, 1150 (1993). 

[92] M. Zrajek, From kaons to neutrinos : quantum mechanics of particle oscillations, 
Acta Phys. Polon. B29, 3925 (1998). 

[93] A. Angelopoulos et al, CPLEAR Collaboration, Measurement of the K^ — Ks mass 
difference using semileptonic decays of tagged neutral kaons, Phys. Lett. B444, 38 
(1998). 

[94] A. Angelopoulos et al, CPLEAR Collaboration, First direct observation of time- 
reversal non-invariance in the neutral-kaon system, Phys. Lett. B444, 43 (1998). 

[95] R. Adler et al, CPLEAR Collaboration, First determination of CP violation param- 
eters from K -K® decay asymmetry, Phys. Lett. B286, 180 (1992) ; A. Apostolakis 
et al, CPLEAR Collaboration, A determination of the CP violation parameter rj + _ 
from the decay of strangeness-tagged neutral kaons, Phys. Lett. B458, 545 (1999). 

[96] R. Adler et al, CPLEAR Collaboration, First observation of a particle-antiparticle 
asymmetry in the decay of neutral kaons into tt°tt , Z. Phys. C70, 211 (1996) ; 
A. Angelopoulos et al, CPLEAR Collaboration, Measurement of the CP-violation 
parameter rj 00 using tagged K° and K°, Phys. Lett. B420, 191 (1998). 

[97] A. Datta et D. Kumbhakar, D°-Pfi mixing : a possible test of physics beyond the 
standard model, Z. Phys. C27, 515 (1985). 

[98] A. B. Carter et A. I. Sanda, CP nonconservation in cascade decays of B mesons, 
Phys. Rev. Lett. 45, 952 (1980) ; I. I. Bigi et A. I. Sanda. Notes on the observability 
of CP violations in B decays, Nucl. Phys. B193, 85 (1981) ; ibid., CP violation in 
heavy flavor decays : predictions and search strategies, B281, 41 (1987) ; voir aussi 
la reference Q . 

[99] D. Delepine, J.-M. Gerard, J. Pestieau et J. Weyers, Final state interaction phases 
in (B^Ktt) decay amplitudes, Phys. Lett. B429, 106 (1998). 

[100] Y. Nir et H. R. Quinn, Theory of CP violation in B decays, in : B decays, Rev. 2nd 
ed., ed. S. Stone (World Scientific, Singapore, 1994) p. 520 ; une synthese figure dans 
la reference @, p. 555 a 562. 

[101] L.-L. Chau, Quark mixing in weak interactions, Phys. Rep. 95, 1 (1983). 

[102] W. Grimus, CP violating phenomena and theoretical results, Fortschr. Phys. 36, 
201 (1988) et references incluses. 

[103] L. Wolfenstein, Models of CP violation, in : Theory and phenomenology in particle 
physics, ed. A. Zichichi (Academic Press, New York, 1969), Part A, p. 218. 

[104] J.-M. Gerard, CP- and T- violations in the standard model in : Cargese, 1992 : 
Quantitative particle physics, p. 149. 



154 



BIBLIO G RAP HIE 



[105] R. E. Marshak, Riazuddin et C. P. Ryan, Theory of weak interactions in particle 
physics, (Wiley-Interscience, New York, 1969), p. 546. 

[106] A. J. Buras et J. M. Gerard, Isospin breaking contributions to e'/e, Phys. Lett. 
B192, 156 (1987). 

[107] N. F. Mott, The wave mechanics of a-ray tracks, Proc. R. Soc. London, A126, 
79 (1929), reprinted in Quantum theory and measurement, ed. J. A. Wheeler et 
W. H. Zurek (Princeton University Press, Princeton, 1983) ; L. I. Schiff, Quantum 
mechanics, 3rd ed. (McGraw-Hill Kogakusha, Tokyo, 1968), p. 335. 

[108] Ya. I. Azimov, K° decays as analyzers of B° decays : How to measure the sign 
of Amg, Phys. Rev. D42, 3705 (1990) ; ibid., Phenomenology of neutral D-meson 
decays and double flavor oscillations, Eur. Phys. J. A4, 21 (1999). 

[109] B. Kayser et L. Stodolsky, Cascade mixing, a new kind of particle mixing phe- 
nomenon, hep-ph/9610522 ; B. Kayser, Cascade mixing and the CP-violating angle 



beta, |hep-ph/9709"382 



